Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,390,235 T→C 50.7% I44T (ATC→ACC)  hfq → global sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,390,2350TC50.7% ‑4.6 / 15.1 16I44T (ATC→ACC) hfqglobal sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication
Reads supporting (aligned to +/- strand):  ref base T (3/5);  new base C (4/4);  total (7/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.07e-01

GAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGC‑CAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAG  >  CP009273/4390094‑4390346
                                                                                                                                             |                                                                                                                
gAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATcc                                                                                                                >  1:228778/1‑144 (MQ=255)
   aaGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGCCTTTTGATCAGTTCGTGACCCTGt                                                                                                             >  1:249761/1‑144 (MQ=255)
      agagaATGGCTAAGGGGCCAGTTTGACAAGATCCGTTCCTGACCGCACTGCGTCGGGAACGTGTTCCAGGTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGa                                                                                                          <  1:94554/144‑1 (MQ=255)
                            tttACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATGTGGTGAATGGTGTTAAGCTGAAAGGGCAAATCGCGTGTTTTGATCAGTTAGTGACCCTGTTGAAAAACACTGGCAGC‑CAGATGGt                                                                                   >  1:44753/1‑144 (MQ=255)
                            tttACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTAATACAGATTCTAATTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGC‑CAGATGGTTTACa                                                                              >  2:3816/1‑149 (MQ=255)
                             ttACAAGATCCGTTCCTGAACGCACTGCGTCGTGAACGTGTTCCAGTTTCTATATATTTGATGAAAGATATTAAGCTGCAAGGGCAAAACGAGTCTTTTGATCAATTCGTGACCCTGTTGAAAAACACGCGCAGC‑CAGATGGTATACaa                                                                             >  2:123517/1‑149 (MQ=255)
                               aCAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGC‑CAGATGGGTTa                                                                                >  1:6833/1‑144 (MQ=255)
                                            cTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGC‑CAGATGGTTTACAAGCACGCGATTTCTAc                                                              <  2:105793/149‑1 (MQ=255)
                                              gAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAGTGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGC‑CAGATGGTTTACAAGCACGCGCTTTCTACtg                                                            >  2:169265/1‑149 (MQ=255)
                                                     cttcGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGC‑CAGATGGTTTACAAGCACGCGATTTCTACtgtt                                                          <  1:151868/141‑1 (MQ=255)
                                                        cGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGC‑CAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGtc                                                  <  2:231317/149‑1 (MQ=255)
                                                                        cAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGC‑CAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCtcat                                  <  2:188834/149‑1 (MQ=255)
                                                                          gtctCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGC‑CAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCtcatca                                <  2:250534/146‑1 (MQ=255)
                                                                                 tttatttGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTCTGATCAGTTCGTGATCCTGTTATGCTGGGCTGGCTGCTGGGATTGTTT‑CTAGCTCGCGATTTCTATTGTTGTCCCGTCTCTCCCGGTTTCTCATCACAGTaac                         <  2:124312/149‑1 (MQ=255)
                                                                                            aaTGGCATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGCGACCCTGTTGAAAAACACGGTCAGC‑CAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGcc                   <  1:61769/144‑1 (MQ=255)
                                                                                                        cTGCAAGGGCAAATCGAGTCTTTTGATCAGTCCGTGATCCTGTTGAAAAACACGGTCAGC‑CAGATGGTTTACAAGCACGCGATTTCTATTGTTGTCTTGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACcagcag  <  2:7541/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                
GAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGC‑CAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAG  >  CP009273/4390094‑4390346

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