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%	153.3 / NA	48	*219* (TAG→TAA)	pphB	serine/threonine‑specific protein phosphatase 2
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/23);  total (25/23)

AATAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAGTCTCATGCTTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACCGTCGTCTGTAATAA  >  NC_000913/2860287‑2860552                                                                                                                                  |                                                                                                                                         aaTAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGct                                                                                                                                 >  2:412964/1‑139 (MQ=255)       cGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCtg                                                                                                                          <  1:134469/139‑1 (MQ=255)                    tttATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATg                                                                                                                                   >  2:255793/1‑118 (MQ=255)                    tttATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATg                                                                                                                                   <  1:255793/118‑1 (MQ=255)                     ttATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTaa                                                                                                            <  1:412964/139‑1 (MQ=255)                        tatTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAAccc                                                                                                         <  2:884315/139‑1 (MQ=255)                                             aTAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAAcc                                                                                                          >  2:1106836/1‑117 (MQ=255)                                             aTAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAAcc                                                                                                          <  1:1106836/117‑1 (MQ=255)                                                   ttCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAg                                                                                                               <  1:85519/106‑1 (MQ=255)                                                   ttCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAg                                                                                                               >  2:85519/1‑106 (MQ=255)                                                   ttCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCAttt                                                                              >  1:1028039/1‑139 (MQ=255)                                                     cAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCAttttt                                                                            >  1:1069872/1‑139 (MQ=255)                                                        aCGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGc                                                                         >  1:923960/1‑139 (MQ=255)                                                          gTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTa                                                                       <  1:992966/139‑1 (MQ=255)                                                           ttCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATg                                                                                                                                   >  1:232169/1‑79 (MQ=255)                                                           ttCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATg                                                                                                                                   <  2:232169/79‑1 (MQ=255)                                                           ttCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTAc                                                                      >  2:91796/1‑139 (MQ=255)                                                             cGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCa                                                                    <  2:739402/139‑1 (MQ=255)                                                              gCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAg                                                                   <  2:1028039/139‑1 (MQ=255)                                                                         ttATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGAt                                                        <  1:521258/139‑1 (MQ=255)                                                                            tatTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCt                                                     >  2:905043/1‑139 (MQ=255)                                                                                gATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCtt                                                 >  1:947028/1‑139 (MQ=255)                                                                                 aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCtgtg                                                                                                                        <  1:177973/67‑1 (MQ=255)                                                                                 aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCtgtg                                                                                                                        >  2:177973/1‑67 (MQ=255)                                                                                 aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGa                                                                                                                <  2:425453/75‑1 (MQ=255)                                                                                 aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGa                                                                                                                >  1:425453/1‑75 (MQ=255)                                                                                 aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTc                                                >  2:842409/1‑139 (MQ=255)                                                                                           gCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATg                                                                                                                                   >  1:487623/1‑47 (MQ=255)                                                                                           gCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATg                                                                                                                                   <  2:487623/47‑1 (MQ=255)                                                                                           gCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAgcgcgc                                      >  2:1376/1‑139 (MQ=255)                                                                                           gCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAgcgcgc                                      >  2:908638/1‑139 (MQ=255)                                                                                           gCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAgcgcgc                                      <  2:923960/139‑1 (MQ=255)                                                                                           gCCGAACAGCGGGCAGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAgcgcgc                                      >  2:854872/1‑139 (MQ=255)                                                                                            ccGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTAcc                                                                     <  1:715814/107‑1 (MQ=255)                                                                                            ccGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTAcc                                                                     >  2:715814/1‑107 (MQ=255)                                                                                               aaCAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACt                                  >  2:250421/1‑139 (MQ=255)                                                                                               aaCAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACt                                  >  2:386031/1‑139 (MQ=255)                                                                                                   gcgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCa                              <  1:908638/139‑1 (MQ=255)                                                                                                    cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCat                                                                                       >  1:1003367/1‑81 (MQ=255)                                                                                                    cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCat                                                                                       <  2:1003367/81‑1 (MQ=255)                                                                                                      ggcggcTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCActc                           >  1:1048954/1‑139 (MQ=255)                                                                                                       gcggcTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGAt                                                        >  2:400572/1‑109 (MQ=255)                                                                                                       gcggcTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGAt                                                        <  1:400572/109‑1 (MQ=255)                                                                                                              tCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAg                                                                                                                     >  2:809904/1‑41 (MQ=255)                                                                                                              tCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAg                                                                                                                     <  1:809904/41‑1 (MQ=255)                                                                                                                      taaaataaaGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACcgtcgt           <  2:1048954/139‑1 (MQ=255)                                                                                                                      taaaataaaGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACcgtcgt           <  1:854872/139‑1 (MQ=255)                                                                                                                      taaaataaaGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACcgtcgt           <  1:386031/139‑1 (MQ=255)                                                                                                                              agtaatCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACCGTCGTCTGtaata   >  2:597432/1‑139 (MQ=255)                                                                                                                               gtaatCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACCGTCGTCTGtaataa  <  2:947028/139‑1 (MQ=255)                                                                                                                                  |                                                                                                                                         AATAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAGTCTCATGCTTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACCGTCGTCTGTAATAA  >  NC_000913/2860287‑2860552

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