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%	85.6 / NA	27	pseudogene (769/2001 nt)	yehQ	pseudogene
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/17);  total (10/17)

GACGATCCGTTTGCCAGCGAAGAAGGCAACGCGTGTCGTCAATATGTTCAGCAATTAAGCCAGACATTATGGCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGACCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTA  >  NC_000913/2209727‑2209999                                                                                                                                              |                                                                                                                                    gACGATCCGTTTGCCAGCGAAGAAGGCAACGCGTGTCGTCAATATGTTCAGCAATTAAGCCAGACATTATGGCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGc                                                                                                                                <  1:672770/147‑1 (MQ=255) gACGATCCGTTTGCCAGCGAAGAAGGCAACGCGTGTCGTCAATATGTTCAGCAATTAAGCCAGACATTATGGCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCaa                                                                                                                              <  2:182310/149‑1 (MQ=255)         gTTTGCCAGCGAAGAAGGCAACGCGTGTCGTCAATATGTTCAGCAATTAAGCCAGACATTATGGCTTGGCGGCATCAGCCAGCCGCCCATCCCTTACGAGGCAGCATTCAACCGCGCATTGCCGGCGGCCGAGGCCTGCAACTGGCGCTgg                                                                                                                    >  1:176467/1‑151 (MQ=255)                     agaagGCAACGCGTGTCGTCAATATGTTCAGCAATTAAGCCAGACATTATGGCTTGGCGGCATCAGCCAGCCCCTCATCCATAACGAGGCAGACTTCAACCGCGCTTTGCAGGAGGCCGCGGCCCGCCACCGGCGCTGGGGGGGTGaa                                                                                                           >  1:538754/1‑148 (MQ=255)                          gCAACGCGTGTCGTCAATATGTTCAGCAATTAAGCCAGACATTATGGCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACg                                                                                                   <  2:754278/151‑1 (MQ=255)                                                       ttAAGCCAGACATTATGGCTTGGCGGCATCAGCCAGCCGCTCCTCCATTACGAGGCAGCATTTAACCGCGCCTTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGGGTGAATCGCTACGGCCACTGCGCGCCAGCGTTGATGCCTTcc                                                                      >  1:453486/1‑151 (MQ=255)                                                         aaGCCAGACATTATGGCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCAc                                                                    <  1:140196/151‑1 (MQ=255)                                                         aaGCCAGACATTATGGCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCAc                                                                    <  1:330322/151‑1 (MQ=255)                                                               gACATTATGGCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCcgc                                                              <  2:575944/151‑1 (MQ=255)                                                                      tGGCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCgccagcc                                                       >  2:246347/1‑151 (MQ=255)                                                                       gtcTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCgccagcca                                                      <  1:494148/149‑1 (MQ=255)                                                                            ggcggcATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACtata                                                 <  1:491974/151‑1 (MQ=255)                                                                            ggcggcATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACtat                                                  <  1:341526/150‑1 (MQ=255)                                                                                 cATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTCCGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGc                                             <  1:161256/150‑1 (MQ=255)                                                                                                   ccATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATc                            >  2:518876/1‑149 (MQ=255)                                                                                                   ccATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAg                          <  2:342968/151‑1 (MQ=255)                                                                                                   ccATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGAGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCa                           <  1:315219/150‑1 (MQ=255)                                                                                                    cATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGc                         <  2:581828/151‑1 (MQ=255)                                                                                                             gcagcaTTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCAt                <  1:643556/151‑1 (MQ=255)                                                                                                              cagcaTTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGTTGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCAtt               >  2:340829/1‑151 (MQ=255)                                                                                                                gcaTTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCCCGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTa              >  1:677741/1‑150 (MQ=255)                                                                                                                 caTTCAACCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTaa             <  1:768028/150‑1 (MQ=255)                                                                                                                      aaCCGCGCATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTc       >  2:510585/1‑151 (MQ=255)                                                                                                                          gcgcATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGAtt   <  1:65366/151‑1 (MQ=255)                                                                                                                           cgcATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGAtt   >  1:64370/1‑150 (MQ=255)                                                                                                                           cgcATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTa  <  1:674529/151‑1 (MQ=255)                                                                                                                            gcATTGCAGGCGGCAGAGGCCTGCAACTGGCGCTGGGGGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTa  >  1:531972/1‑150 (MQ=255)                                                                                                                                              |                                                                                                                                    GACGATCCGTTTGCCAGCGAAGAAGGCAACGCGTGTCGTCAATATGTTCAGCAATTAAGCCAGACATTATGGCTTGGCGGCATCAGCCAGCCGCTCATCCATTACGAGGCAGCATTCAACCGCGCATTGCAGGCGGCAGAGACCTGCAACTGGCGCTGGGTGAGTGAATCGCTACGGCAACTGCGCGCCAGCGTTGATGCCTTCCACGCCCGCGCCAGCCACTATAATGCCGGAGAATGCTTACATCAGCTTGCGGCATTAAACAGTCGATTA  >  NC_000913/2209727‑2209999

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