Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,321,246 A→G 10.0% M1T (ATG→ACG) †
*151R (TGA→CGA) 
phnH ←
phnG ←
carbon‑phosphorus lyase core complex subunit PhnH
carbon‑phosphorus lyase core complex subunit PhnG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,321,2460AG10.0% 90.8 / 2.1 40M1T (ATG→ACG) †
*151R (TGA→CGA) 
phnH
phnG
carbon‑phosphorus lyase core complex subunit PhnH
carbon‑phosphorus lyase core complex subunit PhnG
Reads supporting (aligned to +/- strand):  ref base A (17/19);  new base G (2/2);  total (19/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.91e-01

GCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTAACAGGCGACGAAAACTGTGCTGGGCATCCTGCACGGGAAGCATAAAAGCGGTTTCCAGGGTCATGCGTTGTCTCCGCGAACCATCGTAAAGAAGTCGACCCGGCTGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGA  >  NC_000913/4321118‑4321382
                                                                                                                                |                                                                                                                                        
gCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTAACAGGCGACGAAACCTGTGCTGGGCATCCTGCACGGGAAGCATAAAAGCGGTTTCCAGGGTCATGCGTTGTCTCcgc                                                                                                                            <  1:562859/143‑1 (MQ=255)
   gTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTAACAGGCGACGAAAACTGTGCTGGGCATCCTGCACGGGAAGCATAAAAGCGGTTTCCAGGGTCATGCGTTGTCTCCGCGaa                                                                                                                         <  1:126194/143‑1 (MQ=255)
    gggCCAGCCGCGGTTGGGCGGAGGCGGGGCGACAACCACGCCGGGCTCGCTCATGGCCTTTAACGGGCGACGAAAACTGTGCTGGGCATCCTGCACGGGAAGCATAAAAGCGGTTTCCAGGGTCATGCGTTGTCTCCGCGAAc                                                                                                                        <  2:638013/141‑1 (MQ=255)
     tGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTAACAGGCGACGAAAACTGTGCTGGGCATCCTGCACGGGAAGCATAAAAGCGGTTTCCAGGGTCATGCGTTGTCTCCGCGAAcc                                                                                                                       <  2:569770/143‑1 (MQ=255)
                tttGAGCTGATGCAGGGCGACAATGACGCCCGGCTCGCTCATGGCCTCTAACAGGCGGCGAAAACTGTGCTGGGCATCCTGCACGGGAAGCATAGAAGCGGTTTGCAGGGTCGTGCGTTGTCTCCGCGAACCATCGTAaagaa                                                                                                            >  1:308852/1‑143 (MQ=255)
                tttGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTAACAGGCGACGAAAACTGTGCTGGGCATCCTGCACGGGAAGCATAAAAGCGGTTTCCAGGGACGTGCGTTGCCCCCGCGAACCATCGTAaagaa                                                                                                            >  1:310162/1‑143 (MQ=255)
                     gCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTAACAGGCGACGAAAACTGTGCTGGGCATCCTGCACGGGGAGCATAAAAGCGGTTTCCAGGGTCATGCGTTGTCTCCGCGGACCATCGTAAAGAAGTCGa                                                                                                       >  1:44854/1‑143 (MQ=255)
                     gCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTAACAGGCGACGAAAACTGTGCTGGGCATCCTGCACGGGAAGCATAAAAGCGGTTTCCAGGGTCATGCGTTGTCTCCGCGACCCATCGTAAAGAAGTCGa                                                                                                       >  1:278658/1‑143 (MQ=255)
                     gCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTAACAGGCGACGAAAACTGTGCTGGGCATCCTGCACGGGAAGCATAAAAGCGGTTTCCAGGGTCATGCGTTGTCTCCGCGAACCATCGTAAAGAAGTCGa                                                                                                       <  1:116325/143‑1 (MQ=255)
                      cTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTAACAGGCGACGAAAACTGTGCTGGGCATCCTGCACGGGAAGCATAAAAGCGGATTCCAGGGTCATGCGTTGTCTCCGCGACCCATCGCAAAGAAGTCGAc                                                                                                      >  1:277408/1‑143 (MQ=255)
                          tGCAGGGCGACAATCACGCCCGGCTCGCTCATGACACATCACAGTCTGCGTAAACTGTGCTGGGCATCCTGCACGGGCAGCATAACAGCGGTTTCCATGGTCATGCGTTGTCTCCGCGAACCATCGTAAAGAAGTCGACACgg                                                                                                  >  1:155078/1‑143 (MQ=255)
                                         aCGCCCGGCTCGCTCATGGCCTTTAACAGGCGACGAAAACTGTGCTGGGCATCCTGCACGGGAAGCATAAAAGCGGTTTCCAGGGGCATGCGATGTCTCCGCGAACCATCGTAAAGAAGTCGACCCGGCTGGCGTTCACTTCg                                                                                   >  2:384267/1‑143 (MQ=255)
                                           gCCCGGCTCGCACATGGCACTTAACAGGCGACGAAAACTGTGCTGGGCCTCCCGCACGAAAAGCATAAAAGCGGTTTCCAGGGTCATGCGTTGTCTCCGCGAACCATCGTAAAGAAGTCGACCCGGCTGGCGTTCACTTCGGc                                                                                 <  1:183989/143‑1 (MQ=255)
                                             ccGGCTCGCTCATGGCCTTTAACAGGCGACGAAAACTGTGCTGGGCATCCTGCACGGGAAGCATAAAAGCGGTTTCCAGGGTCATGCGTTGTCTCCGCGAACCATCGTAAAGAAGTCGACACGGCTGGCGTTCACTGCGGccc                                                                               >  2:480399/1‑142 (MQ=255)
                                             ccGGCTCGCTCATGGCCTTTAACAGGCGACGAAAACTGTGCGGGGCATCCTGCGCGGGAAGCATAGAAGCGGTTCCCAGGGTCATGCGTTGGCTCCGCCAACCATCGTCAAGCAGGCGACCCGGCTGCCGTTCACTgccgccc                                                                               >  2:661647/1‑142 (MQ=255)
                                                    ggtcAGGGCCTTTAACAGGCGACGAAAACTGTGCTGGGCATCCTGCACGGGAAGCATAAAAGCGGTTTCCAGGGTCATGCGTTGTCTCCGCGAACCATCGTAAAGAAGTCGACCCGGCTGGCGTTCACTTCGGCCTGGCGTGc                                                                        <  2:771381/141‑1 (MQ=255)
                                                    gctcATGGCCTTTAACAGGCGCCGAAAACTGTGCGGGGCAGCCTGCACGGGAAGCAGAAAAGCGGTTTCCAGGGTCATGCGTTGTCTCCGCGAACCATCGTAAAGAAGTCGACCCGGCTGGCGTTCACTTCGGCCTGGCGTGc                                                                        <  1:35464/143‑1 (MQ=255)
                                                              tttAACAGGCGACGAAAACTGTGCTGGGCATCCTGCACGGGAAGCATAAAAGCGGTTTCCAGGGTCATGCGTTGTCTCCGCGAACCATCGTAAAGAAGTCGACCCGGCTGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGt                                                              >  1:312100/1‑143 (MQ=255)
                                                              tttAACAGGCGACGAAAACTGTGCTGGGCATCCTGCACGGGAAGCATAAAAGCGGTTTCCAGGGTCATGCGTTGTCTCCGCGAACCATCGTAAAGAAGTCGACCCGGATGGCGAACACTTCGGCCTGGCGTGCGGCGCTGCgc                                                              >  2:199331/1‑142 (MQ=255)
                                                                    aGGGGACGAAAACGGTGCTGGGCATCGTGCACGGGAAGCATAAAAGCGGTGTCCAGGGTCATGCGTTGTCTCCGCGAACCATCGTAAAGAAGTCGACCCGGCTGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATa                                                        <  1:587446/143‑1 (MQ=255)
                                                                    aGGCGACGAAAACTGTGCTGGGCATCCTGCACGGGAAGCATAAAAGCGGTTTCCAGGGTCATGCGTTGTCTCCGCGAACCATCGTAAAGAAGTCGACCCGGCTGGCGCTACCCTTGGCCTGGCGTGCGGCAATGCGTGCCATa                                                        >  1:768956/1‑143 (MQ=255)
                                                                       cgacgaAAACTGTGCTGGGCATCCTGCACGGGAAGCATAAAAGCGGTTTCCAGGGTCATGCGTTGTCTCCGCGAACCATCGTAAAGAAGTCGACCCGGCTGGCGTTAACTTCGGCCTGGCGTGCGGCAATGCGTGCCAGACgg                                                     >  1:768004/1‑143 (MQ=255)
                                                                         acgaAAACTGCGCTGGGCATCCTGCACGGGAAGCATAAAAGCGGTTCCGAGGGTCATGCGTTGTCTGCGCTAAGCATCGTAAAGAAGTGGACCCGGCTGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTc                                                   <  2:391128/143‑1 (MQ=255)
                                                                            agcaCTGGGCTGGGCATCCTGCACGGGATGCATAATAGCGGTTTCCAGGGTCATGCGTTGTCTCCGCGAACCAGCGGAAAGAAGGCGACCCGGCTGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGc                                                <  1:138808/140‑1 (MQ=255)
                                                                              aaCTGTGCTGGGCATCCTGCACGGGAAGCATAAAAGCGGTTTCCAGGGTCATGCGTTGTCTCCGCGAACCATCGTAAAGAAGTCGACCCGGCTGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGt                                              <  1:657668/143‑1 (MQ=255)
                                                                                gtgtgCTGGGCATCCTGCACGGGAAGCATAAAAGCGGTTTCCAGGGTCATGCGTTGTCTCCGCGAACCAGCGTAAAGAAGTCGACCCGGCTGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTcc                                            <  1:327123/142‑1 (MQ=255)
                                                                                cTGTGCCGGGCATCCTGCACGGGAAGCACGAACGCGGTCGCCAGGGTCATGCGGGGTCGCTGCGAACCATCGTAAACAAGTCGCCCCGGCTGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTcc                                            <  1:384674/143‑1 (MQ=255)
                                                                                 tgtgCTGGGCATGCTGCACGGGAAGCATAAAAGCGGTTTCCGGGGTCATGCGTTGTCGCCGCGAACCATCGTAAAGAAGTCGACCCGGCTGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTCCa                                           <  1:445056/143‑1 (MQ=255)
                                                                                     cTGGGCATCCTGCACGGGAAGCATAAAAGCGGTTTCCAGGGTCATGCGTTGTCTCCGCGAACCATCGTAAAGAAGTCGACCCGGCTGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTCCCGCgg                                       >  2:149887/1‑143 (MQ=255)
                                                                                       gggCGTCCTTCACGGGACGCATACAAGCGGTTTCCAGGGTCATGCGTTGTCTCCGCGAACCATCGTAAAGAAGTCGACCCGGCTGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTCCAGCgggg                                     <  1:458656/143‑1 (MQ=255)
                                                                                       gggCATCCTGCACGGGAAGCATAACAGCGGTTTCCAGGCTCATGCGTTGTCTCCGCGACCCATCGTAAAGAAGTCGCCCCGGCTGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTCCAGCgggg                                     <  1:481169/143‑1 (MQ=255)
                                                                                                    gggAAGCATAAAAGCGGTTTCCAGGGTCATGCGTTGTCTCCGCGAACCATCGTAAAGAAGTCGACCCGGCTGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTc                        >  1:373696/1‑143 (MQ=255)
                                                                                                         gCATAAAAGCGGTTTCCAGGGTCATGCGTTGTCTCCGCGAACCATCGTAAAGAAGTCGACCCGGCTGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGCCCCGCGGGGCCATCAGGGGTGCTGaga                   >  2:31678/1‑141 (MQ=255)
                                                                                                            tAAAAGCGGTTTCCAGGGTCATGCGTTGTCTCCGCGAACCATCGTAAAGAAGTCGACCCGGCTGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGt                >  1:420565/1‑143 (MQ=255)
                                                                                                              aaaGCGGTTTCCAGGGTCATGCGTTGTCTCCGCGAACCATCGTAAAGAAGTCGACCCGGCTGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGGCAGCGTCCAGCGGGGCCATAAGGGTTTCTGATAGGttt              >  1:56423/1‑143 (MQ=255)
                                                                                                                    gTTTCCAGGGTCGTGCGTTGTCTCCGCGAACCATCGTAAAGAAGTCGACCCGGCTGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTAATAAGTTTTGAAAg        <  2:308852/143‑1 (MQ=255)
                                                                                                                    gTTTCCAGGGTCGTGCGTTGTCTCCGCGAACCATCGTAAAGAAGGCGACCCGGCTGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTAATAAGTTTTGAAAg        <  2:310162/143‑1 (MQ=255)
                                                                                                                     tttCCGGGGTCATGTGTTGTCGTCGCTAACCATAGTAAAGAAGTCGACCCGGCTGGCGTTCACTGCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGt       <  1:387531/143‑1 (MQ=255)
                                                                                                                        ccAGGGTCATGCGTTGTCTCCGCGAACCATCGTAAAGAAGTCGACCCGGCTGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGc    <  1:691730/143‑1 (MQ=255)
                                                                                                                          aGGGTCATGCGTTGTCTCCGCGAACCATCGTAAAGAAGTCGACCCGGCTGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCga  >  1:626778/1‑143 (MQ=255)
                                                                                                                                |                                                                                                                                        
GCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTAACAGGCGACGAAAACTGTGCTGGGCATCCTGCACGGGAAGCATAAAAGCGGTTTCCAGGGTCATGCGTTGTCTCCGCGAACCATCGTAAAGAAGTCGACCCGGCTGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGA  >  NC_000913/4321118‑4321382

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: