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%	159.5 / NA	42	pseudogene (177/2001 nt)	yehQ	pseudogene
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/23);  total (19/23)

TGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTGGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACTTTGTGCCACCACTCAGTCTACGGAAAAAGAAGAAGAGTGGGATCCG  >  NC_000913/2209131‑2209417                                                                                                                                                  |                                                                                                                                              tGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGa                                                                                                                                           >  1:297284/1‑150 (MQ=255)    cacCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACgg                                                                                                                                       >  1:250521/1‑151 (MQ=255)         cAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAgg                                                                                                                                  >  1:293633/1‑151 (MQ=255)          aGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGGCGGTGTCCGTACCCAGCTAGCGAACGGCCAgg                                                                                                                                  >  2:644154/1‑150 (MQ=255)            gCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCAc                                                                                                                               >  1:617957/1‑151 (MQ=255)                 gACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGaaa                                                                                                                          <  2:297317/151‑1 (MQ=255)                 gACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGaa                                                                                                                           <  1:652673/150‑1 (MQ=255)                     gCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAgg                                                                                                                      >  2:402297/1‑151 (MQ=255)                     gCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAgg                                                                                                                      <  2:140882/151‑1 (MQ=255)                      cGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGc                                                                                                                     >  1:559205/1‑151 (MQ=255)                       gTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGc                                                                                                                     >  1:82317/1‑150 (MQ=255)                       gTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCt                                                                                                                    <  1:402248/151‑1 (MQ=255)                          aaGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAg                                                                                                                 <  1:632882/151‑1 (MQ=255)                                tGCCGGTTTTGTTAAGCGCAGTCGTAAGGAACTGGAAAATGGCAACGTTCCGGAGATCAGCCATGAGAAGGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCa                                                                                                             <  1:67316/149‑1 (MQ=255)                                 gCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGtt                                                                                                          >  1:515911/1‑151 (MQ=255)                                     gTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCgg                                                                                                      >  1:182277/1‑151 (MQ=255)                                       tttGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCgggg                                                                                                    >  2:494200/1‑151 (MQ=255)                                                 gcAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTAt                                                                                          >  1:358858/1‑151 (MQ=255)                                                 gcAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTAt                                                                                          >  1:745752/1‑151 (MQ=255)                                                        tAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGt                                                                                   >  2:257555/1‑151 (MQ=255)                                                          aGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTc                                                                                  >  2:216540/1‑150 (MQ=255)                                                          aGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCa                                                                                 <  2:571666/151‑1 (MQ=255)                                                                 ggAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATcgcg                                                                            <  2:218044/149‑1 (MQ=255)                                                                 ggAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTg                                                                          <  1:378389/151‑1 (MQ=255)                                                                        ggCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTg                                                                    <  1:288085/150‑1 (MQ=255)                                                                         gCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGgt                                                                  <  1:715515/151‑1 (MQ=255)                                                                           aaCGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGgtgt                                                                <  2:605242/151‑1 (MQ=255)                                                                                cccGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAg                                                            <  1:15285/150‑1 (MQ=255)                                                                                   ggAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTAt                                                        <  2:259314/151‑1 (MQ=255)                                                                                     agaTCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATc                                                       >  2:330355/1‑150 (MQ=255)                                                                                       aTCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAAc                                                    <  2:298380/151‑1 (MQ=255)                                                                                          aGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGAc                                                 <  2:265827/151‑1 (MQ=255)                                                                                                gagaACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACTTtgtg                                           <  2:284106/151‑1 (MQ=255)                                                                                                         gCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACTTTGTGCCACCACTc                                  <  2:503406/151‑1 (MQ=255)                                                                                                            ttAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACTTTGTGCCACCACTCa                                 <  1:223534/149‑1 (MQ=255)                                                                                                             tAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACTTTGTGCCACCACTCAGTc                              >  2:444817/1‑151 (MQ=255)                                                                                                                         tCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACTTTGTGCCACCACTCAGTCTACGGAAAaag                   <  1:338478/150‑1 (MQ=255)                                                                                                                         tCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGGGTTAAGTTATCAACGACTTTGTGCCACCACTCAGTCTACGGaaaa                     >  2:175737/1‑148 (MQ=255)                                                                                                                                cGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACTTTGTGCCACCACTCAGTCTACGGAAAAAGAAGAagag           <  2:15368/151‑1 (MQ=255)                                                                                                                                 ggTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACTTTGTGCAACCACTCAGTCTACGGAAAAagaagaagatt          <  1:236772/151‑3 (MQ=255)                                                                                                                                      ccGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACTTTGTGCCACCACTCAGTCTACGGAAAAAGAAGAAGAGTGGGa      <  1:740617/150‑1 (MQ=255)                                                                                                                                         tACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACTTTGTGCCACCACTCAGTCTACGGAAAAAGAAGAAGAGTGGGATCCg  >  2:492586/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                              TGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTGGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACTTTGTGCCACCACTCAGTCTACGGAAAAAGAAGAAGAGTGGGATCCG  >  NC_000913/2209131‑2209417

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