Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,048,586 A→T 43.9% V18E (GTG→GAG)  glnG ← fused DNA‑binding response regulator in two‑component regulatory system with GlnL: response regulator/sigma54 interaction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,048,5860AT43.9% 1.9 / 22.0 25V18E (GTG→GAG) glnGfused DNA‑binding response regulator in two‑component regulatory system with GlnL: response regulator/sigma54 interaction protein
Reads supporting (aligned to +/- strand):  ref base A (9/5);  new base T (5/6);  total (14/11)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

GCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCAATTTTGCCTGAATGCTGATC  >  CP009273/4048446‑4048730
                                                                                                                                            |                                                                                                                                                
gcgcCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGCGCGCACGTTCAGGCAcccctcgcc                                                                                                                                         >  1:73965/1‑144 (MQ=255)
gcgcCAGCCCGTCCAGTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAGACGTCTGACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGa                                                                                                                                         <  2:302056/150‑1 (MQ=255)
                 cccGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCCCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCTCCCAACGGATGGAACTATCGTCATCg                                                                                                                        >  2:274547/1‑150 (MQ=255)
                  ccGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGTTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCAACGa                                                                                                                       >  1:125726/1‑150 (MQ=255)
                  ccGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGTTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGAATGGAACTATCGTCATCGa                                                                                                                       >  1:146863/1‑150 (MQ=255)
                               atCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCGGTTCTCAAACGTCGGACAGGTTAAACCTCCCCCCGCGACCGCACGTTCAAGCACCCAACGGATGGAACTATCGACATCGACTACCCGCACTAt                                                                                                          >  1:4426/1‑150 (MQ=255)
                                        cagcaCATCCCGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCTCCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTg                                                                                                 <  2:121495/150‑1 (MQ=255)
                                                         tGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCGCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCTCCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTa                                                                                <  2:66542/150‑1 (MQ=255)
                                                           cTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCTCCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTCtt                                                                              >  1:127492/1‑150 (MQ=255)
                                                                  gcgcCTCCAGCACTTCTGCGCCGTTCCCAAACGTCGTACAGGTTAACCCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGAt                                                                       <  2:73965/150‑1 (MQ=255)
                                                                      cTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGc                                                                   <  2:205315/150‑1 (MQ=255)
                                                                         cAGCCCTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACTCCCCCCAGCGAGCGCACGTTCCAGCTCCCAACGGATGGAACTATCGCCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTCCCTGATAGGCAgg                                                                <  1:176465/150‑1 (MQ=255)
                                                                            cACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCTCCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTag                                                             >  1:227304/1‑149 (MQ=255)
                                                                                             cAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATg                                            >  2:51603/1‑150 (MQ=255)
                                                                                                        cAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGCATGCCCTGGCCAAc                                 >  1:262619/1‑150 (MQ=255)
                                                                                                          ggTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTg                               >  2:133704/1‑150 (MQ=255)
                                                                                                          ggTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGGCATCGACTACCAAGACTAGCCCTCGTTGCATAAACGACACCTTGATTTCGTGATAGGCAGGTAAAGCGAGAACTCGGTTGGCACTAGCCAAccg                               >  1:321413/1‑148 (MQ=255)
                                                                                                               aaCCTGCCCCAGCGAGCGCACGTTCAAGCTCCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTaca                          >  1:292309/1‑148 (MQ=255)
                                                                                                               aaCCGGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTaaa                          <  2:212985/150‑1 (MQ=255)
                                                                                                                       ccAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCAAttt                  >  2:1851/1‑150 (MQ=255)
                                                                                                                        cAGCGAGCGCACGTTCAAGCTCCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGGAAATTCAAttct                 >  1:238918/1‑148 (MQ=255)
                                                                                                                           cggccgcACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCAATTTTGcc              <  2:65780/146‑1 (MQ=255)
                                                                                                                                 cACGTTCAAGCTCCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCAATTTTGCCTGAATg        <  2:63188/150‑1 (MQ=255)
                                                                                                                                     ttCAAGCTCCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCAATTTTGCCTGAATGCTGa    <  2:356269/150‑1 (MQ=255)
                                                                                                                                       cAAGCTCCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCAATTTTGCCTGAATGCTGatc  <  2:183821/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                
GCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCAATTTTGCCTGAATGCTGATC  >  CP009273/4048446‑4048730

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: