Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,963,829	C→T	R134H (CGC→CAC)	gpp ←	guanosine pentaphosphatase/exopolyphosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,963,829	0	C	T	100.0%	161.1 / NA	51	R134H (CGC→CAC)	gpp	guanosine pentaphosphatase/exopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/24);  total (27/24)

GCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGCGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACG  >  NC_000913/3963712‑3963960                                                                                                                      |                                                                                                                                    gCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAgtgg                                                                                                                >  1:839765/1‑139 (MQ=255) gCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAgtgg                                                                                                                >  2:448154/1‑139 (MQ=255) gCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAgtgg                                                                                                                >  2:516552/1‑139 (MQ=255)       tAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCC‑GTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGc                                                                                                         <  1:678267/139‑1 (MQ=255)            ttCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACg                                                                                                     <  1:291055/139‑1 (MQ=255)            ttCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACg                                                                                                     <  2:685119/139‑1 (MQ=255)            ttCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACg                                                                                                     >  2:84324/1‑139 (MQ=255)                               cccATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTg                                                                                  >  2:883661/1‑139 (MQ=255)                                   tCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTc                                                                              <  1:516552/139‑1 (MQ=255)                                       cAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGcac                                                                                                                         <  1:1279104/92‑1 (MQ=255)                                       cAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGcac                                                                                                                         >  2:1279104/1‑92 (MQ=255)                                                caacGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTg                                                                 <  1:992096/139‑1 (MQ=255)                                                caacGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTg                                                                 <  1:448154/139‑1 (MQ=255)                                                caacGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTg                                                                 <  1:883661/139‑1 (MQ=255)                                                caacGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTg                                                                 <  1:84324/139‑1 (MQ=255)                                                   cGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCAc                                                                      <  2:732325/131‑1 (MQ=255)                                                   cGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCAc                                                                      >  1:732325/1‑131 (MQ=255)                                                    gAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCAc                                                                      >  1:1217618/1‑130 (MQ=255)                                                    gAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCAc                                                                      <  2:1217618/130‑1 (MQ=255)                                                    gAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCa                                                             >  2:635237/1‑139 (MQ=255)                                                                    ccGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGc                                                                                                    >  1:1085909/1‑84 (MQ=255)                                                                    ccGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGc                                                                                                    <  2:1085909/84‑1 (MQ=255)                                                                         ccGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGa                                                  <  1:773158/129‑1 (MQ=255)                                                                         ccGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGa                                                  >  2:773158/1‑129 (MQ=255)                                                                         ccGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGaa                                        >  2:1067879/1‑139 (MQ=255)                                                                             ttACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAg                                                                                                                   <  1:123218/60‑1 (MQ=255)                                                                             ttACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAg                                                                                                                   >  2:123218/1‑60 (MQ=255)                                                                                         tACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGc                                                                                                                           >  2:981080/1‑40 (MQ=38)                                                                                         tACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGc                                                                                                                           <  1:981080/40‑1 (MQ=38)                                                                                         tACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCaccacc                                                                                                                     <  1:34087/46‑1 (MQ=255)                                                                                         tACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCaccacc                                                                                                                     <  1:777452/46‑1 (MQ=255)                                                                                         tACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCaccacc                                                                                                                     >  2:777452/1‑46 (MQ=255)                                                                                         tACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCaccacc                                                                                                                     >  2:34087/1‑46 (MQ=255)                                                                                         tACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCAc                                                            <  1:579682/103‑1 (MQ=255)                                                                                         tACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCAc                                                            >  2:579682/1‑103 (MQ=255)                                                                                              gCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGcc                                                                                                   >  1:351668/1‑59 (MQ=255)                                                                                              gCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGcc                                                                                                   <  2:351668/59‑1 (MQ=255)                                                                                                   gCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCg                                                                    >  1:1239262/1‑85 (MQ=255)                                                                                                   gCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCg                                                                    <  2:1239262/85‑1 (MQ=255)                                                                                                   gCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCa              >  1:549512/1‑139 (MQ=255)                                                                                                   gCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCa              >  2:939515/1‑139 (MQ=255)                                                                                                        tatCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGa                                                                                      <  2:263051/62‑1 (MQ=255)                                                                                                        tatCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGa                                                                                      >  1:263051/1‑62 (MQ=255)                                                                                                              ccaccaGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGc                                                                                 >  1:190450/1‑61 (MQ=255)                                                                                                              ccaccaGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGc                                                                                 <  2:190450/61‑1 (MQ=255)                                                                                                               caccaGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGa                                                                >  1:190512/1‑77 (MQ=255)                                                                                                               caccaGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGa                                                                <  2:190512/77‑1 (MQ=255)                                                                                                               caccaGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACg  >  2:795959/1‑139 (MQ=255)                                                                                                               caccaGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTACTGCGCTTTGGCAATAAAATCACCAGCATTGACg  >  1:861308/1‑139 (MQ=255)                                                                                                                  caGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTc                                                                              >  1:245578/1‑60 (MQ=255)                                                                                                                  caGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTc                                                                              <  2:245578/60‑1 (MQ=255)                                                                                                                      |                                                                                                                                    GCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGCGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACG  >  NC_000913/3963712‑3963960

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