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

GCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAGTCTCATGCTTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACCGTCGT  >  NC_000913/2860295‑2860543                                                                                                                          |                                                                                                                                gcgcTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCtgtg                                                                                                               >  1:252390/1‑139 (MQ=255)                tatTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGAGTTATATTGATACCGGATCGCCGAACAGCGTGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑GCTTCTGTGAAGCATGAGTAAccc                                                                                                <  2:683159/139‑1 (MQ=255)                 atTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGctcttc                                                                                                                    <  2:606989/119‑5 (MQ=255)                 atTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGctcttc                                                                                                                    >  1:606989/1‑115 (MQ=255)                                     aTAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGc                                                                                 <  1:174499/133‑1 (MQ=255)                                     aTAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGc                                                                                 >  2:174499/1‑133 (MQ=255)                                     aTAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCattat                                                                           >  2:189363/1‑139 (MQ=255)                                         cATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCAt                                                                       >  2:339713/1‑139 (MQ=255)                                           ttCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACACGGTGTTATTGCAGGCCATTATTCAttt                                                                     >  1:643794/1‑139 (MQ=255)                                                aCGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAAtc                                                                                                                               <  1:830725/77‑1 (MQ=255)                                                  gTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTa                                                              <  1:189363/139‑1 (MQ=255)                                                  gTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTa                                                              <  2:252390/139‑1 (MQ=255)                                                   ttCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATATAGTAATCTCATGC‑TCTTCtgt                                                                                                                >  2:254084/1‑88 (MQ=255)                                                   ttCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCtgt                                                                                                                <  1:254084/88‑1 (MQ=255)                                                   ttCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTAc                                                             >  2:902379/1‑139 (MQ=255)                                                     cGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCa                                                           >  2:1045630/1‑139 (MQ=255)                                                                 ttATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGAt                                               <  2:302556/139‑1 (MQ=255)                                                                 ttATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGAt                                               <  1:37242/139‑1 (MQ=255)                                                                      ttGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATg                                                                                                                          <  2:588118/60‑1 (MQ=255)                                                                      ttGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATg                                                                                                                          >  1:588118/1‑60 (MQ=255)                                                                      ttGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGc                                          >  2:475582/1‑139 (MQ=255)                                                                         aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCtgtg                                                                                                               >  1:371990/1‑67 (MQ=255)                                                                         aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCtgtg                                                                                                               <  2:371990/67‑1 (MQ=255)                                                                         aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCtg                                                                                                                 >  1:798697/1‑65 (MQ=255)                                                                         aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCtg                                                                                                                 <  2:798697/65‑1 (MQ=255)                                                                           aCCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCAc                                     >  1:68486/1‑139 (MQ=255)                                                                             cGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCa                                   >  2:315216/1‑139 (MQ=255)                                                                              ggATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCgg                                                                                              <  1:101153/79‑1 (MQ=255)                                                                              ggATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCgg                                                                                              >  2:101153/1‑79 (MQ=255)                                                                              ggATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCa                                                        >  2:400230/1‑117 (MQ=255)                                                                              ggATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCa                                                        <  1:400230/117‑1 (MQ=255)                                                                                aTCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCa                                                        <  1:620311/115‑1 (MQ=255)                                                                                aTCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCACATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCa                                                        >  2:620311/1‑115 (MQ=255)                                                                                     cGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGa                           <  1:339713/139‑1 (MQ=255)                                                                                           gcgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCa                     <  1:1045630/139‑1 (MQ=255)                                                                                           gcgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCa                     <  2:471719/139‑1 (MQ=255)                                                                                            cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCAc                    >  2:660975/1‑139 (MQ=255)                                                                                            cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCAc                    >  1:693902/1‑139 (MQ=255)                                                                                                              taaaataaaGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACcgtcgt  <  2:693902/139‑1 (MQ=255)                                                                                                                          |                                                                                                                                GCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAGTCTCATGCTTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACCGTCGT  >  NC_000913/2860295‑2860543

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