Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,137,665 A→G 11.6% intergenic (‑64/‑149) PP_4522 ← / → PP_4523 LysR family transcriptional regulator/agmatinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,137,6650AG11.6% 108.0 / 5.2 43intergenic (‑64/‑149)PP_4522/PP_4523LysR family transcriptional regulator/agmatinase
Reads supporting (aligned to +/- strand):  ref base A (19/19);  new base G (2/3);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATCGACCGA  >  NC_002947/5137524‑5137810
                                                                                                                                             |                                                                                                                                                 
ggCGTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTc                                                                                                                                              <  3:230815/147‑1 (MQ=255)
ggCGTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCt                                                                                                                                            <  4:133701/149‑1 (MQ=255)
  cGTTGGCAAAGCCCTGGTGGCGAACCACGCTGATGAAGATGCGAAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTg                                                                                                                                           >  8:159557/1‑148 (MQ=255)
   gTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGaa                                                                                                                                         >  2:137564/1‑149 (MQ=255)
       gCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTaa                                                                                                                                     >  6:144240/1‑149 (MQ=255)
        cAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGGAATTTCTGAACTaaa                                                                                                                                    <  7:62848/149‑2 (MQ=255)
         aaaGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAg                                                                                                                                    >  1:116349/1‑148 (MQ=255)
         aaaGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAg                                                                                                                                    >  6:33873/1‑148 (MQ=255)
         aaaGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAg                                                                                                                                    >  6:134986/1‑148 (MQ=255)
          aaGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTa                                                                                                                                  <  4:40930/149‑1 (MQ=255)
               cTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGt                                                                                                                             <  2:116349/149‑1 (MQ=255)
                  gtggCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGCGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCgg                                                                                                                          <  4:120866/149‑1 (MQ=255)
                           cacGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGAttc                                                                                                                 <  5:134986/149‑1 (MQ=255)
                                     aaGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGgcg                                                                                                       >  8:35514/1‑149 (MQ=255)
                                            gcagcaGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCtt                                                                                                >  2:49783/1‑149 (MQ=255)
                                            gcagcaGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCtt                                                                                                <  5:144240/149‑1 (MQ=255)
                                             cagcagTTTCAGGTCGGGCAAGGGCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTg                                                                                               <  2:1301/149‑1 (MQ=255)
                                              agcagTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGa                                                                                              <  1:137564/149‑1 (MQ=255)
                                                         gTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGc                                                                                   <  7:35514/149‑1 (MQ=255)
                                                                      gAGGCCCTGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTcc                                                                                                             <  6:33784/110‑1 (MQ=255)
                                                                      gAGGCCCTGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTcc                                                                                                             >  5:33784/1‑110 (MQ=255)
                                                                           cATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCccgc                                                                 >  6:127782/1‑149 (MQ=255)
                                                                            aTGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCccgcc                                                                >  2:76003/1‑149 (MQ=255)
                                                                               gTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCgg                                                              >  2:135633/1‑148 (MQ=255)
                                                                                           cacaAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGt                                                  <  1:76003/148‑1 (MQ=255)
                                                                                            acaAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTg                                                 >  4:187137/1‑148 (MQ=255)
                                                                                               aGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATg                                             <  5:33873/149‑1 (MQ=255)
                                                                                                    aTGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGgcgcg                                        >  2:209993/1‑149 (MQ=255)
                                                                                                                tCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAgg                            >  5:71973/1‑149 (MQ=255)
                                                                                                                    aCTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGtt                          <  6:71973/147‑1 (MQ=255)
                                                                                                                        ctctACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGaac                     >  7:30863/1‑148 (MQ=255)
                                                                                                                             cccGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGaacaaa                  <  3:187137/146‑1 (MQ=255)
                                                                                                                              ccGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCCCACGTGAATGGCGCGCGGCGGCACAGGTTCGATCAAacac               >  6:270013/1‑148 (MQ=255)
                                                                                                                               cGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACacaa             <  1:135633/149‑1 (MQ=255)
                                                                                                                                   cttAGTTCAGGAATTTCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCg                        <  2:265113/133‑1 (MQ=255)
                                                                                                                                   cttAGTTCAGGAATTTCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCg                        >  1:265113/2‑134 (MQ=255)
                                                                                                                                      aGTTCAGGAATTTCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGaac                     <  8:184359/134‑1 (MQ=255)
                                                                                                                                      aGTTCAGGAATTTCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTACCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGaac                     >  7:184359/1‑134 (MQ=255)
                                                                                                                                       gTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATcga      <  8:30863/148‑1 (MQ=255)
                                                                                                                                           aGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGGAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATcgaccga  <  6:271899/148‑1 (MQ=255)
                                                                                                                                           aGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCaa                                                           <  1:88155/91‑1 (MQ=255)
                                                                                                                                           aGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCaa                                                           >  2:88155/1‑91 (MQ=255)
                                                                                                                                           aGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCACGCCGGCAAACAACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATcgaccga  >  5:271899/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
GGCGTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATCGACCGA  >  NC_002947/5137524‑5137810

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑