Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,796,612 C→T 22.6% intergenic (‑170/+48) 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,6120CT22.6% 64.1 / 15.1 37intergenic (‑170/+48)pvdH/PP_4224diaminobutyrate‑2‑oxoglutarate transaminase/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base C (14/14);  new base T (3/5);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01

TGGAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTG  >  NC_002947/4796469‑4796758
                                                                                                                                               |                                                                                                                                                   
tGGAAGCGGTGGACATGGGTAAATCCTCGCGAGCATGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAg                                                                                                                                                >  8:117256/1‑149 (MQ=255)
tGGAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAg                                                                                                                                                >  4:80684/1‑149 (MQ=255)
       ggTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACTCTTGAGTGGGGGGAGGATTTAGCGGGTGTTACGCGATTTTATTTTTGTCTCTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTAcgggcggg                                                                                                                                         <  6:277452/149‑1 (MQ=38)
       ggTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCggg                                                                                                                                         <  6:29314/149‑1 (MQ=255)
           gACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCa                                                                                                                                     <  1:75325/149‑1 (MQ=255)
              aTGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGCCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGCGGGTTTAccc                                                                                                                                 <  1:141777/149‑1 (MQ=37)
                 ggTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTTAACCCGCCCGTAGGAGCGGGTTCACCcgc                                                                                                                               >  1:248203/1‑148 (MQ=255)
                     aaTCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                          <  8:51237/149‑1 (MQ=255)
                      aTCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGCGGGTTTACCCGCGAagag                                                                                                                         <  3:200030/149‑1 (MQ=37)
                      aTCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAagag                                                                                                                         <  7:23895/149‑1 (MQ=255)
                         cTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGc                                                                                                                       <  3:80684/148‑1 (MQ=255)
                            gcgAGCACGGGGTGGGGCGGTTATCGCTTTGTAACCGCTTGAGTGGGGGGAGTATTTAGGGGTTGTTACGCGATTTTTTTTTTGTCTGTGCTGGCCTCTTCGCGGTTAAACCCGCCCTTCGCAGCGGGTTCACCCGCTAAGAGGCCa                                                                                                                     <  2:248203/147‑1 (MQ=255)
                                  aaGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGc                                                                                                             <  3:38729/149‑1 (MQ=255)
                                               gTTATCGCTTTGCAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAACCCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAAt                                                                                                <  2:257344/149‑1 (MQ=255)
                                                       tttGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGCCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGCGGGTTTACCcgc                                                                                                                               <  6:180265/110‑1 (MQ=21)
                                                       tttGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGCCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGCGGGTTTACCcgc                                                                                                                               >  5:180265/1‑110 (MQ=21)
                                                             aaCTCTTGAGTGGGGGTAGTATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACgcgc                                                                                   <  2:188772/148‑1 (MQ=255)
                                                                 cTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCGCGTAGGAGCGGGTTCGCCCGCGATGAGGGCGGCACGGGCGAGGTAGATCGTTGCCCTGGCCGCGCCggga                                                                              >  4:192672/1‑147 (MQ=37)
                                                                           gggAGGATTTAGCGGTTGTTGCGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGcaaca                                                                     >  6:233367/1‑148 (MQ=255)
                                                                           gggAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGcaaca                                                                     >  6:83552/1‑148 (MQ=255)
                                                                            ggaggaTTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTTGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACaaa                                                                   <  1:85591/149‑1 (MQ=255)
                                                                                  tttAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCg                                                              >  2:188349/1‑148 (MQ=255)
                                                                                         gttgttACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCAGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCg                                                       >  1:74776/1‑148 (MQ=255)
                                                                                               aCGCGATTTTATTTTTGTCTGTTCTGGCCTCTTAGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACACGCGAAGAGGCCAGAACAGGCGATAAACATCAATGCCCAGGACGCCCCGGTCCCTGCAACAAAACACGCAAGc                                                         >  8:183562/1‑140 (MQ=38)
                                                                                               aCGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGc                                                         <  7:183562/140‑1 (MQ=255)
                                                                                                      tttattttTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGt                                          >  7:85384/1‑148 (MQ=255)
                                                                                                                 cTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCAGCCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGGCACTACGGCTGTCAAAGCGCAgg                               >  1:59417/1‑149 (MQ=255)
                                                                                                                   gtgCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTAc                                                                             <  5:226070/100‑1 (MQ=255)
                                                                                                                   gtgCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAACAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTAc                                                                             >  6:226070/1‑100 (MQ=255)
                                                                                                                             ctTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCaaa                                      >  7:234201/1‑129 (MQ=255)
                                                                                                                             ctTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCaaa                                      <  8:234201/129‑1 (MQ=255)
                                                                                                                                    ggTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCAc           >  1:112735/1‑149 (MQ=255)
                                                                                                                                     gTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCAc           >  5:104290/1‑148 (MQ=255)
                                                                                                                                      tgaaCCCGCTCCTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGc                              <  6:120502/126‑1 (MQ=37)
                                                                                                                                      tgaaCCCGCTCCTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGc                              >  5:120502/3‑128 (MQ=37)
                                                                                                                                       gaaCCCGCTCCTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGAt        >  3:208513/2‑149 (MQ=38)
                                                                                                                                        aaCCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATc       <  2:74776/149‑1 (MQ=255)
                                                                                                                                          cccgctccTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAAAGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGc     >  5:97148/1‑149 (MQ=38)
                                                                                                                                              cccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGAGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTAACTCCGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTg  >  1:207096/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TGGAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTG  >  NC_002947/4796469‑4796758

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: