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%	174.6 / NA	48	L97P (CTG→CCG)	phnC	phosphonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/22);  total (26/22)

TGCCAACGCGGGTCAGCGCCTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCAGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCCGGTAATCAA  >  NC_000913/4324742‑4325009                                                                                                                                       |                                                                                                                                      tGCCAACGCGGGTCAGCGCCTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGtt                                                                                                                                   >  1:1140904/1‑139 (MQ=255)  gCCAACGCGGGTCAGCGCCTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCggg                                                                                                                                     >  2:841211/1‑136 (MQ=255)    cAACGCGGGTCAGCGCCTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGaa                                                                                                                                <  2:247378/139‑1 (MQ=255)           ggTCAGCGCCTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTgg                                                                                                                         >  2:659397/1‑139 (MQ=255)           ggTCAGCGCCTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTgg                                                                                                                         >  2:1246393/1‑139 (MQ=255)             tCAGCGCCTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAAttg                                                                                                                             >  2:616575/1‑133 (MQ=255)             tCAGCGCCTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAAttg                                                                                                                             <  1:616575/133‑1 (MQ=255)                   ccTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCCGAACGGCGTGCTGCCGAGCGCGCCAATCATCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTa                                                                                                                 <  2:1140904/139‑1 (MQ=255)                   ccTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTa                                                                                                                 <  1:659397/139‑1 (MQ=255)                                             ccGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGa                                                                                       >  2:160894/1‑139 (MQ=255)                                             ccGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACATTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGa                                                                                       >  1:1029222/1‑139 (MQ=255)                                                                   tGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTa                                                                                                                 <  1:1050709/91‑1 (MQ=255)                                                                   tGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTa                                                                                                                 >  2:1050709/1‑91 (MQ=255)                                                                           aaCGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGa                                                         <  2:1029222/139‑1 (MQ=255)                                                                                     tgcCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGAGCCAGGCGGCCTTCGCGCTGGACTGTGCGGcc                                               >  2:522337/1‑139 (MQ=255)                                                                                                 cAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGat                                   >  2:1633722/1‑139 (MQ=255)                                                                                                  aaTCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGata                                  <  1:160894/139‑1 (MQ=255)                                                                                                      aGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGa                                                         <  1:1418638/112‑1 (MQ=255)                                                                                                      aGCACGTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGa                                                         >  2:1418638/1‑112 (MQ=255)                                                                                                           gTTCTCCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCg                         <  1:1633722/139‑1 (MQ=255)                                                                                                             tctcCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGa                                                                                                                        <  1:1376331/42‑1 (MQ=255)                                                                                                             tctcCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGa                                                                                                                        >  2:1376331/1‑42 (MQ=255)                                                                                                             tctcCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGa                                                         <  1:538486/105‑1 (MQ=255)                                                                                                             tctcCAGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGa                                                         >  2:538486/1‑105 (MQ=255)                                                                                                                  aGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAggc                                                                        <  2:611622/85‑1 (MQ=255)                                                                                                                  aGTACGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAggc                                                                        >  1:611622/1‑85 (MQ=255)                                                                                                                      cGCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATc              >  1:1202227/1‑139 (MQ=255)                                                                                                                       gCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGc                                                                             >  1:487180/1‑75 (MQ=255)                                                                                                                       gCTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGc                                                                             <  2:487180/75‑1 (MQ=255)                                                                                                                        cTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCgg                                                                                                 <  1:556011/54‑1 (MQ=255)                                                                                                                        cTCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCgg                                                                                                 >  2:556011/1‑54 (MQ=255)                                                                                                                         tCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCgg                                                                                                 >  2:1322652/1‑53 (MQ=255)                                                                                                                         tCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCgg                                                                                                 <  1:1322652/53‑1 (MQ=255)                                                                                                                         tCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGGTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCcagc                                           >  1:1084364/1‑107 (MQ=255)                                                                                                                         tCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGGTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCcagc                                           <  2:1084364/107‑1 (MQ=255)                                                                                                                         tCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTcgcg                                                              <  2:1440808/88‑1 (MQ=255)                                                                                                                         tCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTcgcg                                                              >  1:1440808/1‑88 (MQ=255)                                                                                                                         tCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGc                              <  2:1320056/120‑1 (MQ=255)                                                                                                                         tCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGc                              >  1:1320056/1‑120 (MQ=255)                                                                                                                         tCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGcc           >  1:734897/1‑139 (MQ=255)                                                                                                                         tCAGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGcc           >  2:166662/1‑139 (MQ=255)                                                                                                                           aGGCGGTTCACCGGGTTGAATTGTTGGAATATGTa                                                                                                                 >  2:1706733/1‑35 (MQ=255)                                                                                                                           aGGCGGTTCACCGGGTTGAATTGTTGGAATATGTa                                                                                                                 <  1:1706733/35‑1 (MQ=255)                                                                                                                           aGGCGGTTCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCCgg         >  1:878572/1‑139 (MQ=255)                                                                                                                                 ttCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGcgcg                                                                                                  <  1:811154/44‑1 (MQ=255)                                                                                                                                 ttCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGcgcg                                                                                                  >  2:811154/1‑44 (MQ=255)                                                                                                                                  tCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGcc           <  2:1381502/130‑1 (MQ=255)                                                                                                                                  tCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGcc           >  1:1381502/1‑130 (MQ=255)                                                                                                                                  tCACCGGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCCGGTAATCaa  <  2:878572/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                      TGCCAACGCGGGTCAGCGCCTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCAGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCCGGTAATCAA  >  NC_000913/4324742‑4325009

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