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

CTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAGTTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCC  >  NC_000913/2210817‑2211078                                                                                                                                |                                                                                                                                       cTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGaa                                                                                                                             <  2:680142/139‑1 (MQ=255)           cGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCAcc                                                                                                                   >  1:848456/1‑139 (MQ=255)           cGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCAcc                                                                                                                   <  2:62760/139‑1 (MQ=255)              tttGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACcgcg                                                                                                                >  1:1217266/1‑139 (MQ=255)               ttGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCgg                                                                                                               <  2:1144628/139‑1 (MQ=255)                 gCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTg                                                                                                             <  1:1122969/139‑1 (MQ=255)                                  cTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGc                                                                                            >  1:277557/1‑139 (MQ=255)                                  cTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGc                                                                                            >  1:1373659/1‑139 (MQ=255)                                     atcaGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGcaca                                                                                         >  1:825969/1‑139 (MQ=255)                                               gATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGa                                                                               <  2:277557/139‑1 (MQ=255)                                               gATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGa                                                                               <  2:728061/139‑1 (MQ=255)                                                  aCGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACaa                                                                            <  2:1545766/139‑1 (MQ=255)                                                              gTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACa                                                                <  2:848456/139‑1 (MQ=255)                                                              gTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACa                                                                <  2:1373659/139‑1 (MQ=255)                                                                tACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGAACAATTTCGTAATACAGa                                                              >  1:847138/1‑139 (MQ=255)                                                                           cTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCAc                                                   >  2:1557989/1‑139 (MQ=255)                                                                             gCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACgg                                                 >  2:1480987/1‑139 (MQ=255)                                                                              cTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACggg                                                >  1:571250/1‑139 (MQ=255)                                                                                ccATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTa                                              >  1:115110/1‑139 (MQ=255)                                                                                 cATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAg                                             <  2:115110/139‑1 (MQ=255)                                                                                   tAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGa                                                                               <  2:1342101/103‑1 (MQ=255)                                                                                   tAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGa                                                                               >  1:1342101/1‑103 (MQ=255)                                                                                   tAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGaa                                           >  1:711035/1‑139 (MQ=255)                                                                                   tAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGaa                                           >  2:1177813/1‑139 (MQ=255)                                                                                        gCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAAt                                                                                                                            <  1:771145/53‑1 (MQ=255)                                                                                        gCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAAt                                                                                                                            >  2:771145/1‑53 (MQ=255)                                                                                         cTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGaa                                                                                                                             <  2:171987/51‑1 (MQ=255)                                                                                         cTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGaa                                                                                                                             >  1:171987/1‑51 (MQ=255)                                                                                         cTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGacac                                                                                          >  2:169993/1‑86 (MQ=255)                                                                                         cTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGacac                                                                                          <  1:169993/86‑1 (MQ=255)                                                                                           ggCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTc                                                                                                                      <  1:1227101/56‑1 (MQ=255)                                                                                           ggCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTc                                                                                                                      >  2:1227101/1‑56 (MQ=255)                                                                                              gcTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGTGTTTTTATCGGCTGGCACATGTGTTGTGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATaa                                <  1:519057/139‑1 (MQ=255)                                                                                                tATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTGGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACg                              <  2:1053088/139‑1 (MQ=255)                                                                                                       aaCAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGtttt                       >  1:1443467/1‑139 (MQ=255)                                                                                                            agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCAc                                                                                                                    >  2:1129907/1‑41 (MQ=255)                                                                                                            agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCAc                                                                                                                    <  1:1129907/41‑1 (MQ=255)                                                                                                            agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTgg                                                                                             <  2:1426096/64‑1 (MQ=255)                                                                                                            agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTgg                                                                                             >  1:1426096/1‑64 (MQ=255)                                                                                                            agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAAtt                                                                          <  2:1192712/83‑1 (MQ=255)                                                                                                            agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAAtt                                                                          >  1:1192712/1‑83 (MQ=255)                                                                                                            agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTc                                                                        >  2:490274/1‑85 (MQ=255)                                                                                                            agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTc                                                                        <  1:490274/85‑1 (MQ=255)                                                                                                            agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTtttgcttt                  >  2:1303082/1‑139 (MQ=255)                                                                                                            agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTtttgcttt                  >  1:864654/1‑139 (MQ=255)                                                                                                               gCTATTGGTCAAGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCg               <  1:1480987/139‑1 (MQ=255)                                                                                                                 tATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCgg                                                                                                 <  1:564001/55‑1 (MQ=255)                                                                                                                 tATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCgg                                                                                                 >  2:564001/1‑55 (MQ=255)                                                                                                                 tATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAg                                             <  1:466696/107‑1 (MQ=255)                                                                                                                 tATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAg                                             >  2:466696/1‑107 (MQ=255)                                                                                                                 tATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGaa             <  2:409784/139‑1 (MQ=255)                                                                                                                 tATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGaa             <  2:571250/139‑1 (MQ=255)                                                                                                                 tATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGaa             <  2:1443467/139‑1 (MQ=255)                                                                                                                      gTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAatt        <  2:825969/139‑1 (MQ=255)                                                                                                                        cAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAattat      >  1:708020/1‑139 (MQ=255)                                                                                                                            cATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTccc  <  1:1303082/139‑1 (MQ=255)                                                                                                                                |                                                                                                                                       CTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAGTTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCC  >  NC_000913/2210817‑2211078

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