Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,786,729	G→A	pseudogene (333/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,786,729	0	G	A	100.0%	206.3 / NA	65	pseudogene (333/2253 nt)	ygaQ	uncharacterized protein, putative enzyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (38/27);  total (38/27)

TAATGATCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCT  >  NC_000913/2786606‑2786858                                                                                                                            |                                                                                                                                  taATGATCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGAt                                                                                                                    >  1:156182/1‑139 (MQ=255)   aTGATCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGTTACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATaa                                                                                                                  <  1:1301635/139‑2 (MQ=255)       tCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCaa                                                                                                                          >  2:749361/1‑127 (MQ=255)       tCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCaa                                                                                                                          <  1:749361/127‑1 (MQ=255)       tCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGa                                                                                                             >  2:66480/1‑139 (MQ=255)                     aTTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGa                                                                                                         <  1:517208/129‑1 (MQ=255)                     aTTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGa                                                                                                         >  2:517208/1‑129 (MQ=255)                              tttCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGgt                                                                                      <  2:156182/139‑1 (MQ=255)                                tCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGgtgt                                                                                    >  1:1555914/1‑139 (MQ=255)                                 cccTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGtt                                                                                   >  1:112461/1‑139 (MQ=255)                                  ccTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGttt                                                                                  <  1:66480/139‑1 (MQ=255)                                            gATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATgg                                                                        <  2:1555914/139‑1 (MQ=255)                                                 aaCTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGaaa                                                                                                       <  2:476397/103‑1 (MQ=255)                                                 aaCTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGaaa                                                                                                       >  1:476397/1‑103 (MQ=255)                                                  aCTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGc                                                                  >  1:755593/1‑139 (MQ=255)                                                      aCCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGttt                                                                                  >  2:100048/1‑119 (MQ=255)                                                      aCCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGttt                                                                                  <  1:100048/119‑1 (MQ=255)                                                      aCCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGaa                                                              >  1:897350/1‑139 (MQ=255)                                                      aCCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGaa                                                              >  1:517471/1‑139 (MQ=255)                                                      aCCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGaa                                                              >  2:105338/1‑139 (MQ=255)                                                      aCCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGaa                                                              >  2:127938/1‑139 (MQ=255)                                                             gTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCa                                                       >  1:322898/1‑139 (MQ=255)                                                               tGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                     <  2:517471/139‑1 (MQ=255)                                                                 aaTAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAg                                                                                          >  1:559181/1‑100 (MQ=255)                                                                 aaTAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAg                                                                                          <  2:559181/100‑1 (MQ=255)                                                                      agcagcACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATta                                              >  2:386932/1‑139 (MQ=255)                                                                           cacaGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCg                                                                                        >  2:1066572/1‑92 (MQ=255)                                                                           cacaGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCg                                                                                        <  1:1066572/92‑1 (MQ=255)                                                                           cacaGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCt                                         >  2:656068/1‑139 (MQ=255)                                                                              aGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGAc                                      <  1:656068/139‑1 (MQ=255)                                                                                cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACGTAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTc                                                        <  1:1022640/119‑1 (MQ=255)                                                                                cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTc                                                        >  2:1022640/1‑119 (MQ=255)                                                                                cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                     >  1:1546055/1‑122 (MQ=255)                                                                                cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                     <  2:1546055/122‑1 (MQ=255)                                                                                cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACaa                                    >  2:1381010/1‑139 (MQ=255)                                                                                cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGCTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACaa                                    >  1:111498/1‑139 (MQ=255)                                                                                   aGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAAcc                                 <  2:112461/139‑1 (MQ=255)                                                                                          cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATaa                                                                                                                  >  2:223234/1‑51 (MQ=255)                                                                                          cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATaa                                                                                                                  <  1:223234/52‑2 (MQ=255)                                                                                             cgcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCcaca                                                 >  1:791092/1‑113 (MQ=255)                                                                                             cgcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCcaca                                                 <  2:791092/113‑1 (MQ=255)                                                                                             cgcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCaaaa                       >  2:388709/1‑139 (MQ=255)                                                                                             cgcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCaa                         <  2:989376/137‑1 (MQ=255)                                                                                             cgcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCaa                         >  1:989376/1‑137 (MQ=255)                                                                                               ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGaaa                                                                                                       <  2:927041/57‑1 (MQ=255)                                                                                               ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGaaa                                                                                                       >  1:927041/1‑57 (MQ=255)                                                                                               ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAg                                                                                          >  2:833897/1‑70 (MQ=255)                                                                                               ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAg                                                                                          <  1:833897/70‑1 (MQ=255)                                                                                               ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGa                     >  2:695494/1‑139 (MQ=255)                                                                                                       gATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                     >  2:406530/1‑99 (MQ=255)                                                                                                       gATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                     <  1:406530/99‑1 (MQ=255)                                                                                                       gATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATg             >  2:1258589/1‑139 (MQ=255)                                                                                                            gagaAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCa        >  1:660447/1‑139 (MQ=255)                                                                                                            gagaAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCa        >  1:124052/1‑139 (MQ=255)                                                                                                             agaaagaaAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCat       <  1:388709/139‑1 (MQ=255)                                                                                                              gaaagaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCata      <  1:695494/139‑1 (MQ=255)                                                                                                               aaagaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCt                                                                 <  2:1000094/79‑1 (MQ=255)                                                                                                               aaagaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCt                                                                 >  1:1000094/1‑79 (MQ=255)                                                                                                                 agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGaaa                                                                                                       >  1:354604/1‑39 (MQ=37)                                                                                                                 agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGaaa                                                                                                       <  2:354604/39‑1 (MQ=25)                                                                                                                 agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCa                          <  2:588254/116‑1 (MQ=255)                                                                                                                 agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCa                          >  1:588254/1‑116 (MQ=255)                                                                                                                 agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGc   >  1:908265/1‑139 (MQ=255)                                                                                                                 agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGc   >  2:172507/1‑139 (MQ=255)                                                                                                                  gaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCt  <  2:755593/139‑1 (MQ=255)                                                                                                                            |                                                                                                                                  TAATGATCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCT  >  NC_000913/2786606‑2786858

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: ‑