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

GATTACGATCGGCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGCGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCA  >  NC_000913/3963701‑3963963                                                                                                                                 |                                                                                                                                       gATTACGATCGGCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGc                                                                                                                              <  1:312500/139‑1 (MQ=255)   ttACGATCGGCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGcac                                                                                                                            <  2:316420/139‑1 (MQ=255)   ttACGATCGGCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGcac                                                                                                                            <  1:616011/139‑1 (MQ=255)        aTCGGCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCAccacca                                                                                                                       <  1:824422/139‑1 (MQ=255)         tCGGCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAg                                                                                                                      <  2:865107/139‑1 (MQ=255)            gCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAgtgg                                                                                                                   >  2:253893/1‑139 (MQ=255)                aaTAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGtgtg                                                                                                               >  1:377980/1‑139 (MQ=255)                     cGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGTTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCAccacca                                                                                                                       <  2:427202/126‑1 (MQ=255)                     cGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAAGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCAccacca                                                                                                                       >  1:427202/1‑126 (MQ=255)                        tCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGc                                                                                                       <  2:210096/139‑1 (MQ=255)                        tCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGc                                                                                                       <  1:993925/139‑1 (MQ=255)                                ggTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATaa                                                                                               >  2:944859/1‑139 (MQ=255)                                 gTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATaaa                                                                                              >  2:351594/1‑139 (MQ=255)                                    aCGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCa                                                                                           >  1:520008/1‑139 (MQ=255)                                          cccATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGc                                                                                                                              >  2:775488/1‑98 (MQ=255)                                          cccATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGc                                                                                                                              <  1:775488/98‑1 (MQ=255)                                            cATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGc                                                                                                       <  1:913390/119‑1 (MQ=255)                                            cATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGc                                                                                                       >  2:913390/1‑119 (MQ=255)                                              tCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTc                                                                                 >  2:943136/1‑139 (MQ=255)                                               cGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACAAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTcc                                                                                >  2:676288/1‑139 (MQ=255)                                                      cTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCAc                                                                         <  1:484420/139‑1 (MQ=255)                                                            aacGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGa                                                                   <  1:676288/139‑1 (MQ=255)                                                            aacGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGa                                                                   <  2:498742/139‑1 (MQ=255)                                                               gAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCa                                                                >  2:102606/1‑139 (MQ=255)                                                                aGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCAc                                                               <  1:253893/139‑1 (MQ=255)                                                                   tggtTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTg                                                            <  1:944859/139‑1 (MQ=255)                                                                                  tGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCGCCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCgag                                             >  2:1049086/1‑139 (MQ=255)                                                                                    ccGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCaa                                                                                                          <  1:465189/76‑1 (MQ=255)                                                                                    ccGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCaa                                                                                                          >  2:465189/1‑76 (MQ=255)                                                                                         tACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTcc                                      >  1:190432/1‑139 (MQ=255)                                                                                             aGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTgcg                                  >  2:488490/1‑139 (MQ=255)                                                                                                      cTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGcac                                                                                                                            >  2:502006/1‑40 (MQ=255)                                                                                                      cTGGCTCCGCCTATATCCACCACCAGGTGCTGATCGGcac                                                                                                                            <  1:502006/40‑1 (MQ=38)                                                                                                              gCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGa                                                                   >  1:883599/1‑89 (MQ=255)                                                                                                              gCCTATATCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGa                                                                   <  2:883599/89‑1 (MQ=255)                                                                                                                     tCCACCACCAGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCAt          <  2:190432/139‑1 (MQ=255)                                                                                                                         ccaccaGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGa                                                     <  2:609527/92‑1 (MQ=255)                                                                                                                         ccaccaGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGa                                                     >  1:609527/1‑92 (MQ=255)                                                                                                                          caccaGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGc                                                                                                       >  2:582657/1‑41 (MQ=255)                                                                                                                          caccaGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGc                                                                                                       <  1:582657/41‑1 (MQ=255)                                                                                                                          caccaGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGc                                                                                    >  2:308054/1‑60 (MQ=255)                                                                                                                          caccaGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGc                                                                                    <  1:308054/60‑1 (MQ=255)                                                                                                                          caccaGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACg     >  2:1033622/1‑139 (MQ=255)                                                                                                                             caGGTGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCa  <  2:520008/139‑1 (MQ=255)                                                                                                                                 |                                                                                                                                       GATTACGATCGGCAAAATAGCGTTCCAGCCAGGTGACGCAGCCCATCGACAGGCTGAACAACGAGGTGGTTTGTGCACCCGTGCCGGTTACCAGTTCAGTACTGGCTCCGCCTATATCCACCACCAGGCGCTGATCGGCACCACCAGTGGTGTGAGCAACGCCCTGATAAATCAGACGTGCTTCCTCTTCACCGCTGATCACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCA  >  NC_000913/3963701‑3963963

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