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_0009134,317,4090CT96.9% 96.9 / ‑6.7 32G163S (GGC→AGC) phnLribophosphonate triphosphate synthase subunit, putative ABC transporter‑related ATPase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (14/17);  minor base A (1/0);  total (15/17)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TCATGGAAGATGCCTACGATGGCTGCGCCACGGGTTTTGGCTTCGCGAATCAGTTCCACCACCGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCCGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGCACGTTCAGGCGGGTGAGAAGACGCGCGGCTTTAGCGGCGCAGGCTTCACGCGGAACGCCG  >  NC_000913/4317273‑4317542
                                                                                                                                        |                                                                                                                                     
tCATGGAAGATGCCTACAATGGCTGCGCCACGGGTTTTGGCTTCGCGAATCAGTTCCACCACCGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCtgc                                                                                                                                     >  2:539629/1‑139 (MQ=255)
                cGATGGCTGCGCCACGGGTTTTGGCTTCGCGAATCAGTTCCACCACCGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGAcgcg                                                                                                                     >  2:403188/1‑139 (MQ=255)
                cGATGGCTGCGCCACGGGTTTTGGCTTCGCGAATCAGTTCCACCACCGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGAcgcg                                                                                                                     <  2:1017290/139‑1 (MQ=255)
                   tGGCTGCGCCACGGGTTTTGGCTTCGCGAATCAGTTCCACCACCGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCgctg                                                                                                                  >  1:373454/1‑139 (MQ=255)
                      cTGCGCCACGGGTTTTGGCTTCGCGAATCAGTTCCACCACCGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGctgctg                                                                                                               >  1:1288751/1‑139 (MQ=255)
                           ccACGGGTTTTGGCTTCGCGAATCAGTTCCACCACCGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGc                                                                                                          <  1:1121561/139‑1 (MQ=255)
                             aCGGGTTTTGGCTTCGCGAATCAGTTCCACCACCGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCa                                                                                                        <  1:979714/139‑1 (MQ=255)
                                ggTTTTGGCTTCGCGAATCAGTTCCACCACCGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCg                                                                                                     >  1:1189447/1‑139 (MQ=255)
                                   tttGGCTTCGCGAATCAGTTCCACCACCGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGaa                                                                                                  <  1:1299862/139‑1 (MQ=255)
                                                        caccacCGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATggtgccaggtgc                                                                             <  1:403188/139‑1 (MQ=255)
                                                          ccaccGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATg                                                                                              <  1:386777/120‑1 (MQ=255)
                                                          ccaccGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATg                                                                                              >  2:386777/1‑120 (MQ=255)
                                                             ccGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTg                                                                                                                          <  1:843826/89‑1 (MQ=255)
                                                             ccGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCAGCGGGCGATGTTg                                                                                                                          >  2:843826/1‑89 (MQ=255)
                                                              cgcgGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAgg                                                                       <  2:1288751/139‑1 (MQ=255)
                                                               gcggcgCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGt                                                                                                                            >  1:1156898/1‑85 (MQ=255)
                                                               gcggcgCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGt                                                                                                                            <  2:1156898/85‑1 (MQ=255)
                                                                 ggcgCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGt                                                                    <  2:496569/139‑1 (MQ=255)
                                                                       gTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGc                                                              <  1:1295125/139‑1 (MQ=255)
                                                                                  ccAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGCACGTTCAGGCg                                                   >  1:1184230/1‑139 (MQ=255)
                                                                                        aGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGCACGTTCAGGCGGGTgag                                             >  2:708432/1‑139 (MQ=255)
                                                                                                     gTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTg                                                                                                                          <  1:245981/49‑1 (MQ=255)
                                                                                                     gTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTg                                                                                                                          >  2:245981/1‑49 (MQ=255)
                                                                                                         agcagcAGAATGGGGTAGTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGCACGTTCAGGCGGGTGAGAAGACGCGCGGCTTTAg                            <  2:1184230/139‑1 (MQ=255)
                                                                                                                          gTCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGCACGTTCAGGCGGGTGAGAAGACGCGCGGCTTTAGCGGCGCAGGCTTCAcgc           >  1:1152841/1‑139 (MQ=255)
                                                                                                                           tCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGCACGTTCAGGCGGGTGAGAAGAcgcgc                                    <  1:1148646/113‑1 (MQ=255)
                                                                                                                           tCGACGATAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGCACGTTCAGGCGGGTGAGAAGAcgcgc                                    >  2:1148646/1‑113 (MQ=255)
                                                                                                                             gacgaTAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGtt                                                                   >  1:771505/1‑80 (MQ=255)
                                                                                                                             gacgaTAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGATGCCAGGTGCCACAGGCGtt                                                                   <  2:771505/80‑1 (MQ=255)
                                                                                                                               cgaTAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACTGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGCACGTTCAGGCGGGTGAGAAGACGCGCGGCTTTAGCGGCGCAGGCtt                >  1:857706/1‑129 (MQ=255)
                                                                                                                               cgaTAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGCACGTTCAGGCGGGTGAGAAGACGCGCGGCTTTAGCGGCGCAGGCtt                <  2:857706/129‑1 (MQ=255)
                                                                                                                               cgaTAAAGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGCACGTTCAGGCGGGTGAGAAGACGCGCGGCTTTAGCGGCGCAGGCTTCACGCGGAAc      >  1:1226700/1‑139 (MQ=255)
                                                                                                                                   aaaGCTGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGCACGTTCAGGCGGGTGAGAAGACGCGCGGCTTTAGCGGCGCAGGCTTCACGCGGAACGCCg  <  1:670391/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                     
TCATGGAAGATGCCTACGATGGCTGCGCCACGGGTTTTGGCTTCGCGAATCAGTTCCACCACCGCGGCGCTGTTTTTGGCGTCCAGCGAGGCGGTAGGTTCGTCAAGCAGCAGAATGGGGTAGTCGACGATAAAGCCGCGGGCGATGTTGACGCGCTGCTGTTCGCCACCGGAAAATGTCGATGGTGCCAGGTGCCACAGGCGTTCCGGCACGTTCAGGCGGGTGAGAAGACGCGCGGCTTTAGCGGCGCAGGCTTCACGCGGAACGCCG  >  NC_000913/4317273‑4317542

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: