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

GAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAG  >  CP009273/4390094‑4390346
                                                                                                                                             |                                                                                                                
gAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATcc                                                                                                                >  1:228778/1‑144 (MQ=255)
   aaGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGCCTTTTGATCAGTTCGTGACCCTGt                                                                                                             >  1:249761/1‑144 (MQ=255)
      agagaATGGCTAAGGGGCCAGTTTGACAAGATCCGTTCCTGACCGCACTGCGTCGGGAACGTGTTCCAGGTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGa                                                                                                          <  1:94554/144‑1 (MQ=255)
                            tttACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATGTGGTGAATGGTGTTAAGCTGAAAGGGCAAATCGCGTGTTTTGATCAGTTAGTGACCCTGTTGAAAAACACTGGCAGCCAGATGGt                                                                                   >  1:44753/1‑144 (MQ=255)
                            tttACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTAATACAGATTCTAATTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACa                                                                              >  2:3816/1‑149 (MQ=255)
                             ttACAAGATCCGTTCCTGAACGCACTGCGTCGTGAACGTGTTCCAGTTTCTATATATTTGATGAAAGATATTAAGCTGCAAGGGCAAAACGAGTCTTTTGATCAATTCGTGACCCTGTTGAAAAACACGCGCAGCCAGATGGTATACaa                                                                             >  2:123517/1‑149 (MQ=255)
                               aCAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGGTTa                                                                                >  1:6833/1‑144 (MQ=255)
                                            cTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAc                                                              <  2:105793/149‑1 (MQ=255)
                                              gAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAGTGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGCTTTCTACtg                                                            >  2:169265/1‑149 (MQ=255)
                                                     cttcGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACtgtt                                                          <  1:151868/141‑1 (MQ=255)
                                                        cGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGtc                                                  <  2:231317/149‑1 (MQ=255)
                                                                        cAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCtcat                                  <  2:188834/149‑1 (MQ=255)
                                                                          gtctCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCtcatca                                <  2:250534/146‑1 (MQ=255)
                                                                                 tttatttGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTCTGATCAGTTCGTGATCCTGTTATGCTGGGCTGGCTGCTGGGATTGTTTCTAGCTCGCGATTTCTATTGTTGTCCCGTCTCTCCCGGTTTCTCATCACAGTaac                         <  2:124312/149‑1 (MQ=255)
                                                                                            aaTGGCATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGCGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGcc                   <  1:61769/144‑1 (MQ=255)
                                                                                                        cTGCAAGGGCAAATCGAGTCTTTTGATCAGTCCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTATTGTTGTCTTGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACcagcag  <  2:7541/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                
GAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAG  >  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: