Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,625,316	T→C	C68R (TGT→CGT)	ppx →	exopolyphosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,625,316	0	T	C	100.0%	188.0 / NA	53	C68R (TGT→CGT)	ppx	exopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/28);  total (25/28)

AGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCC  >  NC_000913/2625195‑2625443                                                                                                                          |                                                                                                                                aGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCg                                                                                                               <  2:73134/140‑1 (MQ=255)    ccGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAAc                                                                                                            <  1:716804/140‑1 (MQ=255)    ccGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAAc                                                                                                            <  1:445804/140‑1 (MQ=255)       tgtgGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGc                                                                                                         >  2:9896/1‑140 (MQ=255)       tgtgGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGc                                                                                                         >  2:429036/1‑140 (MQ=255)       tgtgGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGc                                                                                                         >  1:695227/1‑140 (MQ=255)       tgtgGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGc                                                                                                         >  1:458530/1‑140 (MQ=255)       tgtgGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGc                                                                                                         >  1:420233/1‑140 (MQ=255)           gTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACa                                                                                                     <  2:420233/140‑1 (MQ=255)           gTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACa                                                                                                     <  2:390644/140‑1 (MQ=255)                 ggTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGtt                                                                                               >  2:162483/1‑140 (MQ=255)                   tGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGtttt                                                                                             >  2:327737/1‑140 (MQ=255)                   tGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGtttt                                                                                             >  1:139493/1‑140 (MQ=255)                           aGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGcc                                                                                     <  1:585479/140‑1 (MQ=255)                           aGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGcc                                                                                     <  1:213349/140‑1 (MQ=255)                           aGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGcc                                                                                     <  2:50507/140‑1 (MQ=255)                           aGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGcc                                                                                     <  1:9896/140‑1 (MQ=255)                             attattGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAg                                                                                   >  1:628511/1‑140 (MQ=255)                                 ttGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTc                                                                               <  2:258128/140‑1 (MQ=255)                                    gccgccTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGt                                                                            >  2:32782/1‑140 (MQ=255)                                    gccgccTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGt                                                                            >  1:119572/1‑140 (MQ=255)                                     ccgccTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGta                                                                           <  2:119572/140‑1 (MQ=255)                                       gccTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGcc                                                                                     >  1:303869/1‑128 (MQ=255)                                       gccTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGcc                                                                                     <  2:303869/128‑1 (MQ=255)                                       gccTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGtata                                                                         >  1:416542/1‑140 (MQ=255)                                         cTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGt                                                                       <  2:139493/140‑1 (MQ=255)                                         cTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGt                                                                       >  2:288428/1‑140 (MQ=255)                                           gAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTg                                                                     >  2:632757/1‑140 (MQ=255)                                            aaaCAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTgg                                                                    <  1:327737/140‑1 (MQ=255)                                                 gCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTAccc                                                               <  2:628511/140‑1 (MQ=255)                                                  cGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCa                                                              >  2:135335/1‑140 (MQ=255)                                                               cGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAg                                                 <  2:332109/140‑1 (MQ=255)                                                                   cGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGgcgc                                             >  1:515551/1‑140 (MQ=255)                                                                       cTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGaa                                         <  1:32782/140‑1 (MQ=255)                                                                       cTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGaa                                         <  1:429036/140‑1 (MQ=255)                                                                          gggCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGc                                      <  2:823681/140‑1 (MQ=255)                                                                              cAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACt                                  >  2:456138/1‑140 (MQ=255)                                                                                           gAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAc                     >  2:441861/1‑140 (MQ=255)                                                                                           gAGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAc                     >  2:211338/1‑140 (MQ=255)                                                                                            aGTGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAcg                    >  2:219716/1‑140 (MQ=255)                                                                                              tGAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAcgcg                  <  2:416542/140‑1 (MQ=255)                                                                                               gAAGAGGCAATGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAcgcgc                 <  1:211338/140‑1 (MQ=255)                                                                                                        aTGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGt        <  2:695227/140‑1 (MQ=255)                                                                                                        aTGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGt        <  2:515551/140‑1 (MQ=255)                                                                                                        aTGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACCTTCTGAAACGCGCGGAAAAGGt        <  2:458530/140‑1 (MQ=255)                                                                                                         tGACGCGCGGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTc       <  1:135335/140‑1 (MQ=255)                                                                                                            cgcgcgGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCAtt    <  1:162483/140‑1 (MQ=255)                                                                                                              cgcgGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGaa             <  2:425079/129‑1 (MQ=255)                                                                                                              cgcgGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGaa             >  1:425079/1‑129 (MQ=255)                                                                                                              cgcgGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTcc  >  1:697205/1‑140 (MQ=255)                                                                                                              cgcgGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTcc  >  1:847099/1‑140 (MQ=255)                                                                                                               gcgGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTc       >  1:18262/1‑134 (MQ=255)                                                                                                               gcgGTTTAAACCGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTc       <  2:18262/134‑1 (MQ=255)                                                                                                                          |                                                                                                                                AGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCC  >  NC_000913/2625195‑2625443

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