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

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

ATCATCTGACGATCCGTTTGCCAGCGAAGAAGGCAACGCGTGTCGTCAATATGTTCAGCAATTAAGCCAGACATTATGGCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGACCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTAAATTGCGCACAAGAGA  >  NC_000913/2209720‑2210015                                                                                                                                                     |                                                                                                                                                    atcatcTGACGATCCGTTTGCCAGCGAAGAAGGCAACGCGTGTCGTCAATATGTTCAGCAATTAAGCCAGACATTATGGCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGcc                                                                                                                                                   <  2:276042/151‑1 (MQ=255)      cTGACGATCCGTTTGCCAGCGAAGAAGGCAACGCGTGTCGTCAATATGTTCAGCAATTAAGCCAGACATTATGGCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCa                                                                                                                                               >  2:349541/1‑150 (MQ=255)        gACGATCCGTTTGCCAGCGAAGAAGGCAACGCGTGTCGTCAATATGTTCAGCAATTAAGCCAGACATTATGGCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGc                                                                                                                                                <  1:303535/147‑1 (MQ=255)        gACGATCCGTTTGCCAGCGAAGAAGGCAACGCGTGTCGTCAATATGTTCAGCAATTAAGCCAGACATTATGGCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAAACGCGCATTGCAGGCGGCAGAGGCCCTCAAc                                                                                                                                             >  1:887427/1‑150 (MQ=255)           gATCCGTTTGCCAGCGAAGAAGGCAACGCGTGTCGTCAATATGTTCAGCAATTAAGCCAGACATTATGGCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTgg                                                                                                                                          >  1:131379/1‑150 (MQ=255)               cGTTTGCCAGCGAAGAAGGCAACGCGTGTCGTCAATATGTTCAGCAATTAAGCCAGACATTATGGCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTg                                                                                                                                     <  1:168972/151‑1 (MQ=255)                       aGCGAAGAAGGCAACGCGTGTCGTCAATATGTTCAGCAATTAAGCCAGACATTATGGCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGgtgagtg                                                                                                                             <  1:310803/151‑1 (MQ=255)                              aagGCAACGCGTGTCGTCAATATGTTCAGCAATTAAGCCAGACATTATGGCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGGGAGTGAATCGc                                                                                                                       >  1:347893/1‑150 (MQ=255)                                  cAACGCGTGTCGTCAATATGTTCAGCAATTAAGCCAGACATTATGGCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACgg                                                                                                                  >  1:466359/1‑151 (MQ=255)                                  cAACGCGTGTCGTCAATATGTTCAGCAATTAAGCCAGACATTATGGCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACgg                                                                                                                  <  1:509161/151‑1 (MQ=255)                                   aaCGCGTGTCGTCAATATGTTCAGCAATTAAGCCAGACATTATGGCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACgg                                                                                                                  >  2:765830/1‑150 (MQ=255)                                    aCGCGTGTCGTCAATATGTTCAGCAATTAAGCCAGACATTATGGCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGGGAGTGAATCGCTACGGc                                                                                                                 >  1:871101/1‑150 (MQ=255)                                         tgtCGTCAATATGTTCAGCAATTAAGCCAGACATTATGGCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGGGAGTGAATCGCTACGGCAAc                                                                                                              >  2:673488/1‑148 (MQ=255)                                          gtcgtcAATATGTTCAGCAATTAAGCCAGACATTATGGCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTgcg                                                                                                          >  1:379633/1‑151 (MQ=255)                                                aaTATGTTCAGCAATTAAGCCAGACATTATGGCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGcc                                                                                                      >  1:886991/1‑149 (MQ=255)                                                    tGTTCAGCAATTAAGCCAGACATTATGGCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGtt                                                                                                >  1:380173/1‑151 (MQ=255)                                                     gTTCAGCAATTAAGCCAGACATTATGGCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTg                                                                                               <  1:594112/151‑1 (MQ=255)                                                                     aGACATTATGGCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCCACTGCGCGCCAGCGTTGATGCCTTCCACGccc                                                                                >  2:753639/1‑150 (MQ=255)                                                                      gACATTATGGCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCcgc                                                                              <  2:106301/151‑1 (MQ=255)                                                                              ggCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGACGCCTGCAACTGGCGCTGGGCGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCCCGCCCGCgccaac                                                                        >  2:711482/1‑147 (MQ=255)                                                                                     cggcATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAAt                                                               >  1:472123/1‑151 (MQ=255)                                                                                       gcATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGc                                                             <  2:60703/151‑1 (MQ=255)                                                                                       gcATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGc                                                             <  1:243726/151‑1 (MQ=255)                                                                                         aTCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCg                                                           <  2:114898/151‑1 (MQ=255)                                                                                           cagccagcCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCgg                                                          <  2:95530/150‑1 (MQ=255)                                                                                           cagccagcCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGga                                                         <  1:610524/151‑1 (MQ=255)                                                                                              ccagccGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGaa                                                      <  2:887356/151‑1 (MQ=255)                                                                                               cagccGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAAt                                                     <  1:655307/151‑1 (MQ=255)                                                                                                  ccGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGc                                                   <  1:662461/150‑1 (MQ=255)                                                                                                  ccGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCt                                                  >  1:705266/1‑151 (MQ=255)                                                                                                  ccGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCt                                                  >  1:651647/1‑151 (MQ=255)                                                                                                      tCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTAc                                               >  1:596396/1‑150 (MQ=255)                                                                                                         tCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGACGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATc                                            >  2:254019/1‑150 (MQ=255)                                                                                                          ccATTACGAGGCAGCATTCACCCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCa                                           <  1:322782/150‑1 (MQ=255)                                                                                                          ccATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAg                                          <  1:861939/151‑1 (MQ=255)                                                                                                           cATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGc                                         <  2:407272/151‑1 (MQ=255)                                                                                                                cGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGCTTGCg                                    <  1:409739/151‑1 (MQ=255)                                                                                                                 gAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGGGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTACCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGCtgg                                      <  2:399019/148‑3 (MQ=255)                                                                                                                  aGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGCTTGCGGc                                  >  2:887993/1‑151 (MQ=255)                                                                                                                    gcagcaTTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCAt                                >  2:607748/1‑151 (MQ=255)                                                                                                                     cagcaTTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCAt                                >  2:240279/1‑150 (MQ=255)                                                                                                                        caTTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGCGAATGCTTACATCAGCCTGCGGCATTaa                             >  2:228829/1‑150 (MQ=255)                                                                                                                        caTTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTaaa                            <  1:621347/151‑1 (MQ=255)                                                                                                                             aaCCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTc                       >  1:291352/1‑151 (MQ=255)                                                                                                                               ccGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTc                       <  2:705217/149‑1 (MQ=255)                                                                                                                                 gcgcATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGAtt                   <  1:853317/151‑1 (MQ=255)                                                                                                                                           ggcggcAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTAAATTGCGca         <  2:781407/151‑1 (MQ=255)                                                                                                                                              ggcAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTAAATTGCGCACaa      <  1:571524/151‑1 (MQ=255)                                                                                                                                               gcAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTAAATTGCGCACAag     <  2:851232/151‑1 (MQ=255)                                                                                                                                                 agagGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTAAATTGCGCACAaga    <  1:383428/150‑1 (MQ=255)                                                                                                                                                  gaggCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGTTTGCGGCATTAAACAGTCGATTAAATTGCGCACAAgaga  <  1:554524/151‑1 (MQ=255)                                                                                                                                                  gaggCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTAAATTGCGCACAag     >  2:857149/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                    ATCATCTGACGATCCGTTTGCCAGCGAAGAAGGCAACGCGTGTCGTCAATATGTTCAGCAATTAAGCCAGACATTATGGCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGACCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTAAATTGCGCACAAGAGA  >  NC_000913/2209720‑2210015

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