Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 2,681,289 C→A D292Y (GAT→TAT)  glrR ← response regulator regulating glmY sRNA in two‑component system with sensor protein GlrK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,681,2890CA97.1% 108.0 / ‑3.9 34D292Y (GAT→TAT) glrRresponse regulator regulating glmY sRNA in two‑component system with sensor protein GlrK
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (22/11);  total (23/11)
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.56e-01

GTCTTCTGTGCGCTCCGCCAGTGCCGGAATTTTCAGGCTGACAACGTTGAGGCGGTAATAGAGATCTTCGCGGAATTCCCCGCGCGCCATCGCTTTTGGCAGATCACGGTGAGTGGCAGAAATAATCCGC  >  CP009273/2681224‑2681353
                                                                 |                                                                
gTCTTCTGTGCGCTCCGCCAGTGCCGGAATTTTCAGGCTGACAACGTTGAGGCGGTAATAGAGATAtt                                                                >  1:3075733/1‑68 (MQ=255)
   ttctGTGCGCTCCGCCAGTGCCGGAATTTTCAGGCTGACAACGTTGAGGCGGTAATAGAGATATTcgc                                                             >  1:2124014/1‑68 (MQ=255)
       gtgCGCTCCGCCAGTGCCGGAATTTTCAGGCTGACAACGTTGAGGCGGTAATAGAGATATTCGCGGaa                                                         <  1:1379031/68‑1 (MQ=255)
       gtgCGCTCCGCCAGTGCCGGAATTTTCAGGCTGACAACGTTGAGGCGGTAATAGAGATATTCGCGGaa                                                         >  1:1686393/1‑68 (MQ=255)
         gcgcTCCGCCAGTGCCGGAATTTTCAGGCTGACAACGTTGAGGCGGTAATAGAGATATTCGCGGAAtt                                                       >  1:2094554/1‑68 (MQ=255)
           gcTCCGCCAGTGCCGGAATTTTCAGGCTGACAACGTTGAGGCGGTAATAGAGATAATCGCGGAATTcc                                                     >  1:4884435/1‑68 (MQ=255)
                       ccGGAATTTTCAGGCTGACAACGTTGAGGCGGTAATAGAGATATTCGCGGAATTCCCCGCGCGCCATc                                         >  1:5997659/1‑68 (MQ=255)
                              tttCAGGCTGACAACGTTGAGGCGGTAATAGAGATATTCGCGGAATTCCCCGCGCGCCATCGCTTTTg                                  <  1:5222580/68‑1 (MQ=255)
                                tCAGGCTGACAACGTTGAGGCGGTAATAGAGATATTCGCGGAATTCCCCGCGCGCCATCGCTTTTGGc                                <  1:953549/68‑1 (MQ=255)
                                   ggCTGACAACGTTGAGGCGGTAATAGAGATATTCGCGGAATTCCCCGCGCGCCATCGCTTTTGGCAGa                             >  1:4733748/1‑68 (MQ=255)
                                    gCTGACAACGTTGAGGCGGTAATAGAGATATTCGCGGAATTCCCCGCGCGCCATCGCTTTTGGCAGAt                            <  1:5124059/68‑1 (MQ=255)
                                      tGACAACGTTGAGGCGGTAATAGAGATATTCGCGGAATTCCCCGCGCGCCATCGCTTTTGGCAGATCa                          <  1:583569/68‑1 (MQ=255)
                                          aaCGTTGAGGCGGTAATAGAGATATTCGCGGAATTCCCCGCGCGCCATCGCTTTTGGCAGATCACGGt                      <  1:2543517/68‑1 (MQ=255)
                                             gTTGAGGCGGTAATAGAGATATTCGCGGAATTCCCGCGCGCCATCGCTTTTGGCAGATCACGGTGAGt                  <  1:2044314/68‑1 (MQ=255)
                                              ttGAGGCGGTAATAGAGATATTCGCGGAATTCCCCGCGCGCCATCGCTTTTGGCAGATCACGGTGAGt                  >  1:4295254/1‑68 (MQ=255)
                                              ttGAGGCGGTAATAGAGATATTCGCGGAATTCCCCGCGCGCCATCGCTTTTCGCAGATCACGGTGAGt                  >  1:2837916/1‑68 (MQ=255)
                                               tGAGGCGGTAATAGAGATATTCGCGGAATTCCCCGCGCGCCATCGCTTTTGGCAGATCACGGTGAGTg                 >  1:812644/1‑68 (MQ=255)
                                                  ggCGGTAATAGAGATATTCGCGGAATTCCCCGCGCGCCATCGCTTTTGGCAGATCACGGTGAGTGGCa              <  1:5178410/68‑1 (MQ=255)
                                                  ggCGGTAATAGAGATATTCGCGGAATTCCCCGCGCGCCATCGCTTTTGGCAGATCACGGTGACTgtct              >  1:771439/1‑65 (MQ=255)
                                                   gCGGTAATAGAGATATTCGCGGAATTCCCCGCGCGCCATCGCTTTTGGCAGATCACGGTGAGTGGCAg             >  1:3922796/1‑68 (MQ=255)
                                                   gCGGTAATAGAGATATTCGCGGAATTCCCCGCGCGCCATCGCTTTTGGCAGATCACGGTGAGTGGCAg             >  1:4628597/1‑68 (MQ=255)
                                                    cGGTAATAGAGATATTCGCGGAATTCCCCGCGCGCGATCGCTTTTGGCAGATCACGGTGAGTGGCAGa            <  1:5517917/68‑1 (MQ=255)
                                                    cGGTAATAGAGATATTCGCGGAAATCCCCGCGCGCCATCGCTTTTGGCAGATCACGGTGAGTGGCAGa            <  1:5223023/68‑1 (MQ=255)
                                                    aggTAATAGAGATATTCGCGGAATTCCCCGCGCGCCATCGCTTTTGGCAGATCACGGTGAGTGGCAGa            <  1:1113383/67‑1 (MQ=255)
                                                          tAGAGATATTCGCGGAATTCCCCGCGCGCCATCGCTTTTGGCAGATCACGGTGAGTGGCAGAaataat      >  1:830072/1‑68 (MQ=255)
                                                             agaTATTCGCGGAATTCCCCGCGCGCCATCGCTTTTGGCAGATCACGGTGAGTGGCTGAAATAATCCg   >  1:5097169/1‑68 (MQ=255)
                                                             agaTATTCGCGGAATTCCCCGCGCGCCATCGCTTTTGGCAGATCACGGTGAGTGGCAGAAATAATCCg   >  1:1808392/1‑68 (MQ=255)
                                                             agaTATTCGCGGAATTCCCCGCGCGCCATCGCTTTTGGCAGATCACGGTGAGTGGCAGAAATAATCCg   >  1:5405506/1‑68 (MQ=255)
                                                             agaTATTCGCGGAATTCCCCGCGCGCCATCGCTTTTGGCAGATCACGGTGAGTGGCAGAAATAATCCg   >  1:1453049/1‑68 (MQ=255)
                                                             agaTATTCGCGGAATTCCCCGCGCGCCATCGCTTTTGGCAGATCACGGTGAGTGGCAGAAATAATCCg   >  1:654702/1‑68 (MQ=255)
                                                             agaTATTCGCGGAATTCCCCGCGCGCCATCGCTTTTGGCAGATCACGGTGAGTGGCAGAAATAATCCg   >  1:2003000/1‑68 (MQ=255)
                                                             agaTATTCGCGGAATTCCCCGCGCGCCATCGCTTTTGGCAGATCACGGTGAGTGGCAGAAATAATCCg   >  1:4237936/1‑68 (MQ=255)
                                                             agaTATTCGCGGAATTCCCCGCGCGCCATCGCTTTTGGCAGATCACGGTGAGTGGCAGAAATAATCCg   >  1:3364261/1‑68 (MQ=255)
                                                              gaTCTTCGCGGAATTCCCCGCGCGCCATCGCTTTTGGCAGATCACGGTGAGTGGCAGAAATAATCCGc  >  1:2209134/1‑68 (MQ=255)
                                                                 |                                                                
GTCTTCTGTGCGCTCCGCCAGTGCCGGAATTTTCAGGCTGACAACGTTGAGGCGGTAATAGAGATCTTCGCGGAATTCCCCGCGCGCCATCGCTTTTGGCAGATCACGGTGAGTGGCAGAAATAATCCGC  >  CP009273/2681224‑2681353

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: