Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,786,746	Δ1 bp	pseudogene (350/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,746	0	C	.	100.0%	253.4 / NA	59	pseudogene (350/2253 nt)	ygaQ	uncharacterized protein, putative enzyme
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (30/29);  total (30/29)

TTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCTATTAACATTCCTGAC  >  NC_000913/2786616‑2786873                                                                                                                                   |                                                                                                                                ttAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGa                                                                                                                        >  2:477737/1‑139 (MQ=255) ttAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGa                                                                                                                        >  2:346162/1‑139 (MQ=255)   aaaaaaCGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGaaa                                                                                                                      <  1:733553/139‑1 (MQ=255)     aaaaCGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAaga                                                                                                                    <  2:603778/139‑1 (MQ=255)           aTTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGAt                                                                                                              >  2:656531/1‑139 (MQ=255)           aTTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGAt                                                                                                              >  2:114288/1‑139 (MQ=255)                        ccTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGttt                                                                                                 <  1:346162/139‑1 (MQ=255)                              tcCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGt                                                                                           >  2:702061/1‑139 (MQ=255)                                        aCTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCg                                                                                                       >  2:479303/1‑117 (MQ=255)                                        aCTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCg                                                                                                       <  1:479303/117‑1 (MQ=255)                                        aCTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGgtg                                                                                                    <  2:695826/120‑1 (MQ=255)                                        aCTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGgtg                                                                                                    >  1:695826/1‑120 (MQ=255)                                            aCCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGaa                                                                             >  2:684166/1‑139 (MQ=255)                                             ccGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAAc                                                                            >  2:95501/1‑139 (MQ=255)                                             ccGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAAc                                                                            >  2:135964/1‑139 (MQ=255)                                                  agTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTc                                                                       >  2:448362/1‑139 (MQ=255)                                                   gTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCa                                                                      >  1:97922/1‑139 (MQ=255)                                                    ttGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAc                                                                     >  1:579893/1‑139 (MQ=255)                                                     tGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                                    <  1:565875/139‑1 (MQ=255)                                                     tGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                                    <  1:95501/139‑1 (MQ=255)                                                              cagcaCAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATtata                                                           <  2:673008/139‑1 (MQ=255)                                                                 cacaGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCt                                                        >  2:635688/1‑139 (MQ=255)                                                                      cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAc                                                                     <  2:361646/121‑1 (MQ=255)                                                                      cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAc                                                                     >  1:361646/1‑121 (MQ=255)                                                                      cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACaa                                                   >  1:786400/1‑139 (MQ=255)                                                                         aGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTAGTGCTGAACATTTCACCCACATTATATCTGACAAAcc                                                >  1:240928/1‑139 (MQ=255)                                                                                cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCca                                                                  <  2:599793/114‑1 (MQ=255)                                                                                cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCca                                                                  >  1:599793/1‑114 (MQ=255)                                                                                 aGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCaa                                        <  1:684166/139‑1 (MQ=255)                                                                                  gcgcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTatat                                                          >  1:314735/1‑120 (MQ=255)                                                                                  gcgcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTatat                                                          <  2:314735/120‑1 (MQ=255)                                                                                  gcgcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCaaa                                       >  1:330968/1‑139 (MQ=255)                                                                                   cgcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCaaaa                                      >  1:410077/1‑139 (MQ=255)                                                                                   cgcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCaaaa                                      >  2:217734/1‑139 (MQ=255)                                                                                     ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCa                                         <  1:481008/134‑1 (MQ=255)                                                                                     ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCa                                         >  2:481008/1‑134 (MQ=255)                                                                                     ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGa                                    >  1:122481/1‑139 (MQ=255)                                                                                       aGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATa                                  <  1:476548/139‑1 (MQ=255)                                                                                        gACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATaa                                 <  1:217734/139‑1 (MQ=255)                                                                                        gACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATaa                                 <  2:240928/139‑1 (MQ=255)                                                                                            ggATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAatat                             <  1:674890/139‑1 (MQ=255)                                                                                                aGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATg                            >  1:197708/1‑136 (MQ=255)                                                                                                aGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATg                            <  2:197708/136‑1 (MQ=255)                                                                                                  gagaAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCa                       >  1:539424/1‑139 (MQ=255)                                                                                                   agaaagaaAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCat                      <  1:114288/139‑1 (MQ=255)                                                                                                   agaaagaaAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCat                      <  2:75170/139‑1 (MQ=255)                                                                                                   agaaagaaAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCat                      <  2:116072/139‑1 (MQ=255)                                                                                                    gaaagaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCata                     <  1:135964/139‑1 (MQ=255)                                                                                                    gaaagaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCata                     <  1:702061/139‑1 (MQ=255)                                                                                                     aaagaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCatat                    <  1:635688/139‑1 (MQ=255)                                                                                                       agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAat                               >  2:400431/1‑126 (MQ=255)                                                                                                       agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAat                               <  1:400431/126‑1 (MQ=255)                                                                                                       agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGc                  >  2:207179/1‑139 (MQ=255)                                                                                                             aCATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCTATTaa            <  2:330968/139‑1 (MQ=255)                                                                                                               ataagTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCTATTAACa          <  2:539424/135‑1 (MQ=255)                                                                                                                 aagTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCTATTAACAtt        <  2:786400/137‑1 (MQ=255)                                                                                                                       aTCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCTATTAACATTCCTGAc  >  2:564799/1‑139 (MQ=255)                                                                                                                             gagaTA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATtat                                                            <  2:402332/75‑1 (MQ=255)                                                                                                                             gagaTA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATtat                                                            >  1:402332/1‑75 (MQ=255)                                                                                                                                   |                                                                                                                                TTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCTATTAACATTCCTGAC  >  NC_000913/2786616‑2786873

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