Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,796,604 A→G 11.9% intergenic (‑162/+56) 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,6040AG11.9% 102.9 / 5.9 42intergenic (‑162/+56)pvdH/PP_4224diaminobutyrate‑2‑oxoglutarate transaminase/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base A (18/19);  new base G (3/2);  total (21/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 = 6.90e-01

CGGGCAAGGCTGGAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCG‑CGGG‑TAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCAC  >  NC_002947/4796459‑4796749
                                                                                                                                                   |                                                                                                                                                 
cGGGCAAGGCTGGAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCG‑CGGG‑TAAAc                                                                                                                                                <  5:85626/149‑1 (MQ=255)
            gAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCG‑CGGG‑Taaa                                                                                                                                                 <  3:27581/136‑1 (MQ=255)
            gAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCG‑CGGG‑Taaa                                                                                                                                                 >  4:27581/1‑136 (MQ=255)
                     gACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCG‑CGGG‑TAAACCCGCCCGTAGGAGCGGGTTc                                                                                                                            <  5:266452/148‑1 (MQ=255)
                        aTGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCG‑CGGG‑TAAACCCGCCCGTAGGAGCGGGTTCAccc                                                                                                                        <  5:266095/149‑1 (MQ=255)
                          gggTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCG‑CGGG‑TAAACCCGCCCGTAGGAGCGGGTTCACCcgc                                                                                                                      >  6:233662/1‑149 (MQ=255)
                           ggTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCG‑CGGG‑TAAACCCGCCCGTAGGAGCGGGTTCACCcgc                                                                                                                      >  5:206967/1‑148 (MQ=255)
                            gTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCG‑CGGG‑TAAACCCGCCCGTAGGAGCGGGTTCACCcgcg                                                                                                                     >  7:22193/1‑148 (MQ=255)
                             tAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCG‑CGGG‑TGAACCCGCTCCTACGGGCGGGTTTACCCGCGaa                                                                                                                   <  6:219648/149‑1 (MQ=37)
                             tAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCG‑CGGG‑TAAACCCGCCCGTAGGAGCGGGTTCACCcgcg                                                                                                                     >  4:106243/1‑147 (MQ=255)
                              aaaTCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCG‑CGGG‑TAAACCCGCCCGTAGGAGCGGGTTCACCCGCGaa                                                                                                                   <  3:106243/148‑1 (MQ=255)
                                aTCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCG‑CGGG‑TAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAagag                                                                                                                <  4:55163/149‑1 (MQ=255)
                                   cTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCG‑CGGG‑TAAACCCGCCCGTAGGAGCgg                                                                                                                                >  3:342241/1‑130 (MQ=255)
                                   cTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCG‑CGGG‑TAAACCCGCCCGTAGGAGCgg                                                                                                                                <  4:342241/130‑1 (MQ=255)
                                           cAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCG‑CGGG‑TAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAgg                                                                                                     <  8:22193/149‑1 (MQ=255)
                                              gCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCG‑CGGG‑TAAACCCGCCCGTAGGAGCGGGTTCAccc                                                                                                                        <  8:245433/127‑1 (MQ=255)
                                              gCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCG‑CGGG‑TAAACCCGCCCGTAGGAGCGGGTTCAccc                                                                                                                        >  7:245433/1‑127 (MQ=255)
                                                                     gAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCG‑CGGG‑TAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGAcgcg                                                                           >  5:115301/1‑149 (MQ=255)
                                                                       aaCGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCG‑CGGG‑TAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACgcgc                                                                          <  5:233662/148‑1 (MQ=255)
                                                                        aCGCTTGTGTGGGGGGAGGATTTAGCGGTTGTTAGGCTTTTTTATTTTTGTCTTTGCTGGCCTCTTCG‑CGGG‑TAAACCGGCCGGTAGTCGGGGTTTAACCCGGGATTTGGCAAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCg                                                                        <  6:115301/149‑1 (MQ=255)
                                                                        aCGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCCGGT‑TAAACCCCCCCTTA‑GAGCGGGTTGAACCGCCCAGAGGGCCGCACCCGGGCAAAACAACAAATGCCACGACgcgcgg                                                                        >  6:224940/1‑147 (MQ=255)
                                                                          gCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCG‑CGGG‑TGAACCCGCTCCTAcgggcggg                                                                                                                               >  3:94491/1‑92 (MQ=25)
                                                                          gCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCG‑CGGG‑TGAACCCGCTCCTAcgggcggg                                                                                                                               <  4:94491/92‑1 (MQ=25)
                                                                                       gaggaTTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCG‑CGGG‑TAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAg                                                                                 <  3:242616/125‑1 (MQ=255)
                                                                                       gaggaTTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCG‑CGGG‑TAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAg                                                                                 >  4:242616/1‑125 (MQ=255)
                                                                                         gtaTTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTCCTGGCCTCTTCG‑CGGG‑TAAACCCGCCGTTTGCAGCGGGTTCACCCGCTAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAAcc                                                       <  5:224940/147‑1 (MQ=255)
                                                                                          gaTTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTAG‑CGGGTTAAACCCCCCCG‑AGGAGCGGGGTGACCCGCGAAGAGGGCACCACACGCGATAAAAAAAAAAGCCCCAGACGCGCCGGGAACAGCAACAAAAacc                                                      >  2:161587/1‑149 (MQ=255)
                                                                                                   gttgttACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCG‑CGGG‑TAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGt                                             >  6:334429/1‑149 (MQ=255)
                                                                                                         aCGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCG‑CGGG‑TAAACCCGCCCGGACGAGCGGGTTAACCCGCGAAGAGGCCCGCACAGGCGAAAAAACTCAATGCCCAGGCCGCGCCGGTACCTGCCACAAAACCCCCAAGCCGTCACAAc                                       >  4:203850/1‑149 (MQ=25)
                                                                                                             gATTTTATTTTTGTCTGTGCTGGCCTCTTCG‑CGGG‑TAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGccc                                                                                   <  4:69547/101‑1 (MQ=255)
                                                                                                             gATTTTATTTTTGTCTGTGCTGGCCTCTTCG‑CGGG‑TAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGccc                                                                                   >  3:69547/1‑101 (MQ=255)
                                                                                                                    tttttGTCTGTGCTGGCCTCTTCG‑CGGG‑TAAACCCGCCCGTAGTAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATACACATCAATGCCCa                                                                                  <  4:51163/95‑1 (MQ=255)
                                                                                                                    tttttGTCTGTGCTGGCCTCTTCG‑CGGG‑TAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATACACATCAATGCCCa                                                                                  >  3:51163/1‑95 (MQ=255)
                                                                                                                        tGTCTGTGCTGGCCTCTTCG‑CGGG‑TAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAgcgc                         <  1:161587/148‑1 (MQ=255)
                                                                                                                         gtgtgtgGTGCCCTCTGCG‑CGGT‑AAAACCCCCCCTTCGCCGGGGTTTCACCGGCAAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCa                        <  6:59862/145‑1 (MQ=255)
                                                                                                                                 tGGCCTCTTCG‑CGGG‑TAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCt                    >  7:2564/1‑144 (MQ=255)
                                                                                                                                 tGGCCTCTTCG‑CGGG‑TAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCt                    <  8:2564/144‑1 (MQ=255)
                                                                                                                                     ctctTCG‑CGGG‑TAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTa                                        <  2:312944/120‑1 (MQ=255)
                                                                                                                                     ctctTCG‑CGGG‑TAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTa                                        >  1:312944/1‑120 (MQ=255)
                                                                                                                                     ctctTCG‑CGGG‑TAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGcaca           >  1:107849/1‑149 (MQ=255)
                                                                                                                                           g‑cgGG‑TAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCt      <  2:107849/148‑1 (MQ=255)
                                                                                                                                              ggg‑TGAACCCGCTCCTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCa   >  2:336835/1‑149 (MQ=38)
                                                                                                                                               gg‑TGAACCCGCTCCTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCAc  >  4:331601/1‑149 (MQ=38)
                                                                                                                                               gg‑TAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCAc  >  7:171205/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CGGGCAAGGCTGGAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCG‑CGGG‑TAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCAC  >  NC_002947/4796459‑4796749

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