Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,796,581 T→C 11.7% intergenic (‑139/+79) 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,5810TC11.7% 82.4 / 3.3 34intergenic (‑139/+79)pvdH/PP_4224diaminobutyrate‑2‑oxoglutarate transaminase/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (13/17);  new base C (2/2);  total (15/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 = 7.54e-01

GTTGCGGCACGCTGGCCGGTACAACCCGGGCAAGGCTGGAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGT  >  NC_002947/4796433‑4796718
                                                                                                                                                    |                                                                                                                                         
gTTGCGGCACGCTGGCCGGTACAACCCGGGCAAGGCTGGAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGt                                                                                                                                           <  5:141273/149‑1 (MQ=255)
gTTGCGGCACGCTGGCCGGTACAACCCGGGCAAGGCTGGAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTAGAGTGGGGGGAGGATTTAGCGGTTGTTGCGCGATTTTATTTTTGt                                                                                                                                           <  4:51344/149‑1 (MQ=255)
              gCCGGTACAACCCGGGCAAGGCTGGAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCctct                                                                                                                             >  2:85591/1‑149 (MQ=255)
                 ggTACAACCCGGGCAAGGCTGGAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTcg                                                                                                                          <  5:8911/149‑1 (MQ=255)
                  gTACAACCCGGGCAAGGCTGGAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTcgc                                                                                                                         >  4:38729/1‑149 (MQ=255)
                          cGGGCAAGGCTGGAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAAc                                                                                                                 <  8:141829/149‑1 (MQ=255)
                            ggCAAGGCTGGAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAAccc                                                                                                               <  3:135132/149‑1 (MQ=255)
                                    tGGAAGCGGTGGACATGGGTAAATCCTCGCGAGCATGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAg                                                                                                       >  8:117256/1‑149 (MQ=255)
                                    tGGAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAg                                                                                                       >  4:80684/1‑149 (MQ=255)
                                    tGGAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGCCTGTGCTGGCCTCTTCGCGGGTGAAcccgctcctac                                                                                                       >  2:141777/1‑143 (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)
                                                                gcgAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCct                                                                                                                               <  2:37233/97‑1 (MQ=255)
                                                                gcgAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCct                                                                                                                               >  1:37233/1‑97 (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)
                                                                                                                                                    |                                                                                                                                         
GTTGCGGCACGCTGGCCGGTACAACCCGGGCAAGGCTGGAAGCGGTGGACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAAACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGT  >  NC_002947/4796433‑4796718

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: