Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,357,790 T→C 9.7% intergenic (+92/‑162) PP_3831 → / → PP_3832 hypothetical protein/CsrA‑like carbon storage regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,357,7900TC 9.7% 106.3 / 2.7 41intergenic (+92/‑162)PP_3831/PP_3832hypothetical protein/CsrA‑like carbon storage regulator
Reads supporting (aligned to +/- strand):  ref base T (20/17);  new base C (2/2);  total (22/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

GGAGGAAGTGGCCTTGCTGTTGTTTCTTCAGGCACTTGAGGAACAGGAGGACCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAGAAGGAAGAGGGACAAGCTCCCGCCGCCGCGGCAACTAGCCCAGCGGACTTATA  >  NC_002947/4357641‑4357932
                                                                                                                                                     |                                                                                                                                              
ggaggaAGTGGCCTTGCTGTTGTTTCTTCTGGCACTTGAGGAACAGGAGGACCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGtt                                                                                                                                                <  1:781466/150‑1 (MQ=255)
      aGTGGCCTTGCTGTTGTTTCTTCAGGCACTTGAGGAACAGGAGGACCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTCAACGGt                                                                                                                                          <  2:112858/150‑1 (MQ=255)
            cTTGCTGTTGTTTCTTCAGGCACTTGAGGAACAGGAGGACCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCtt                                                                                                                                    <  1:228127/150‑1 (MQ=255)
              tGCTGTTGTTTCTTCAGGCACTTGAGGAACAGGAGGACCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTg                                                                                                                                  <  2:342763/150‑1 (MQ=255)
                      tttCTTCAGGCACTTGAGGAACAGGAGGACCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGc                                                                                                                          >  2:59287/1‑150 (MQ=255)
                             aGGCACTTGAGGAACAGGAGGACCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATCCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCCCCTTCCACTGaa                                                                                                                   >  2:192775/1‑150 (MQ=255)
                                 aCTTGAGGAACAGGAGGACCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTCAACGGTTTgt                                                                                                                                      >  2:203770/1‑126 (MQ=255)
                                 aCTTGAGGAACAGGAGGACCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTCAACGGTTTgt                                                                                                                                      <  1:203770/127‑2 (MQ=255)
                                   ttGAGGAACAGGAGGACCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACg                                                                                                             >  1:637847/1‑150 (MQ=255)
                                             ggaggaCCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGCAGACACGCAACAGAAtt                                                                                                   >  2:155198/1‑150 (MQ=255)
                                                                 tATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTa                                                                               >  2:771028/1‑150 (MQ=255)
                                                                 tATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTa                                                                               >  2:537902/1‑150 (MQ=255)
                                                                  aTTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCAc                                                                                                                              >  2:801894/1‑102 (MQ=255)
                                                                  aTTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCAc                                                                                                                              <  1:801894/102‑1 (MQ=255)
                                                                  aTTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACg                                                                                                             <  2:413693/119‑1 (MQ=255)
                                                                  aTTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACg                                                                                                             >  1:413693/1‑119 (MQ=255)
                                                                  aTTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTAt                                                                              >  2:629281/1‑150 (MQ=255)
                                                                   ttCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGa                                                                                                                 <  1:256919/114‑1 (MQ=255)
                                                                   ttCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGa                                                                                                                 >  2:256919/1‑114 (MQ=255)
                                                                    tCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATcc                                                                            >  1:156289/1‑150 (MQ=255)
                                                                       gCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGtt                                                                                                                                                >  1:7269/1‑79 (MQ=255)
                                                                       gCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGtt                                                                                                                                                <  2:7269/79‑1 (MQ=255)
                                                                              atgCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGATCACCTTTGCCGACACACCGTTGACGGTTTTCTTTGGCACCTGCAACTAAAGACCCGCAACAGAATTTGTCTGTCGGGGAATTAATATCCTCCGAAAGcc                                                                  >  2:66788/1‑149 (MQ=255)
                                                                                  gaagaaATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCtt                                                              <  1:276094/150‑1 (MQ=255)
                                                                                     gaaATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTaa                                                           <  1:516909/150‑1 (MQ=255)
                                                                                      aaaTTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAAt                                                          <  1:629281/150‑1 (MQ=255)
                                                                                         tttCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATgag                                                       <  2:156289/150‑1 (MQ=255)
                                                                                               aTTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCAc                                                                                                                              >  2:150853/1‑73 (MQ=255)
                                                                                               aTTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCAc                                                                                                                              <  1:150853/73‑1 (MQ=255)
                                                                                                ttAGCCGCCAACTCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTa                                                                                  <  1:624510/116‑1 (MQ=255)
                                                                                                ttAGCCGCCAACTCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTa                                                                                  >  2:624510/1‑116 (MQ=255)
                                                                                                ttAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAgaaggaag                                                >  2:42395/1‑150 (MQ=255)
                                                                                                ttAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAgaaggaag                                                >  1:456245/1‑150 (MQ=255)
                                                                                                  aGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAGAAGGAagag                                              <  1:537902/150‑1 (MQ=255)
                                                                                                        cAACCCCCCTATTCTGGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAGAAGGAAGAGGGACaa                                        <  1:155198/150‑1 (MQ=255)
                                                                                                              cccTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTtgtct                                                                                              >  2:458531/1‑90 (MQ=255)
                                                                                                              cccTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTtgtct                                                                                              <  1:458531/90‑1 (MQ=255)
                                                                                                              cccTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAGAAGGAAGAGGGACAAGCTccc                                  >  2:73960/1‑150 (MQ=255)
                                                                                                              cccTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACACAATTTGTCTGTCGGGGAATTACTATCCGACGCAAACGTCCTTAATCAGGAGGAAGAGGGCCAGGCcccc                                  >  1:252419/1‑150 (MQ=255)
                                                                                                                          cAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAGAAGGAAGAGGGACAAGCTCCCGCCGCCGCGGCa                      >  2:599536/1‑150 (MQ=255)
                                                                                                                              gAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAGAAGGAAGAGGGACAAGCTCCCGCCGCCGCGGCAACTa                  <  2:456245/150‑1 (MQ=255)
                                                                                                                                 cACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAGAAGGAAGAGGGACAAGCTCCCGCCGCCGCGGCAACTAGcc               <  1:59287/150‑1 (MQ=255)
                                                                                                                                          aagcAAACCGTCAACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAGAAGGAAGAGGGACAAGCTCCCGCCGCCGCGGCAACTAGCCCAGCGGACt      >  2:527654/4‑150 (MQ=255)
                                                                                                                                              aaaCCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAGAAGGAAGAGGGACAAGCTCCCGCCGCCGCGGCAACTAGCCCAGCGGACTtata  <  1:42395/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
GGAGGAAGTGGCCTTGCTGTTGTTTCTTCAGGCACTTGAGGAACAGGAGGACCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAGAAGGAAGAGGGACAAGCTCCCGCCGCCGCGGCAACTAGCCCAGCGGACTTATA  >  NC_002947/4357641‑4357932

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