Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,210,811	C→T	pseudogene (1712/2001 nt)	yehQ →	pseudogene

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,210,811	0	C	T	100.0%	279.6 / NA	77	pseudogene (1712/2001 nt)	yehQ	pseudogene
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (36/41);  total (36/41)

TCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTG  >  NC_000913/2210715‑2210907                                                                                                 |                                                                                                 tCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCt                                                                                               <  2:966490/100‑1 (MQ=255)    gTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTAcc                                                                                            >  1:331030/1‑100 (MQ=255)        tgttgAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCt                                                                                        <  1:2209546/100‑1 (MQ=255)         gttgAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTg                                                                                       >  1:551488/1‑100 (MQ=255)         gttgAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTg                                                                                       <  1:396450/100‑1 (MQ=255)          ttgAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGc                                                                                      <  2:336953/100‑1 (MQ=255)           tgAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGcc                                                                                     <  2:1107925/100‑1 (MQ=255)           tgAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGcc                                                                                     >  1:2319793/1‑100 (MQ=255)            gAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCa                                                                                    >  1:1792193/1‑100 (MQ=255)            gAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCa                                                                                    >  2:562155/1‑100 (MQ=255)               gggCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCg                                                                                 >  1:406344/1‑100 (MQ=255)                  cAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGttt                                                                              <  1:1833345/100‑1 (MQ=255)                    aTTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTg                                                                            <  1:2132302/100‑1 (MQ=255)                     ttGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGc                                                                           <  1:2086156/100‑1 (MQ=255)                       gACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCt                                                                         <  1:859553/100‑1 (MQ=255)                       gACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCt                                                                         >  2:2233089/1‑100 (MQ=255)                        acacTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCtg                                                                        <  2:2066772/100‑1 (MQ=255)                           cTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGc                                                                     <  1:1150235/100‑1 (MQ=255)                             ggAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCg                                                                   >  2:2031420/1‑100 (MQ=255)                             ggAACCACGGGTGATGATGACAAAAACCCGTGCCTGCGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCg                                                                   >  1:1651559/1‑100 (MQ=255)                              gAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCg                                                                   >  1:1154736/1‑99 (MQ=255)                               aaCCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTc                                                                 <  2:731853/100‑1 (MQ=255)                                aCCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTc                                                                 >  2:1118598/1‑99 (MQ=255)                                  cACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTAc                                                              <  1:562156/100‑1 (MQ=255)                                  cACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTAc                                                              <  2:1069951/100‑1 (MQ=255)                                  cACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTAc                                                              >  2:25946/1‑100 (MQ=255)                                     tggTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGc                                                           <  1:2199797/99‑1 (MQ=255)                                     gggTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGc                                                           <  1:2233097/100‑1 (MQ=255)                                     gggTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGCGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGc                                                           <  2:1714768/100‑1 (MQ=255)                                      ggTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGca                                                          >  2:2199789/1‑99 (MQ=255)                                       gTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCtc                                                         <  2:535358/100‑1 (MQ=255)                                         gatgatgaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCtcat                                                       >  1:1061825/1‑100 (MQ=255)                                          atgatgaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCtcatc                                                      >  1:1033132/1‑100 (MQ=255)                                           tgatgaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCtcatca                                                     >  2:1680859/1‑100 (MQ=255)                                            gatgaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAg                                                    <  1:1696857/100‑1 (MQ=255)                                             atgaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGt                                                   >  2:244759/1‑100 (MQ=255)                                              tgaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGtt                                                  >  2:571188/1‑100 (MQ=255)                                              tgaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGtt                                                  >  2:256126/1‑100 (MQ=255)                                               gaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTg                                                 >  1:940295/1‑100 (MQ=255)                                                aCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGc                                                <  1:1740508/100‑1 (MQ=255)                                                aCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGc                                                >  2:889559/1‑100 (MQ=255)                                                aCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGc                                                >  2:1554270/1‑100 (MQ=255)                                                aCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGc                                                >  2:1740503/1‑100 (MQ=255)                                                 cAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCt                                               <  1:889561/100‑1 (MQ=255)                                                 cAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCt                                               >  1:1814703/1‑100 (MQ=255)                                                 cAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCt                                               <  2:1651554/100‑1 (MQ=255)                                                  aaaaaCCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTg                                              <  2:331028/100‑1 (MQ=255)                                                  aaaaaCCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTg                                              <  1:1492877/100‑1 (MQ=255)                                                    aaaCCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGAt                                            <  1:2080676/100‑1 (MQ=255)                                                    aaaCCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGAt                                            <  2:9409/100‑1 (MQ=255)                                                      aCCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATAc                                          <  1:13008/100‑1 (MQ=255)                                                       cccGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACg                                         <  1:571189/100‑1 (MQ=255)                                                         cGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCt                                       <  2:435184/100‑1 (MQ=255)                                                           tGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGc                                     >  1:1900485/1‑100 (MQ=255)                                                             ccTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCa                                   <  2:2319784/100‑1 (MQ=255)                                                               tGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAgg                                 >  2:1144312/1‑100 (MQ=255)                                                                gggCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGc                                >  2:1967368/1‑100 (MQ=255)                                                                      tGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACt                          >  2:2283498/1‑100 (MQ=255)                                                                       ggACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACtt                         <  1:244761/100‑1 (MQ=255)                                                                        gACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTa                        <  1:570292/100‑1 (MQ=255)                                                                        gACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTa                        <  2:1073247/100‑1 (MQ=255)                                                                        gACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTG‑CAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTAt                       >  2:943925/1‑100 (MQ=255)                                                                         aCGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTAt                       <  1:1554274/100‑1 (MQ=255)                                                                          cGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTAtt                      >  2:1673262/1‑100 (MQ=255)                                                                            cAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCa                    >  1:195064/1‑100 (MQ=255)                                                                            cAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCa                    >  2:2333567/1‑100 (MQ=255)                                                                            cAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCa                    >  1:332388/1‑100 (MQ=255)                                                                                 aCTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGc               <  1:758650/100‑1 (MQ=255)                                                                                  cTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCt              <  2:50579/100‑1 (MQ=255)                                                                                  cTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCt              <  2:50632/100‑1 (MQ=255)                                                                                    aCGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTcc            >  2:18901/1‑100 (MQ=255)                                                                                      gCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCAt          >  2:1666804/1‑100 (MQ=255)                                                                                         tgtgGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCCACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAAc       <  1:1118600/100‑1 (MQ=255)                                                                                          gtgGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACg      <  2:1347998/100‑1 (MQ=255)                                                                                           tgGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACgg     <  2:1825998/100‑1 (MQ=255)                                                                                           tgGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACgg     <  2:1814698/100‑1 (MQ=255)                                                                                              cACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTg  >  2:1282266/1‑100 (MQ=255)                                                                                                 |                                                                                                 TCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTG  >  NC_000913/2210715‑2210907

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