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%	106.9 / NA	34	R134H (CGC→CAC)	gpp	guanosine pentaphosphatase/exopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/16);  total (18/16)

CCCCAGATTACGATCGGCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGCGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCA  >  NC_000913/3963696‑3963963                                                                                                                                      |                                                                                                                                       ccccAGATTACGATCGGCAAAATAGCGTTCCAGCCAGGTGACGCATCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGa                                                                                                                                   <  2:338505/139‑1 (MQ=255)      gATTACGATCGGCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGc                                                                                                                              <  1:247897/139‑1 (MQ=255)           cGATCGGCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGcaccac                                                                                                                         >  2:174901/1‑139 (MQ=255)                 gCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAgtgg                                                                                                                   >  1:1007436/1‑139 (MQ=255)                 gCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAgtgg                                                                                                                   >  1:326356/1‑139 (MQ=255)                            ttCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCAccacca                                                                                                                       >  1:185805/1‑124 (MQ=255)                            ttCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCAccacca                                                                                                                       <  2:185805/124‑1 (MQ=255)                             tCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGc                                                                                                       <  2:594847/139‑1 (MQ=255)                                 gccagGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTg                                                                                                   >  2:214460/1‑139 (MQ=255)                                           gCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGa                                                                                         <  1:1007299/139‑1 (MQ=255)                                               cccATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGTTAAATCAGACGTg                                                                                     >  2:867746/1‑139 (MQ=255)                                                         ggCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGt                                                                                      >  2:817621/1‑128 (MQ=255)                                                         ggCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGt                                                                                      <  1:817621/128‑1 (MQ=255)                                                           cTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCAc                                                                         <  2:1007436/139‑1 (MQ=255)                                                                     aGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCAc                                                               <  1:174901/139‑1 (MQ=255)                                                                                     cGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAgt                                                                                                                     >  2:986919/1‑69 (MQ=255)                                                                                     cGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAgt                                                                                                                     <  1:986919/69‑1 (MQ=255)                                                                                     cGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCg                                               >  2:341506/1‑139 (MQ=255)                                                                                         ccGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGaa                                           >  2:187329/1‑139 (MQ=255)                                                                                                  aGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTgcg                                  >  1:284702/1‑139 (MQ=255)                                                                                                     tCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCtt                               <  1:341506/139‑1 (MQ=255)                                                                                                       aGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAg                                                                                                                      <  2:1116100/50‑1 (MQ=255)                                                                                                       aGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAg                                                                                                                      >  1:1116100/1‑50 (MQ=255)                                                                                                        gTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTgg                            <  2:43245/139‑1 (MQ=255)                                                                                                              gCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTgg                            >  1:6059/1‑133 (MQ=255)                                                                                                              gCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTgg                            <  2:6059/133‑1 (MQ=255)                                                                                                               cTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATaaa                     <  1:187329/139‑1 (MQ=255)                                                                                                                   gCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATaaa                     >  1:140869/1‑135 (MQ=255)                                                                                                                   gCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATaaa                     <  2:140869/135‑1 (MQ=255)                                                                                                                    ccTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCAc                >  2:967654/1‑139 (MQ=255)                                                                                                                              ccaccaGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTcc                                      >  1:683821/1‑107 (MQ=255)                                                                                                                              ccaccaGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTcc                                      <  2:683821/107‑1 (MQ=255)                                                                                                                               caccaGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACg     >  2:866690/1‑139 (MQ=255)                                                                                                                                  caGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTTCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCa  >  1:777378/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                       CCCCAGATTACGATCGGCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGCGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCA  >  NC_000913/3963696‑3963963

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