Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,317,409	C→T	G163S (GGC→AGC)	phnL ←	ribophosphonate triphosphate synthase subunit, putative ABC transporter‑related ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,317,409	0	C	T	100.0%	150.4 / NA	47	G163S (GGC→AGC)	phnL	ribophosphonate triphosphate synthase subunit, putative ABC transporter‑related ATPase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/22);  total (25/22)

ATGGAAGATGCCTACGATGGCTGCGCCACGGGTTTTGGCTTCGCGAATCAGTTCCACCACCGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCCGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGCACGTTCAGGCGGGTGAGAAGACGCGCGGCTTTAGCGGCGCAGGCTTCACGCGGAACGC  >  NC_000913/4317275‑4317540                                                                                                                                       |                                                                                                                                    aTGGAAGATGCCTACGATGGCTGCGCCACGGGTTTTGGCTTCGCGAATCAGTTCCACCACCGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCgg                                                                                                                                 <  1:1512243/139‑1 (MQ=255)     aaGATGCCTACGATGGCTGCGCCACGGGTTTTGGCTTCGCGAATCAGTTCCACCACCGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGa                                                                                                                             <  1:868349/139‑1 (MQ=255)          gCCTACGATGGCTGCGCCACGGGTTTTGGCTTCGCGAATCAGTTCCACCACCGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTg                                                                                                                        >  1:766968/1‑139 (MQ=255)               cGATGGCTGCGCCACGGGTTTTGGCTTCGCGAATCAGTTCCACCACCGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCggg                                                                                                                                >  2:1199739/1‑126 (MQ=255)               cGATGGCTGCGCCACGGGTTTTGGCTTCGCGAATCAGTTCCACCACCGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGAcgcg                                                                                                                   >  2:320394/1‑139 (MQ=255)               cGATGGCTGCGCCACGGGTTTTGGCTTCGCGAATCAGTTCCACCACCGCGGCGCTGATCTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCggg                                                                                                                                <  1:1199739/126‑1 (MQ=255)                     cTGCGCCACGGGTTTTGGCTTCGCGAATCAGTTCCACCACCGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCggg                                                                                                                                >  2:609668/1‑120 (MQ=255)                     cTGCGCCACGGGTTTTGGCTTCGCGAATCAGTTCCACCACCGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCggg                                                                                                                                <  1:609668/120‑1 (MQ=255)                     cTGCGCCACGGGTTTTGGCTTCGCGAATCAGTTCCACCACCGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGctgctg                                                                                                             >  2:1522790/1‑139 (MQ=255)                     cTGCGCCACGGGTTTTGGCTTCGCGAATCAGTTCCACCACCGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGctgctg                                                                                                             >  2:1913472/1‑139 (MQ=255)                                     ggCTTCGCGAATCAGTTCCACCACCGCGGCGCTGTTTTGGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAAt                                                                                             >  2:373749/1‑139 (MQ=255)                                         tCGCGAATCAGTTCCACCACCGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCg                                                                                         <  1:1746641/139‑1 (MQ=255)                                          cgcgAATCAGTTCCACCACCGCGGCGCTGTTTTGGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGa                                                                                        <  1:373749/139‑1 (MQ=255)                                                  aGTTCCACCACCGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATggtgcca                                                                                <  1:1522790/139‑1 (MQ=255)                                                    ttCCACCACCGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATggtgcc                                                                                 <  1:58683/136‑1 (MQ=255)                                                    ttCCACCACCGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATggtgcc                                                                                 >  2:58683/1‑136 (MQ=255)                                                         ccaccGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATggtgccaggtgcca                                                                         >  2:1791746/1‑139 (MQ=255)                                                              gcggcgCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTg                                                                                                                        <  1:61725/87‑1 (MQ=255)                                                              gcggcgCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTg                                                                                                                        >  2:61725/1‑87 (MQ=255)                                                                      gTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGc                                                            <  1:1791746/139‑1 (MQ=255)                                                                               gTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGACGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGCACGTTCAgg                                                   <  1:1335940/139‑1 (MQ=255)                                                                                 ccAGCGAGGCGGTAGGTTGGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATggtg                                                                                   >  2:1355910/1‑105 (MQ=255)                                                                                 ccAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATATAGCTGCGGGCGATGCTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATggtg                                                                                   <  1:1355910/105‑1 (MQ=255)                                                                                 ccAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGCACGTTCAgg                                                   <  1:1566546/137‑1 (MQ=255)                                                                                 ccAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGCACGTTCAgg                                                   >  2:1566546/1‑137 (MQ=255)                                                                                    gCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGCACGTTCAGGCGGGt                                              <  1:320394/139‑1 (MQ=255)                                                                                       aGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGtt                                                                                                           <  1:938192/75‑1 (MQ=255)                                                                                       aGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGtt                                                                                                           >  2:938192/1‑75 (MQ=255)                                                                                                  tCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCg                                                                                                                              <  1:1226446/45‑1 (MQ=255)                                                                                                  tCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCg                                                                                                                              >  2:1226446/1‑45 (MQ=255)                                                                                                       aaGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGCACGTTCAGGCGGGTGAGAAGACGCGCGGCTTTa                           >  1:473091/1‑139 (MQ=255)                                                                                                         gcagcaGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGCACGTTCAGGCGGGTGAGAAGACGCGCGGCTTTAgc                         >  1:1245554/1‑139 (MQ=255)                                                                                                           agcagAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGCACGTTCAGGCGGGTGAGAAGACGCGCGGCTTTAgcgg                       >  2:950543/1‑139 (MQ=255)                                                                                                            gcagAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGCACGTTCAGGCGGGTGAGAAGACGCGCGGCTTTAgcggc                      >  1:1689224/1‑139 (MQ=255)                                                                                                              agAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGCACGTTCAGGCGGGTGAGAAGACGCGCGGCTTTAGCGgcgc                    <  1:950543/139‑1 (MQ=255)                                                                                                              agAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGCACGTTCAGGCGGGTGAGAAGACGCGCGGCTTTAGCGgcgc                    <  2:473091/139‑1 (MQ=255)                                                                                                               gAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGCACGTTCAGGCGGGTGAGAAGACGCGCGGCTTTAGCGGCGCa                   >  2:1165440/1‑139 (MQ=255)                                                                                                                     ggTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGCACGTTCAGGCGGGTGAGAAGACGCGCGGCTTTAGCGGCGCAGGCTTc             >  2:1158306/1‑139 (MQ=255)                                                                                                                          tCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCgg                                                             >  1:740098/1‑86 (MQ=255)                                                                                                                          tCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCgg                                                             <  2:740098/86‑1 (MQ=255)                                                                                                                          tCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGACGTTCCGGCACGTTCAGGCGGGTGAGAAGACGCGCGGCTTTAGCGGCGCAGGCTTCAcgcg        >  2:1268233/1‑139 (MQ=255)                                                                                                                            gacgaTAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGCACGTTCAGGCGGGTGAGAAGACGCGCGGCTTTAGCGGCGCAGGCTTCACGCGGa      >  1:1404834/1‑139 (MQ=255)                                                                                                                              cgaTAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATggtgccagg                                                                              <  1:919248/65‑1 (MQ=255)                                                                                                                              cgaTAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATggtgccagg                                                                              >  2:919248/1‑65 (MQ=255)                                                                                                                               gaTAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGCACGTTCAGGCGGGTGAGAAGACGCGCGGCTTTAGCGGCGCAGGCTTCACGCGGAACg   >  2:1313091/1‑139 (MQ=255)                                                                                                                                aTAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGCACGTTCAGGCGGGTGAGAAGACGCGCGGCTTTAGCGGCGCAGGCTTCACGCGGAACGc  <  1:1313091/139‑1 (MQ=255)                                                                                                                                    aGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCAcc                                                                                                    <  2:411773/37‑1 (MQ=255)                                                                                                                                       |                                                                                                                                    ATGGAAGATGCCTACGATGGCTGCGCCACGGGTTTTGGCTTCGCGAATCAGTTCCACCACCGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCCGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGCACGTTCAGGCGGGTGAGAAGACGCGCGGCTTTAGCGGCGCAGGCTTCACGCGGAACGC  >  NC_000913/4317275‑4317540

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: ‑