Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,796,619 A→G 21.6% intergenic (‑177/+41) 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,6190AG21.6% 63.5 / 14.7 37intergenic (‑177/+41)pvdH/PP_4224diaminobutyrate‑2‑oxoglutarate transaminase/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base A (14/15);  new base G (4/4);  total (18/19)
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.02e-01

GGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGCACAAT  >  NC_002947/4796476‑4796764
                                                                                                                                               |                                                                                                                                                  
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                                            <  8:234201/129‑1 (MQ=255)
                                                                                                                      ctTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCaaa                                            >  7:234201/1‑129 (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)
                                                                                                                                             ggAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCATGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGCACAAt  <  7:117256/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGCACAAT  >  NC_002947/4796476‑4796764

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: