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

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

TTATGCCGTTGAACGCATGGCGGCGCTTTTGCAACAAACAGACGACCCCGTGTGTCTGGTTTCTGGCTTTGTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCT  >  NC_000913/2210644‑2210912                                                                                                                                        |                                                                                                                                      ttATGCCGTTGAACGCATGGCGGCGCTTTTGCAACAAACAGACGACCCCGTGTGTCTGGTTTCTGGCTTTGTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAgcg                                                                                                                                    >  2:1000391/1‑139 (MQ=255)         ttGAACGCATGGCGGCGCTTTTGCAACAAACAGACGACCCCGTGTGTCTGGTTTCTGGCTTTGTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGc                                                                                                                            >  1:276944/1‑139 (MQ=255)          tGAACGCATGGCGGCGCTTTTGCAACAAACAGACGACCCCGTGTGTCTGGTTTCTGGCTTTGTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCa                                                                                                                           <  2:1028874/139‑1 (MQ=255)          tGAACGCATGGCGGCGCTTTTGCAACAAACAGACGACCCCGTGTGTCTGGTTTCTGGCTTTGTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCa                                                                                                                           >  2:679276/1‑139 (MQ=255)          tGAACGCATGGCGGCGCTTTTGCAACAAACAGACGACCCCGTGTGTCTGGTTTCTGGCTTTGTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCa                                                                                                                           <  1:189068/139‑1 (MQ=255)          tGAACGCATGGCGGCGCTTTTGCAACAAACAGACGACCCCGTGTGTCTGGTTTCTGGCTTTGTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCa                                                                                                                           >  2:189068/1‑139 (MQ=255)              cGCATGGCGGCGCTTTTGCAACAAACAGACGACCCCGTGTGTCTGGTTTCTGGCTTTGTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGaaa                                                                                                                       >  1:342509/1‑139 (MQ=255)                  tGGCGGCGCTTTTGCAACAAACAGACGACCCCGTGTGTCTGGTTTCTGGCTTTGTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGTGCTGGACGCAGAAACTAc                                                                                                                   <  2:1030164/139‑1 (MQ=255)                        cgcTTTTGCAACAAACAGACGACCCCGTGTGTCTGGTTTCTGGCTTTGTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCtgt                                                                                                             >  2:600340/1‑139 (MQ=255)                          cTTTTGCAACAAACAGACGACCCCGTGTGTCTGGTTTCTGGCTTTGTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTgg                                                                                                           <  2:342509/139‑1 (MQ=255)                            tttGCAACAAACAGACGACCCCGTGTGTCTGGTTTCTGGCTTTGTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCa                                                                                                         <  1:290844/139‑1 (MQ=255)                                          acgacCCCGTGTGTCTGGTTTCTGGCTTTGTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGc                                                                                           >  2:263516/1‑139 (MQ=255)                                           cgacCCCGTGTGTCTGGTTTCTGGCTTTGTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGcc                                                                                          <  1:645133/139‑1 (MQ=255)                                               cccGTGTGTCTGGTTTCTGGCTTTGTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCg                                                                                      >  1:427914/1‑139 (MQ=255)                                                       tCTGGTTTCTGGCTTTGTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCt                                                                              >  2:293077/1‑139 (MQ=255)                                                          ggTTTCTGGCTTTGTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCtgtg                                                                           <  2:276944/139‑1 (MQ=255)                                                          ggTTTCTGGCTTTGTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCtgtg                                                                           <  1:263516/139‑1 (MQ=255)                                                             ttCTGGCTTTGTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGc                                                                                <  2:865346/131‑1 (MQ=255)                                                             ttCTGGCTTTGTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGc                                                                                >  1:865346/1‑131 (MQ=255)                                                               cTGGCTTTGTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGa                                                                                                                         <  1:900621/88‑1 (MQ=255)                                                               cTGGCTTTGTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGa                                                                                                                         >  2:900621/1‑88 (MQ=255)                                                               cTGGCTTTGTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTc                                                                      >  1:1035960/1‑139 (MQ=255)                                                                 ggCTTTGTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAgcg                                                                                                                                    >  1:193777/1‑75 (MQ=255)                                                                 ggCTTTGTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAgcg                                                                                                                                    <  2:193777/75‑1 (MQ=255)                                                                 ggCTTTGTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGt                                                                                     <  2:609116/122‑1 (MQ=255)                                                                 ggCTTTGTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGt                                                                                     >  1:609116/1‑122 (MQ=255)                                                                 ggCTTTGTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTa                                                                    >  2:33737/1‑139 (MQ=255)                                                                        tCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCtca                                                             <  1:1000391/139‑1 (MQ=255)                                                                                gttgAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGt                                                        >  1:824123/1‑136 (MQ=255)                                                                                gttgAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGt                                                        <  2:824123/136‑1 (MQ=255)                                                                                 ttgAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCt                                                    >  2:920017/1‑139 (MQ=255)                                                                                         cAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCtcat                                                            <  2:968013/123‑1 (MQ=255)                                                                                         cAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCtcat                                                            >  1:968013/1‑123 (MQ=255)                                                                                         cAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCt                                            >  2:898386/1‑139 (MQ=255)                                                                                          aaTTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTg                                           <  1:77622/139‑1 (MQ=255)                                                                                           aTTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGc                                          <  2:1163570/139‑1 (MQ=255)                                                                                                cacTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGc                                     <  2:487051/139‑1 (MQ=255)                                                                                                   tGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGtt                                  <  2:625734/139‑1 (MQ=255)                                                                                                   tGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGtt                                  <  1:920017/139‑1 (MQ=255)                                                                                                     gAACCACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTAc                                <  2:331019/139‑1 (MQ=255)                                                                                                        ccACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTg                                           >  2:1048785/1‑125 (MQ=255)                                                                                                        ccACGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTg                                           <  1:1048785/125‑1 (MQ=255)                                                                                                           cGGGTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTc                          >  2:1104369/1‑139 (MQ=255)                                                                                                             ggTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTa                             >  2:495339/1‑134 (MQ=255)                                                                                                             ggTGATGATGACAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTa                             <  1:495339/134‑1 (MQ=255)                                                                                                               tgatgatgaCAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACt                      >  2:29527/1‑139 (MQ=255)                                                                                                                gatgatgaCAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAGCTg                     >  1:798747/1‑139 (MQ=255)                                                                                                                   gatgaCAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGtt                                  >  1:763989/1‑123 (MQ=255)                                                                                                                   gatgaCAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGtt                                  <  2:763989/123‑1 (MQ=255)                                                                                                                      gaCAAAAACCCGTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCAt               <  2:689062/139‑1 (MQ=255)                                                                                                                                 gTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGgcg    <  1:293077/139‑1 (MQ=255)                                                                                                                                 gTGCCTGAGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGgcg    <  1:1104369/139‑1 (MQ=255)                                                                                                                                  tGCCTGAGCGCTGGACGCAGAAGCTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATC‑GTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCt  <  2:1005562/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                      TTATGCCGTTGAACGCATGGCGGCGCTTTTGCAACAAACAGACGACCCCGTGTGTCTGGTTTCTGGCTTTGTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCT  >  NC_000913/2210644‑2210912

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