Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,860,416	G→A	*219* (TAG→TAA)	pphB →	serine/threonine‑specific protein phosphatase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,860,416	0	G	A	100.0%	179.1 / NA	56	*219* (TAG→TAA)	pphB	serine/threonine‑specific protein phosphatase 2
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/29);  total (27/29)

AATAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAGTCTCATGCTTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACCGTCGTCTGTAATA  >  NC_000913/2860287‑2860551                                                                                                                                  |                                                                                                                                        aaTAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGct                                                                                                                                <  1:61370/139‑1 (MQ=255) aaTAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGct                                                                                                                                <  1:574032/139‑1 (MQ=255)      aCGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑tcttct                                                                                                                          >  1:1713855/1‑139 (MQ=255)                             ggACATATGATGTTTGATAACATTCAGGCGTTAGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGgtgt                                                                                                   >  1:1224958/1‑139 (MQ=255)                             ggACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCtgtg                                                                                                                       <  2:462031/119‑1 (MQ=255)                             ggACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCtgtg                                                                                                                       >  1:462031/1‑119 (MQ=255)                             ggACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATATATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGgtgt                                                                                                   >  1:1453584/1‑139 (MQ=255)                               aCATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTa                                                                                                 <  2:1713855/139‑1 (MQ=255)                                             aTAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGc                                                                                                                   <  2:1916212/107‑1 (MQ=255)                                             aTAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGc                                                                                                                   >  1:1916212/1‑107 (MQ=255)                                                 cATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATg                                                                                                                                  <  2:382221/89‑1 (MQ=255)                                                 cATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATg                                                                                                                                  >  1:382221/1‑89 (MQ=255)                                                  aTTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAAtctc                                                                                                                                     >  2:879308/1‑85 (MQ=255)                                                  aTTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAAtctc                                                                                                                                     <  1:879308/85‑1 (MQ=255)                                                   ttCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAAtc                                                                                                                                       <  1:1992837/82‑1 (MQ=255)                                                   ttCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAAtc                                                                                                                                       >  2:1992837/1‑82 (MQ=255)                                                          gTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTa                                                                      <  2:428527/139‑1 (MQ=255)                                                             cGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCtg                                                                                                                         <  1:193856/85‑1 (MQ=255)                                                             cGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCtg                                                                                                                         >  2:193856/1‑85 (MQ=255)                                                                   ccAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAag                                                             <  2:1453584/139‑1 (MQ=255)                                                                    cAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAaga                                                            <  2:86198/139‑1 (MQ=255)                                                                      gATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGc                                                                                                                   >  2:817406/1‑82 (MQ=255)                                                                      gATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGc                                                                                                                   <  1:817406/82‑1 (MQ=255)                                                                        tttATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGATACCAGCAAAgagaga                                                        >  1:549827/1‑139 (MQ=255)                                                                             atTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAAccc                                                                                                        >  2:1748986/1‑86 (MQ=255)                                                                             atTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAAccc                                                                                                        <  1:1748986/86‑1 (MQ=255)                                                                              ttGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGc                                                  >  2:1856543/1‑139 (MQ=255)                                                                                 aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGct                                                                                                                                >  2:811795/1‑59 (MQ=255)                                                                                 aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGct                                                                                                                                <  1:811795/59‑1 (MQ=255)                                                                                 aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTc                                               >  1:1108511/1‑139 (MQ=255)                                                                                     cGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCa                                           >  1:1254245/1‑139 (MQ=255)                                                                                      ggATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTg                                                                                              <  2:254118/87‑1 (MQ=255)                                                                                      ggATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTg                                                                                              >  1:254118/1‑87 (MQ=255)                                                                                        aTCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAgcg                                        >  2:684050/1‑139 (MQ=255)                                                                                           gCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAAccc                                                                                                        >  2:1236559/1‑72 (MQ=255)                                                                                           gCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAAccc                                                                                                        <  1:1236559/72‑1 (MQ=255)                                                                                            ccGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCtgt                                                                                                                        <  1:1890170/55‑1 (MQ=255)                                                                                            ccGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCtgt                                                                                                                        >  2:1890170/1‑55 (MQ=255)                                                                                                   gcgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCa                             <  2:1108511/139‑1 (MQ=255)                                                                                                    cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGct                                                                                                                                <  1:610112/40‑1 (MQ=255)                                                                                                    cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGct                                                                                                                                >  2:610112/1‑40 (MQ=255)                                                                                                    cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGt                                                                                                             >  2:75092/1‑58 (MQ=255)                                                                                                    cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGt                                                                                                             <  1:75092/58‑1 (MQ=255)                                                                                                    cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGAt                                                       >  2:700236/1‑112 (MQ=255)                                                                                                    cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGAt                                                       <  1:700236/112‑1 (MQ=255)                                                                                                    cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAgcgcgc                                     >  1:1717780/1‑130 (MQ=255)                                                                                                    cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAgcgcgc                                     <  2:1717780/130‑1 (MQ=255)                                                                                                      ggcggcTGTCATTTTATAAAATAAAGTAATCTCATGctctt                                                                                                                             >  1:337613/1‑38 (MQ=255)                                                                                                      ggcggcTGTCATTTTATAAAATAAAGTAATCTCATGctctt                                                                                                                             <  2:337613/41‑4 (MQ=255)                                                                                                      ggcggcTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGAt                                                       <  1:1487647/110‑1 (MQ=255)                                                                                                      ggcggcTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGAt                                                       >  2:1487647/1‑110 (MQ=255)                                                                                                               cATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCaaa                 <  1:684050/139‑1 (MQ=255)                                                                                                               cATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGACATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCaaa                 <  2:1254245/139‑1 (MQ=255)                                                                                                                  tttATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACcg              <  2:1731441/139‑1 (MQ=255)                                                                                                                           taaaGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACcgtc            <  2:668641/132‑1 (MQ=255)                                                                                                                           taaaGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACcgtc            >  1:668641/1‑132 (MQ=255)                                                                                                                              agtaatCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACcgtcgt          >  1:670145/1‑131 (MQ=255)                                                                                                                              agtaatCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACcgtcgt          <  2:670145/131‑1 (MQ=255)                                                                                                                              agtaatCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACCGTCGTCTGtaata  >  2:1811906/1‑139 (MQ=255)                                                                                                                                  |                                                                                                                                        AATAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAGTCTCATGCTTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACCGTCGTCTGTAATA  >  NC_000913/2860287‑2860551

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