Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,788,238	G→A	pseudogene (1842/2253 nt)	ygaQ →	uncharacterized protein, putative enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,788,238	0	G	A	100.0%	247.8 / NA	75	pseudogene (1842/2253 nt)	ygaQ	uncharacterized protein, putative enzyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (36/39);  total (36/39)

AACAATTAGCAATGGCATATATAGCCGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAGCTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAGCGCATTTGAGGGTGTATCTGCAACTAAAAACGTTCTGGCCTT  >  NC_000913/2788105‑2788360                                                                                                                                      |                                                                                                                           aaCAATTAGCAATGGCATATATAGCCGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTtaa                                                                                                                       >  1:438861/1‑139 (MQ=255)     aTTAGCAATGGCATATATAGCCGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATaaa                                                                                                                   >  2:48528/1‑139 (MQ=255)        aGCAATGGCATATATAGCCGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGt                                                                                                                <  1:1009044/139‑1 (MQ=255)        aGCAATGGCATATATAGCCGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGt                                                                                                                <  2:762291/139‑1 (MQ=255)          cAATGGCATATATAGCCGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTg                                                                                                              <  2:438861/139‑1 (MQ=255)                cATATATAGCCGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTAc                                                                                                        >  2:442765/1‑139 (MQ=255)                        gCCGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATcagc                                                                                                <  1:48528/139‑1 (MQ=255)                         ccGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATaaa                                                                                                                   >  2:1658290/1‑119 (MQ=255)                         ccGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATaaa                                                                                                                   <  1:1658290/119‑1 (MQ=255)                         ccGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATcagca                                                                                               >  2:926686/1‑139 (MQ=255)                            ctctGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGctc                                                                                                           >  2:2117/1‑124 (MQ=255)                            ctctGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGctc                                                                                                           <  1:2117/124‑1 (MQ=255)                            ctctGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATcag                                                                                                 >  1:361333/1‑134 (MQ=255)                            ctctGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATcag                                                                                                 <  2:361333/134‑1 (MQ=255)                              ctGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGctct                                                                                          >  2:879767/1‑139 (MQ=255)                              ctGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGct                                                                                            >  2:1311223/1‑137 (MQ=255)                              ctGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGct                                                                                            <  1:1311223/137‑1 (MQ=255)                                       agcagGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAAt                                                                                 >  2:341021/1‑139 (MQ=255)                                            gCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCAc                                                                            <  1:879767/139‑1 (MQ=255)                                              cTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGc                                                                          <  2:351663/139‑1 (MQ=255)                                              cTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGc                                                                          <  2:837490/139‑1 (MQ=255)                                               tttAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCa                                                                         <  1:442765/139‑1 (MQ=255)                                               tttAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCa                                                                         <  2:396615/139‑1 (MQ=255)                                               tttAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCa                                                                         <  1:341021/139‑1 (MQ=255)                                                       aCTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGACGTAACTTAATAAAGGTTGCTCTACCTATcag                                                                                                 <  1:1223992/107‑1 (MQ=255)                                                       aCTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATGATGATGGGTCCACGCGTGTCGGCGGTGAGACGTAACTTAATAAAGGTTGCTCTACCTATcag                                                                                                 >  2:1223992/1‑107 (MQ=255)                                                          gggATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGa                                                              <  2:1089346/139‑1 (MQ=255)                                                           ggATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGAt                                                             >  1:1565872/1‑139 (MQ=255)                                                             aTGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATaaa                                                                                                                   >  1:1579604/1‑83 (MQ=255)                                                             aTGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATaaa                                                                                                                   <  2:1579604/83‑1 (MQ=255)                                                             aTGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATc                                                                                                   <  2:237683/99‑1 (MQ=255)                                                             aTGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATc                                                                                                   >  1:237683/1‑99 (MQ=255)                                                               ggCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCAGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAg                                                         <  1:880376/139‑1 (MQ=255)                                                                gCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGc                                                        >  1:735153/1‑139 (MQ=255)                                                                 ccAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCa                                                                         >  2:1346703/1‑121 (MQ=255)                                                                 ccAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCa                                                                         <  1:1346703/121‑1 (MQ=255)                                                                           aaaTCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGt                                                                                                                >  1:92635/1‑72 (MQ=255)                                                                           aaaTCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGt                                                                                                                <  2:92635/72‑1 (MQ=255)                                                                               cAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAgcgc                                         <  2:1565872/139‑1 (MQ=255)                                                                                aaCGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAGCGCa                                        >  2:806576/1‑139 (MQ=255)                                                                                  cGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCa                                                                             <  1:1262208/100‑1 (MQ=255)                                                                                  cGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCa                                                                             >  2:1262208/1‑100 (MQ=255)                                                                                  cGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCa                                                                         >  1:710885/1‑104 (MQ=255)                                                                                  cGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCa                                                                         <  2:710885/104‑1 (MQ=255)                                                                                        aaaTTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGctct                                                                                          <  2:1197435/81‑1 (MQ=255)                                                                                        aaaTTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGctct                                                                                          >  1:1197435/1‑81 (MQ=255)                                                                                                    atgatgGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAAtt                                                                                <  2:289680/79‑1 (MQ=255)                                                                                                    atgatgGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAAtt                                                                                >  1:289680/1‑79 (MQ=255)                                                                                                    atgatgGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAGCGCAttt                                     <  1:1629553/122‑1 (MQ=255)                                                                                                    atgatgGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAGCGCAttt                                     >  2:1629553/1‑122 (MQ=255)                                                                                                            tCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCa                                                                         <  2:691578/78‑1 (MQ=255)                                                                                                            tCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCa                                                                         >  1:691578/1‑78 (MQ=255)                                                                                                               aCGCGTGTCGGCGGTGAGGGGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAGCGCATTTGAGGGTGTATCTGCAACTAAAAACGTTc         >  1:1155154/1‑139 (MQ=255)                                                                                                                 gcgTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAGCGCATTTGAGGGTGTATCTGCAACTAAAAACGTTCTg       <  2:1155154/139‑1 (MQ=255)                                                                                                                    tgtCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTcc                                                                              <  2:627009/65‑1 (MQ=255)                                                                                                                    tgtCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTcc                                                                              >  1:627009/1‑65 (MQ=255)                                                                                                                    tgtCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCAt                                                                        >  1:726997/1‑71 (MQ=255)                                                                                                                    tgtCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCAt                                                                        <  2:726997/71‑1 (MQ=255)                                                                                                                      tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTa                                                                                                         <  1:1115406/36‑1 (MQ=255)                                                                                                                      tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATcagcag                                                                                              <  2:1447102/47‑1 (MQ=255)                                                                                                                      tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATcagcag                                                                                              >  1:1447102/1‑47 (MQ=255)                                                                                                                      tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATcag                                                                                                 <  2:709963/44‑1 (MQ=255)                                                                                                                      tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATcag                                                                                                 >  1:709963/1‑44 (MQ=255)                                                                                                                      tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGct                                                                                            <  2:662328/49‑1 (MQ=255)                                                                                                                      tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGct                                                                                            >  1:662328/1‑49 (MQ=255)                                                                                                                      tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCAc                                                                            >  2:255820/1‑65 (MQ=255)                                                                                                                      tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCAc                                                                            <  1:255820/65‑1 (MQ=255)                                                                                                                      tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCa                                                                         >  2:635876/1‑68 (MQ=255)                                                                                                                      tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCa                                                                         <  1:635876/68‑1 (MQ=255)                                                                                                                      tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATCGa                                                              >  2:1241224/1‑79 (MQ=255)                                                                                                                      tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGa                                                              <  1:1241224/79‑1 (MQ=255)                                                                                                                      tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAg                                                         <  2:1256537/84‑1 (MQ=255)                                                                                                                      tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAg                                                         >  1:1256537/1‑84 (MQ=255)                                                                                                                      tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAGCGCATTTGAGGGTGTATCTGCAACTAAAAACGTTCTGGCCtt  >  1:1440971/1‑139 (MQ=255)                                                                                                                      tCGACGGTGAGGCGTAACTTAATAAAGGTTGCTCTa                                                                                                         >  2:1115406/1‑36 (MQ=255)                                                                                                                                      |                                                                                                                           AACAATTAGCAATGGCATATATAGCCGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAGCTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAGCGCATTTGAGGGTGTATCTGCAACTAAAAACGTTCTGGCCTT  >  NC_000913/2788105‑2788360

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