Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,796,659 G→C 8.5% intergenic (‑172/+46) 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,6590GC 8.5% 136.0 / 5.5 47intergenic (‑172/+46)pvdH/PP_4224diaminobutyrate‑2‑oxoglutarate transaminase/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (23/20);  new base C (2/2);  total (25/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 = 5.66e-01

GGAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGCAC  >  NC_002947/4796515‑4796806
                                                                                                                                                |                                                                                                                                                   
ggAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGa                                                                                                                                                >  1:484601/1‑150 (MQ=255)
    gCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGctg                                                                                                                                            <  2:484601/150‑3 (MQ=255)
     cGGAGGACAGGGTTAAATCCGTGCGAGCAAGCGGAGTGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCggg                                                                                                                                           <  2:331554/150‑1 (MQ=255)
       gTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCGTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGtt                                                                                                                                         >  1:375677/1‑150 (MQ=255)
            cATGGGTAAATCCTCGCGAGCACACGACGCGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGGTTTAccc                                                                                                                                    <  1:965596/150‑1 (MQ=25)
             aTGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCcg                                                                                                                                   >  1:955921/1‑150 (MQ=255)
                    aaTCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCa                                                                                                                                       >  2:587173/1‑139 (MQ=255)
                    aaTCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCa                                                                                                                                       <  1:587173/139‑1 (MQ=255)
                     aTCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAgg                                                                                                                           >  2:325480/1‑150 (MQ=255)
                      tCCTCGCGAGCAAGCGGATGGGCGGTTATCTCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGc                                                                                                                          <  1:325480/150‑1 (MQ=255)
                      tCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGc                                                                                                                          >  2:982727/1‑150 (MQ=255)
                         tCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGa                                                                                                                                                <  2:1286699/125‑1 (MQ=255)
                         tCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGa                                                                                                                                                >  1:1286699/1‑125 (MQ=255)
                           gcgAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCTAAGTGGCCAGca                                                                                                                     <  2:955921/150‑1 (MQ=255)
                                 aaGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGCGGGTTTAc                                                                                                                                      <  2:512573/127‑1 (MQ=25)
                                   gCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                             >  1:1130935/1‑134 (MQ=255)
                                   gCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                             <  2:1130935/134‑1 (MQ=255)
                                           gCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAgg                                                                                                                 <  1:206820/138‑1 (MQ=255)
                                           gCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAgg                                                                                                                 >  2:206820/1‑138 (MQ=255)
                                               ttATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGACTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGc                                                                                                 >  2:153382/1‑150 (MQ=255)
                                                 aTCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTAGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGccc                                                                                               >  1:456354/1‑150 (MQ=255)
                                                           aaaCGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGcc                                                                                     >  1:432811/1‑150 (MQ=255)
                                                           aaaCGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGcc                                                                                     >  1:1089841/1‑150 (MQ=255)
                                                             aCGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCcgc                                                                                                                                  >  1:6368/1‑103 (MQ=255)
                                                             aCGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCcgc                                                                                                                                  <  2:6368/103‑1 (MQ=255)
                                                                    aGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCa                                                                                              <  2:397609/132‑1 (MQ=255)
                                                                    aGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCa                                                                                              >  1:397609/1‑132 (MQ=255)
                                                                                  ttAGCGGTTGTTACGCGATTTTATTTTTGGCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCTATAAACATCAATGCCCAGGACACGCCGGTACCTGCAACAAAACCCGCaa                                                              <  1:623627/150‑1 (MQ=255)
                                                                                    aGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAagag                                                                                                                            <  1:570791/86‑1 (MQ=255)
                                                                                    aGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAagag                                                                                                                            >  2:570791/1‑86 (MQ=255)
                                                                                     gCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGcc                                                           >  2:1153635/1‑150 (MQ=255)
                                                                                        gttgttACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGcaa                                                                          <  2:671865/132‑1 (MQ=255)
                                                                                        gttgttACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGcaa                                                                          >  1:671865/1‑132 (MQ=255)
                                                                                        gttgttACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTc                                                        <  1:1153635/150‑1 (MQ=255)
                                                                                           gttACGCGATTTTATTTTTGTCTGTAAAGGCTTCTTCGCAACGACAATCACCCGTAGGTGAAGGGTCACCTGGGAAGCCGACGATGTCTTCGAGAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGGAAGCCGTCACt                                                     >  1:1254409/1‑150 (MQ=255)
                                                                                            ttACGCGATTTTATTTTTGTCTGTGCTGGACTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTa                                                    <  1:153382/150‑1 (MQ=255)
                                                                                                   attttattttTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGt                                             <  1:982727/150‑1 (MQ=255)
                                                                                                   attttattttTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGt                                             <  1:749751/150‑1 (MQ=255)
                                                                                                           tttGTCTGTGCTGGCCTCTTAGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAgcgc                                     <  2:323984/150‑1 (MQ=255)
                                                                                                                      tGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGc                          >  1:320699/1‑150 (MQ=255)
                                                                                                                           tctTCGCGGGTGAACCCGCTCCTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACAcc                     >  1:418668/1‑150 (MQ=37)
                                                                                                                             ttCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCt                                               <  2:1226129/122‑1 (MQ=255)
                                                                                                                             ttCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCt                                               >  1:1226129/1‑122 (MQ=255)
                                                                                                                                     tAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGAt           >  1:435555/1‑150 (MQ=255)
                                                                                                                                       aaCCCGCTCCTACGGGCGGGTCGACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCg         >  2:1077903/1‑150 (MQ=38)
                                                                                                                                        aCCCGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAga                                                                                                                 <  1:719201/45‑2 (MQ=11)
                                                                                                                                        aCCCGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAga                                                                                                                 >  2:719201/1‑44 (MQ=11)
                                                                                                                                         cccgctccTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTAc                                                   <  1:131607/106‑1 (MQ=37)
                                                                                                                                         cccgcccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGcc                                                                                                <  2:1083818/61‑1 (MQ=255)
                                                                                                                                         cccgcccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGcc                                                                                                >  1:1083818/1‑61 (MQ=255)
                                                                                                                                          ccgcccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACgg                                                 <  1:567327/107‑1 (MQ=255)
                                                                                                                                          ccgcccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACgg                                                 >  2:567327/1‑107 (MQ=255)
                                                                                                                                          ccgcccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCt      >  2:1134498/1‑150 (MQ=255)
                                                                                                                                             cccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGca   <  2:1089841/150‑1 (MQ=255)
                                                                                                                                             cccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGCTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGca   <  2:232701/150‑1 (MQ=255)
                                                                                                                                              ccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCTCAGGCTACCGGCACACCGCTGCACGATCGCCTGcac  >  2:318688/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GGAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGCAC  >  NC_002947/4796515‑4796806

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