Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,563,410 C→G 9.9% G775A (GGA→GCA)  dosP ← oxygen‑sensing c‑di‑GMP phosphodiesterase DosP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,563,4100CG 9.9% 101.6 / 3.9 41G775A (GGA→GCA) dosPoxygen‑sensing c‑di‑GMP phosphodiesterase DosP
Reads supporting (aligned to +/- strand):  ref base C (20/16);  new base G (2/2);  total (23/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

ATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTC  >  NC_000913/1563268‑1563546
                                                                                                                                              |                                                                                                                                        
aTCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGGAAAGAAATATc                                                                                                                                          >  1:895619/1‑143 (MQ=255)
aTCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATc                                                                                                                                          >  1:311869/1‑143 (MQ=255)
          ttATGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCATCGGGAGCACGCTGCTCATCCAGCCCGGCATTTCTTCGGCGGGTGGGGGTCGGGAGAGGCAATAGCCCTGCACAAc                                                                                                                                >  2:360823/1‑143 (MQ=255)
           tAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAAcg                                                                                                                               <  1:750238/143‑1 (MQ=255)
           tAAGGGAAAGGGGCACGTTAAATGAAAGCCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAACATCTACACAGAGTAAAATGTAGGAAAAGAAATATCCCTGAATAAcg                                                                                                                               <  1:440414/143‑1 (MQ=255)
               tggcggggtCACGTTAAAGGAAAACCCGGGGGGGCGGGGGAGAGGAATGTGTCAGATTTTCAGCGGTAACACGCTGCGCAGCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGAc                                                                                                                           <  2:150867/136‑1 (MQ=255)
                  cAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGt                                                                                                                        <  1:870929/143‑1 (MQ=255)
                  cAGTGTCACGTGAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGt                                                                                                                        <  1:18843/143‑1 (MQ=255)
                          cGTTAGATGAAAGTGCGGGAGTGCGGGCGAGGGGAGGTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTg                                                                                                                <  2:624767/143‑1 (MQ=255)
                            ttAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATGCCTGAATAACGCGCCAGTGGATCTTGCg                                                                                                              >  1:649345/1‑143 (MQ=255)
                            ttAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATGCCTGAATAACGCGACAGTGGATCTTGCg                                                                                                              >  1:648314/1‑143 (MQ=255)
                                  gAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATACCTGAATAAGGCGACAGTGGAGCTTGCgaaccac                                                                                                        >  2:646930/1‑137 (MQ=255)
                                                 cGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGAAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGACTCAAACAGCGGTGTATCTTGCGTAGCATCTCAAATTGCTCttt                                                                                         <  2:822810/143‑1 (MQ=255)
                                                        gagGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGCCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTGCg                                                                                  >  1:881475/1‑143 (MQ=255)
                                                                 gTCAGAGTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTcc                                                                         <  1:94632/143‑1 (MQ=255)
                                                                  tCAGATTTTCAGCGGTAACACGCTTCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCg                                                                        >  1:440198/1‑143 (MQ=255)
                                                                    aGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAACAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCCTCCGcc                                                                      >  2:471561/1‑143 (MQ=255)
                                                                          tCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGTCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCacgacg                                                                >  1:484946/1‑143 (MQ=255)
                                                                              cGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTa                                                            >  1:331026/1‑143 (MQ=255)
                                                                              cGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCGCTTTGGTTTCGACGCCCTCCGCCACGACGGTTa                                                            >  2:20079/1‑143 (MQ=255)
                                                                              cGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATATTGCGAAGCAGCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTa                                                            >  1:684850/1‑143 (MQ=255)
                                                                                   acacGCTGCTCATCCAGCCTGGACGTTCTGCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTg                                                       <  1:783522/143‑1 (MQ=255)
                                                                                    cacGCTGCTCATCCAGCCTTGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGa                                                      <  2:368947/143‑1 (MQ=255)
                                                                                      cGCTGCTCATCCCGCCTGGAATGTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAgg                                                    <  1:276253/143‑1 (MQ=255)
                                                                                        ccgctCATGCGACCTGGAATGTCGTCGGGGGGTAGGGGGCGGGAAAAGAAATATCCTGGAATAACGCGCTAGGGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTGTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCt                                                  <  1:8572/141‑1 (MQ=255)
                                                                                         tgctCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCtt                                                 >  1:445711/1‑143 (MQ=255)
                                                                                           ctCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCGCTGAATAACGCGGCAGTGGAGCTTGCGTAGCATCTCAAATTGCTCTTCGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCttgg                                               >  2:357788/1‑141 (MQ=255)
                                                                                            tCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAAGAACGCGACAGTGGATCTTGCGTAGCAGCTCAAATTGCTCGTTGGGTTCGACGCCTTCCGCCACGACGGGTAAATTGAGGCTTTGc                                              >  1:116221/1‑143 (MQ=255)
                                                                                                 cAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACCGGGGATCTTGCGTAGCATCTCAAGTTGCTCTTTGGTTTCGGCGCCTTCCGCCACGACGTTTATCTCGAGGCTTTCCCCaag                                         >  1:9926/1‑142 (MQ=255)
                                                                                                  gcccTTGAATTTCTGCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATCACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATg                                        <  2:805960/141‑1 (MQ=255)
                                                                                                     cgggAATTTCTTCGGCGTGTAGGGGGCGAGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCGCAAATTGCTCTTGGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTg                                     <  2:598020/141‑1 (MQ=255)
                                                                                                        gAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGATCTGTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCATTGCCAAATGcgcgta                                  >  2:801878/1‑138 (MQ=255)
                                                                                                               ttcGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCCCGACGGTTAAATTGAGGCTTCGCCCAGTGGTGGGAACGGCtt                           >  2:440919/1‑143 (MQ=255)
                                                                                                               ttcGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCGCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCtt                           >  1:144113/1‑143 (MQ=255)
                                                                                                               ttcGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCAGCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCAACGACGGTTCAATTGAGGCTTTGCCCAATGCTGGCAGTGGCtt                           >  2:365257/1‑143 (MQ=255)
                                                                                                                          gggggCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCa                <  1:748425/143‑1 (MQ=255)
                                                                                                                           tgcgctggAAAAGAAATATCCCTGAATAACGCGACAGTGGATCGTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCaa               <  1:95065/137‑1 (MQ=255)
                                                                                                                             ggCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAAGTGCTCGTTGGTTTCGGCGCCTTCCGCCCCGACGGTGAAATTGGGGGTGTGCCCACTGCTGGTCCAGGCGCCAACTAAGGCagtg             >  2:333244/1‑140 (MQ=255)
                                                                                                                              gCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGCCGCCTTCCGCCACGACGGTTCAATTGAGGCTTTGCCCAATGCTGGACCTGGCTTCAACTAACGCAACGa            >  2:893443/1‑143 (MQ=255)
                                                                                                                                  gAAAAGAAATATGCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCg        <  2:648314/143‑1 (MQ=255)
                                                                                                                                  gAAAAGAAATACGCCTGAATAACGCGACAGTGTATCTTGCGTAGCATCTCAAATTGCTATTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCg        <  2:649345/143‑1 (MQ=255)
                                                                                                                                   aaaaGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCCATGCTGGTAATGGCTTCAAGTAAGGCAAGGCTGCGt       >  2:292386/1‑143 (MQ=255)
                                                                                                                                        aaaTATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTGTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTc  >  1:362716/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                        
ATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTC  >  NC_000913/1563268‑1563546

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