Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,796,627 C→T 19.8% intergenic (‑185/+33) pvdH ← / ← PP_4224 diaminobutyrate‑2‑oxoglutarate transaminase/sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,796,6270CT19.8% 91.7 / 14.8 46intergenic (‑185/+33)pvdH/PP_4224diaminobutyrate‑2‑oxoglutarate transaminase/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base C (19/17);  new base T (5/4);  total (25/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 = 2.02e-01

GACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAG‑‑GCGATAAACAT‑‑‑CAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGCACAATCGCCAGCCC  >  NC_002947/4796480‑4796773
                                                                                                                                                   |                                                                                                                                                       
gACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTc                                                                                                                                                         <  7:67657/148‑1 (MQ=255)
 aCATGGGTAAATCCTCCCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGCGGGTTTAc                                                                                                                                                       <  6:251067/149‑1 (MQ=38)
   aTGGGTATATGCTCGCGAGCAAGCGGATGGGCGGTTTTCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGCGGGTTTAccc                                                                                                                                                     <  8:73390/149‑1 (MQ=37)
      ggTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCcgc                                                                                                                                                   <  4:21401/148‑1 (MQ=255)
      ggTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCcgc                                                                                                                                                   <  4:242600/148‑1 (MQ=255)
        tAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCcgcg                                                                                                                                                  >  8:186667/1‑147 (MQ=255)
        tAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCcgcg                                                                                                                                                  >  6:150099/1‑147 (MQ=255)
        tAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCGTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGaa                                                                                                                                                <  8:351853/149‑1 (MQ=255)
              cTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGc                                                                                                                                           <  5:150099/148‑1 (MQ=255)
              cTCGCGAGCAAGAGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGc                                                                                                                                           <  4:131950/148‑1 (MQ=255)
               tCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGc                                                                                                                                           >  2:122020/1‑147 (MQ=255)
                      cAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAg‑‑g                                                                                                                                <  8:287925/149‑1 (MQ=255)
                                   ggTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAG‑‑Gc                                                                                                                               >  4:83504/1‑137 (MQ=255)
                                   ggTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAG‑‑Gc                                                                                                                               <  3:83504/137‑1 (MQ=255)
                                       aTCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAG‑‑GCGATAAACAT‑‑‑CAATGcc                                                                                                            >  6:284067/1‑149 (MQ=255)
                                              tGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAG‑‑GCGATAAACAT‑‑‑CAATGCCCAGGa                                                                                                       <  2:103185/147‑1 (MQ=255)
                                               ggAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAG‑‑GCGATAAACAT‑‑‑CAATGCCCAGGAcgc                                                                                                    >  4:366899/1‑149 (MQ=255)
                                                   aCGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGAGGGTAAACCCGCCCGTAGGATCGGGTTAACCCGAGAAGACGCCAGCAAAG‑‑GCGATAAACAT‑‑‑CAAAGCCCATGACGCGCCg                                                                                                >  5:85342/1‑149 (MQ=255)
                                                     gCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGCGGGTTTACCCGCGAaga                                                                                                                                              >  4:75250/1‑106 (MQ=18)
                                                     gCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGCGGGTTTACCCGCGAaga                                                                                                                                              <  3:75250/106‑1 (MQ=18)
                                                                  gaggagTTAGCGGTTGTTACGCGATTTAATTTTGGTCTGTGCTGCCCTCTTCGCGGGTAAACCCGCCCGTCGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAG‑‑GCGATAAACAT‑‑‑CAATGCCCAGGACGCGCCGGTACCTGCAACaaaa                                                                                 <  4:191611/149‑1 (MQ=38)
                                                                      aTTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAAACCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAG‑‑GCGa                                                                                                                             >  1:169276/1‑104 (MQ=255)
                                                                                        gATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAG‑‑GCGATAAACAT‑‑‑CAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCt                                                           <  5:284067/149‑1 (MQ=255)
                                                                                           tttattttTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAG‑‑GCGATAAACAT‑‑‑CAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGt                                                         >  2:238468/1‑148 (MQ=255)
                                                                                               tgtttGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAG‑‑GCGATAAACAT‑‑‑CAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCaaa                                                     <  8:265120/146‑1 (MQ=255)
                                                                                                atttGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAG‑‑GCGATAAACCT‑‑‑CAATGCCCAGGACGCGCCGGTACCTGCAACCAAACCCGCAAGCCGTCACTACGGCTGTCAAAgc                                                   >  7:44262/2‑149 (MQ=255)
                                                                                                     tCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAG‑‑GCGATAAACAT‑‑‑CAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAg                                               <  7:186667/148‑1 (MQ=255)
                                                                                                         tgCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAG‑‑GCGATAAACAT‑‑‑CAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGc                                             <  3:366899/146‑1 (MQ=255)
                                                                                                         tgCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAG‑‑GCGATAAACAT‑‑‑CAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTAc                                          >  5:299712/1‑149 (MQ=255)
                                                                                                             ggCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAG‑‑GCGATAAACAT‑‑‑CAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCgg                                       >  8:99483/1‑148 (MQ=255)
                                                                                                                ctctTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGCCGCGATAAACATACCCCTTGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCg                                        <  1:122020/149‑1 (MQ=37)
                                                                                                                 tctTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAG‑‑GCGATAAACAT‑‑‑CAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCacac                                  >  8:194003/1‑149 (MQ=255)
                                                                                                                 tacgggcgGGTTTACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAG‑‑GCGATAAACAT‑‑‑CAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCacac                                  >  7:353813/6‑149 (MQ=255)
                                                                                                                        gggTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAG‑‑GCGATAAACAT‑‑‑CAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCa                           <  5:6722/149‑1 (MQ=255)
                                                                                                                          gtgaaCCCGCTCCTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAG‑‑GCGATAAACAT‑‑‑CAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCAc                          >  5:73550/4‑148 (MQ=38)
                                                                                                                          gTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAG‑‑GCGATAAACAT‑‑‑CAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACAGCTGCAc                          >  6:37386/1‑148 (MQ=255)
                                                                                                                           tgaaCCCGCTCCTACGGGCGGGTTTACCCGCGAAGAGGCCGGCACAG‑‑GCGATAAACAT‑‑‑CAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGa                        >  1:279036/3‑149 (MQ=38)
                                                                                                                               cccgctccTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAG‑‑GCGATAAACAT‑‑‑CAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGc                    >  2:289575/1‑149 (MQ=38)
                                                                                                                               cccgctccTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAG‑‑GCGATAAACAT‑‑‑CAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGc                    >  3:205872/1‑149 (MQ=38)
                                                                                                                                ccgcccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAG‑‑GCGATAAACAT‑‑‑CAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGcc                   >  6:349100/1‑149 (MQ=255)
                                                                                                                                ccgcccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAG‑‑GCGATAAACAT‑‑‑CAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGcc                   >  1:341868/1‑149 (MQ=255)
                                                                                                                                  gcccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAG‑‑GCGATAAACAT‑‑‑CAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCAc                          <  2:249351/140‑1 (MQ=255)
                                                                                                                                  gcccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAG‑‑GCGATAAACAT‑‑‑CAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCAc                          >  1:249351/1‑140 (MQ=255)
                                                                                                                                      ctacgggCGGGTTTACCCGCGAAGAGGCCAGCACAG‑‑GCGATAAACAT‑‑‑CAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGca               <  1:289575/141‑1 (MQ=255)
                                                                                                                                          gAGCGGGTTCACCCGCGAAGAGGCCAGCACAG‑‑GCGATAAACAT‑‑‑CAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCaaa                                                     >  4:294167/1‑105 (MQ=255)
                                                                                                                                          gAGCGGGTTCACCCGCGAAGAGGCCAGCACAG‑‑GCGATAAACAT‑‑‑CAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCaaa                                                     <  3:294167/105‑1 (MQ=255)
                                                                                                                                                  tCACCCGCGAAGAGGCCAGCACAG‑‑GCGATAAACAT‑‑‑CAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGCACAATCgcccgccc  >  4:120167/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                       
GACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAG‑‑GCGATAAACAT‑‑‑CAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGCACAATCGCCAGCCC  >  NC_002947/4796480‑4796773

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