Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,210,944	G→A	pseudogene (1845/2001 nt)	yehQ →	pseudogene

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,210,944	0	G	A	100.0%	193.3 / NA	59	pseudogene (1845/2001 nt)	yehQ	pseudogene
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/28);  total (31/28)

CCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAGTTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCC  >  NC_000913/2210816‑2211078                                                                                                                                 |                                                                                                                                       ccTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGa                                                                                                                              >  1:258688/1‑139 (MQ=255)                          cGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTa                                                                                                     <  2:977060/139‑1 (MQ=255)                                   cTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGc                                                                                            >  2:597965/1‑139 (MQ=255)                                      atcaGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGcaca                                                                                         <  1:617811/139‑1 (MQ=255)                                       tcaGGTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACAt                                                                                        <  2:258688/139‑1 (MQ=255)                                                gATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTGTTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGa                                                                               >  1:922762/1‑139 (MQ=255)                                                   aCGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTAt                                                                                                    <  2:808681/115‑1 (MQ=255)                                                   aCGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTAt                                                                                                    >  1:808681/1‑115 (MQ=255)                                                           aGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATgtgt                                                                                    <  1:217510/123‑1 (MQ=255)                                                           aGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATgtgt                                                                                    >  2:217510/1‑123 (MQ=255)                                                                 tACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAAt                                                                                                                            <  1:1106671/77‑1 (MQ=255)                                                                 tACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAAt                                                                                                                            <  2:903096/77‑1 (MQ=255)                                                                 tACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAAt                                                                                                                            >  2:1106671/1‑77 (MQ=255)                                                                 tACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTAAGCCATAATTGCTGGCGAAt                                                                                                                            >  1:903096/1‑77 (MQ=255)                                                                         cAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAgg                                                      <  1:597965/139‑1 (MQ=255)                                                                           aCTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCa                                                    <  2:922762/139‑1 (MQ=255)                                                                            cTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGaa                                                                                                                             >  1:1078552/1‑65 (MQ=255)                                                                            cTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGaa                                                                                                                             <  2:1078552/65‑1 (MQ=255)                                                                              gCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGAc                                                                                                                         <  1:654279/67‑1 (MQ=255)                                                                              gCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGAc                                                                                                                         >  2:654279/1‑67 (MQ=255)                                                                              gCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGt                                                                                                              <  2:1221944/78‑1 (MQ=255)                                                                              gCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGt                                                                                                              >  1:1221944/1‑78 (MQ=255)                                                                              gCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACaa                                                                            <  1:945815/112‑1 (MQ=255)                                                                              gCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACaa                                                                            >  2:945815/1‑112 (MQ=255)                                                                              gCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACgg                                                 >  2:631325/1‑139 (MQ=255)                                                                               cTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGc                                                                                                                                    <  2:1061380/55‑1 (MQ=255)                                                                               cTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGc                                                                                                                                    >  1:1061380/1‑55 (MQ=255)                                                                               cTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACggg                                                >  2:620990/1‑139 (MQ=255)                                                                               cTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACggg                                                <  2:1065650/139‑1 (MQ=255)                                                                                 ccATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACaa                                                                            <  1:1275166/109‑1 (MQ=255)                                                                                 ccATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACaa                                                                            >  2:1275166/1‑109 (MQ=255)                                                                                    tAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGaa                                           >  2:440301/1‑139 (MQ=255)                                                                                            ggCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGaa                                   >  1:703308/1‑139 (MQ=255)                                                                                             gcgcTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGc                                                     >  1:618450/1‑120 (MQ=255)                                                                                             gcgcTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGc                                                     <  2:618450/120‑1 (MQ=255)                                                                                             gcgcTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAAt                                  >  2:416749/1‑139 (MQ=255)                                                                                              cgcTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCg                                                                                                                               <  1:480000/45‑1 (MQ=255)                                                                                              cgcTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCg                                                                                                                               >  2:480000/1‑45 (MQ=255)                                                                                               gcTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATaa                                <  2:1020093/139‑1 (MQ=255)                                                                                                 tATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACg                              <  1:416749/139‑1 (MQ=255)                                                                                                   tCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTAt                                                                                                    <  2:1249199/67‑1 (MQ=255)                                                                                                   tCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTAt                                                                                                    >  1:1249199/1‑67 (MQ=255)                                                                                                             agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCgg                                                                                                               >  1:22674/1‑46 (MQ=255)                                                                                                             agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCgg                                                                                                               <  2:22674/46‑1 (MQ=255)                                                                                                             agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGcaca                                                                                         >  1:799086/1‑68 (MQ=255)                                                                                                             agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGcaca                                                                                         <  2:799086/68‑1 (MQ=255)                                                                                                             agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACaa                                                                            <  2:976136/81‑1 (MQ=255)                                                                                                             agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACaa                                                                            >  1:976136/1‑81 (MQ=255)                                                                                                             agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTtttgcttt                  >  1:346469/1‑139 (MQ=255)                                                                                                               tGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTtttgctttgc                >  1:304097/1‑139 (MQ=255)                                                                                                                 cTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGa              >  1:742970/1‑139 (MQ=255)                                                                                                                  tATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGaa             <  1:440301/139‑1 (MQ=255)                                                                                                                       gTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAAt                                      >  1:272082/1‑109 (MQ=255)                                                                                                                       gTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAAt                                      <  2:272082/109‑1 (MQ=255)                                                                                                                         cAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGTACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAattat      >  1:388835/1‑139 (MQ=255)                                                                                                                           ggCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTc    >  1:216378/1‑139 (MQ=255)                                                                                                                            gCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTcc   >  2:738977/1‑139 (MQ=255)                                                                                                                             cATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTccc  <  2:703308/139‑1 (MQ=255)                                                                                                                              aTAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCgg                                                                                                 <  1:988300/43‑1 (MQ=255)                                                                                                                                 |                                                                                                                                       CCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAGTTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCC  >  NC_000913/2210816‑2211078

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