Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,324,876	A→G	L97P (CTG→CCG)	phnC ←	phosphonate ABC transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,324,876	0	A	G	100.0%	125.4 / NA	35	L97P (CTG→CCG)	phnC	phosphonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/19);  total (16/19)

CATGCCAACGCGGGTCAGCGCCTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCAGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCCGGTAATCA  >  NC_000913/4324740‑4325008                                                                                                                                         |                                                                                                                                     cATGCCAACGCGGGTCAGCGCCTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCggg                                                                                                                                    >  1:373391/1‑139 (MQ=255)   tGCCAACGCGGGTCAGCGCCTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGtt                                                                                                                                  >  2:1606058/1‑139 (MQ=255)   tGCCAACGCGGGTCAGCGCCTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGtt                                                                                                                                  >  1:973883/1‑139 (MQ=255)   tGCCAACGCGGGTCAGCGCCTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGtt                                                                                                                                  >  1:731586/1‑139 (MQ=255)    gCCAACGCGGGTCAGCGCCTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTg                                                                                                                                 >  2:1019052/1‑139 (MQ=255)      cAACGCGGGTCAGCGCCTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGaa                                                                                                                               <  1:57346/139‑1 (MQ=255)         cgcgGGTCAGCGCCTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAAttg                                                                                                                            <  2:1490411/139‑1 (MQ=255)         cgcgGGTCAGCGCCTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAAttg                                                                                                                            <  1:1606058/139‑1 (MQ=255)            gggTCAGCGCCTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAAttgttg                                                                                                                         <  1:1466213/139‑1 (MQ=255)             ggTCAGCGCCTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTg                                                                                                                                 <  1:1179453/130‑1 (MQ=255)             ggTCAGCGCCTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTg                                                                                                                                 >  2:1179453/1‑130 (MQ=255)             ggTCAGCGCCTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTgg                                                                                                                        <  1:851167/139‑1 (MQ=255)                                cgcTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATggg                                                                                                     >  2:1015842/1‑139 (MQ=255)                                     tttCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCgg                                                                                                <  2:613830/139‑1 (MQ=255)                                                          gCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTg                                                                                                                                 <  1:984642/85‑1 (MQ=255)                                                          gCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTg                                                                                                                                 >  2:984642/1‑85 (MQ=255)                                                              aaaCAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAggc                                                                       <  2:373391/139‑1 (MQ=255)                                                                    gTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATg                                                                                                                  >  1:1411856/1‑90 (MQ=255)                                                                    gTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATg                                                                                                                  <  2:1411856/90‑1 (MQ=255)                                                                     tGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTc                                                                >  2:1726896/1‑139 (MQ=255)                                                                      gcgcCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCgg                                                                                       <  1:665360/115‑1 (MQ=255)                                                                      gcgcCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCgg                                                                                       >  2:665360/1‑115 (MQ=255)                                                                             aaCGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGa                                                        <  2:460695/139‑1 (MQ=255)                                                                                  cGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACtgtg                                                   >  1:977552/1‑139 (MQ=255)                                                                                  cGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACtgtg                                                   >  1:20875/1‑139 (MQ=255)                                                                                                        aGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGc                             <  2:731586/139‑1 (MQ=255)                                                                                                        aGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGc                             <  1:1019052/139‑1 (MQ=255)                                                                                                            cGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTAcc                         <  1:1015842/139‑1 (MQ=255)                                                                                                                       aCGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTAt              <  1:1726896/139‑1 (MQ=255)                                                                                                                         gCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCg            >  2:1258270/1‑139 (MQ=255)                                                                                                                           tCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCg                                                                                        <  2:1088714/61‑1 (MQ=255)                                                                                                                           tCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCg                                                                                        >  1:1088714/1‑61 (MQ=255)                                                                                                                           tCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGccc                                             >  2:1098745/1‑104 (MQ=255)                                                                                                                           tCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGccc                                             <  1:1098745/104‑1 (MQ=255)                                                                                                                           tCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGcc          <  2:973883/139‑1 (MQ=255)                                                                                                                                   ttCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCCGGTAATCa  >  1:616459/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                     CATGCCAACGCGGGTCAGCGCCTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCAGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCCGGTAATCA  >  NC_000913/4324740‑4325008

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