Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,048,586 A→T 100% 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,5860AT65.0% 14.2 / 13.7 20V18E (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 (6/1);  new base T (10/3);  total (16/4)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01
Rejected as polymorphism: Not supported by required number of reads on both strands.

CAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCAATTTTGCCTGAATGC  >  CP009273/4048444‑4048725
                                                                                                                                              |                                                                                                                                           
cagcGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCTCCCAAc                                                                                                                                       <  1:38106/149‑1 (MQ=255)
 agcGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCATGCTCAAGCTcc                                                                                                                                           >  2:289784/1‑144 (MQ=255)
                     cGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAACCGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTc                                                                                                                       >  2:83190/1‑144 (MQ=255)
                               atatCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCTCCCAACGGATGGAACTATCGTCATCGACTACCCAGAc                                                                                                        <  1:274863/149‑1 (MQ=255)
                                    tGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCTCCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATcc                                                                                                   >  1:16654/1‑149 (MQ=255)
                                                aTCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCTAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTg                                                                                            >  2:15254/1‑144 (MQ=255)
                                                   cGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTAAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACg                                                                                    >  1:147408/1‑149 (MQ=255)
                                                    ggCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCACGGCGAGCCCACGCCCAACCACCCAACGGATGGAACGATCACCATCGACTACCCAGACTATCCCTCGCTACATa                                                                                        >  2:126917/1‑144 (MQ=255)
                                                              tCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCTCCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAACCGTCACCt                                                                              >  2:243604/1‑144 (MQ=255)
                                                              tCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCTCCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCt                                                                              >  2:314802/1‑144 (MQ=255)
                                                                  cAGCGCCTCCCGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCt                                                                     <  1:146372/149‑1 (MQ=255)
                                                                             gCACTTCTGCGCCGTTCTCAAACGTGGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTAGCATAAACGTCACCTTTATTTCCTGATAGGCAGGt                                                          >  1:221812/1‑149 (MQ=255)
                                                                                                        tACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCTCCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGcc                               >  1:254412/1‑149 (MQ=255)
                                                                                                        tACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCTCCCAACGGATGGAACTACCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCCGGTCAACCGAGAACTCGGTATGccc                                    >  2:303706/1‑144 (MQ=255)
                                                                                                           aGGTTAATCCTGCCCCAGCGTGCGCACGTTCAAGCTCCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAAc                            <  1:145301/149‑1 (MQ=255)
                                                                                                             gTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCTCCCAACGGATGGAACTATCGTCATCGACTACCCAAACTATCCCTCGTTGCATAAACGCCCCCTTTATTTTCTAATGGGCAGGTAAACCGAGAACTCGGGATGCCCTGGcc                               >  2:37435/1‑144 (MQ=255)
                                                                                                             gTTAAACCTGCCCCAGCGAGCGCACGGTCAAGCTCCCAACGTTTGGCACGATCGGCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCCCCTTTATTTCCTGAGAGGCAGGCAAACCGAGAACCCCATATGCCCTGGcc                               >  2:96872/1‑144 (MQ=255)
                                                                                                                aaaCCTGCCCCAGCGAGCGCACGTTCAAGCTCCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTa                       >  1:2993/1‑149 (MQ=255)
                                                                                                                     tGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCACCTGGTa                       >  2:54171/1‑144 (MQ=255)
                                                                                                                                          aaGCTCCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTGAATTCAATTTTGCCTGAATGc  >  2:130521/1‑144 (MQ=255)
                                                                                                                                              |                                                                                                                                           
CAGCGCCAGCCCGTCCATTCCCGGCATACGGATATCTGAAAGCAGCACATCCGGCGTTTTGCTCGCCAGCGCCTCCAGCACTTCTGCGCCGTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCAATTTTGCCTGAATGC  >  CP009273/4048444‑4048725

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: