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

CAACAAATAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAGTCTCATGCTTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACCGTCGTCTGTAATA  >  NC_000913/2860282‑2860551                                                                                                                                       |                                                                                                                                        caacaaATAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATATCATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCAGCTGTCATTTTATAAAATAAAGTAAtctc                                                                                                                                     >  1:1171527/1‑139 (MQ=255)    caaATAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATg                                                                                                                                  <  1:1473467/139‑1 (MQ=255)       aTAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGctc                                                                                                                               >  2:1264626/1‑138 (MQ=255)            aggCGCTGATTATTTTATATTTGGAGATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCtg                                                                                                                         <  2:1171527/138‑1 (MQ=255)                      tattTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGa                                                                                                               <  1:257524/139‑1 (MQ=255)                         tttATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAAtc                                                                                                                                       <  2:1503253/113‑1 (MQ=255)                         tttATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAAtc                                                                                                                                       >  1:1503253/1‑113 (MQ=255)                         tttATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTa                                                                                                            <  1:1328405/139‑1 (MQ=255)                            atatTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAAcc                                                                                                         <  1:824753/139‑1 (MQ=255)                             tatTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAAccc                                                                                                        <  1:1264626/139‑1 (MQ=255)                                  ggACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGgtgt                                                                                                   <  2:8395/139‑1 (MQ=255)                                       tatGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGct                                                                                                                                >  2:1255645/1‑106 (MQ=255)                                       tatGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGct                                                                                                                                <  1:1255645/106‑1 (MQ=255)                                              tttGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATg                                                                                                                                  >  1:1516421/1‑97 (MQ=255)                                              tttGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATg                                                                                                                                  <  2:1516421/97‑1 (MQ=255)                                                        ttCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCAttt                                                                             >  2:1614535/1‑139 (MQ=255)                                                                 tCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTa                                                                      >  2:756992/1‑137 (MQ=255)                                                                 tCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTa                                                                      <  1:756992/137‑1 (MQ=255)                                                                 tCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTAcc                                                                    >  2:670684/1‑139 (MQ=255)                                                                  cGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGgtgt                                                                                                   <  2:800995/107‑1 (MQ=255)                                                                  cGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGgtgt                                                                                                   >  1:800995/1‑107 (MQ=255)                                                                  cGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTAcc                                                                    <  2:354963/138‑1 (MQ=255)                                                                  cGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTAcc                                                                    >  1:354963/1‑138 (MQ=255)                                                                              ttATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGAt                                                       <  1:1614535/139‑1 (MQ=255)                                                                                 tatTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAAcc                                                                                                         <  2:780087/86‑1 (MQ=255)                                                                                 tatTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAAcc                                                                                                         >  1:780087/1‑86 (MQ=255)                                                                                 tatTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGAt                                                       >  1:1431749/1‑136 (MQ=255)                                                                                 tatTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGAt                                                       <  2:1431749/136‑1 (MQ=255)                                                                                      aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGAt                                                       <  2:73121/131‑1 (MQ=255)                                                                                      aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGAt                                                       >  1:73121/1‑131 (MQ=255)                                                                                          cGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCtg                                                                                                                         >  1:325253/1‑61 (MQ=255)                                                                                          cGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCtg                                                                                                                         <  2:325253/61‑1 (MQ=255)                                                                                          cGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCa                                           >  1:1499525/1‑139 (MQ=255)                                                                                               cGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATg                                                                                                                                  >  1:519884/1‑48 (MQ=255)                                                                                               cGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATg                                                                                                                                  <  2:519884/48‑1 (MQ=255)                                                                                               cGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGAt                                                       <  1:224417/122‑1 (MQ=255)                                                                                               cGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGAt                                                       >  2:224417/1‑122 (MQ=255)                                                                                                         cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAAcc                                                                                                         >  1:1127656/1‑62 (MQ=255)                                                                                                         cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAAcc                                                                                                         <  2:1127656/62‑1 (MQ=255)                                                                                                         cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGgtgt                                                                                                   >  2:1480866/1‑68 (MQ=255)                                                                                                         cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGgtgt                                                                                                   <  1:1480866/68‑1 (MQ=255)                                                                                                         cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGAc                                  <  1:1003976/133‑1 (MQ=255)                                                                                                         cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCAc                            >  1:1533532/1‑139 (MQ=255)                                                                                                         cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCAc                            >  2:120999/1‑139 (MQ=255)                                                                                                         cgcgcggcTGTAATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCCTGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGAc                                  >  2:1003976/4‑133 (MQ=255)                                                                                                           ggcggcTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACTCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCActc                          >  1:1698608/1‑139 (MQ=255)                                                                                                           ggcggcTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCActc                          >  2:237938/1‑139 (MQ=255)                                                                                                           ggcggcTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCActc                          >  1:254001/1‑139 (MQ=255)                                                                                                                   tCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTgg                               <  2:1477582/126‑1 (MQ=255)                                                                                                                   tCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTgg                               >  1:1477582/1‑126 (MQ=255)                                                                                                                       tttATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACcg              <  1:237938/139‑1 (MQ=255)                                                                                                                           taaaataaaGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACcgtcgt          <  1:120999/139‑1 (MQ=255)                                                                                                                                   agtaatCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACcg              >  2:1312328/1‑127 (MQ=255)                                                                                                                                   agtaatCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACcg              <  1:1312328/127‑1 (MQ=255)                                                                                                                                   agtaatCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACCGTCGTCTGtaata  >  2:1406028/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                        CAACAAATAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAGTCTCATGCTTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACCGTCGTCTGTAATA  >  NC_000913/2860282‑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: ‑