Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,390,178 A→C 50.3% Y25S (TAT→TCT)  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,1780AC50.3% ‑5.4 / 12.0 15Y25S (TAT→TCT) hfqglobal sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication
Reads supporting (aligned to +/- strand):  ref base A (5/2);  new base C (5/2);  total (10/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAG  >  CP009273/4390031‑4390319
                                                                                                                                                   |                                                                                                                                             
tGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCAAGGTTCtatttatt                                                                                                                                             >  1:32723/1‑150 (MQ=255)
   gaCGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGGGTTCCAGTTTCTATCTCTTTgg                                                                                                                                          >  2:275121/1‑150 (MQ=255)
   gaCGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTAAAAATAAGCATATAAGGAAAGAGGGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTCTTTggg                                                                                                                                         >  2:344421/1‑149 (MQ=255)
               gcAATTTGTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACCAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTCTTTGGTGAATGGTATTa                                                                                                                              <  2:167460/150‑1 (MQ=255)
                                gAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATc                                                                                                             >  2:339326/1‑150 (MQ=255)
                                               cAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTCTTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAg                                                                                              >  2:253536/1‑150 (MQ=255)
                                                            aaGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCGGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATGCTGt                                                                                 >  2:139028/1‑150 (MQ=255)
                                                                                        aTCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTCTTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTa                                                     <  1:55848/150‑1 (MQ=255)
                                                                                                       gTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTCTTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCTGATGGTTTACAAGCACGCGATTTc                                      >  2:87421/1‑150 (MQ=255)
                                                                                                           cTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACt                                  <  1:98746/150‑1 (MQ=255)
                                                                                                                  cACTGCGTCGTGAACGTGTTCCAGTTTCTATTTCGTTGGTGAATGGTATTAAGCTGCAAGGGCAAATAGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGCTTTCTACTGTTGTcc                           >  2:140856/1‑150 (MQ=255)
                                                                                                                            ggAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGccc                 >  2:109974/1‑150 (MQ=255)
                                                                                                                            ggAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCAAGCGATTTCTACTGTTGTCCCGTCTCGccc                 >  1:309507/1‑150 (MQ=255)
                                                                                                                                 gtgtTCCAGGTTCTATTTTTTTGGTGAATGGTGTTGGGCTGCAAGGGCAACTCGAGTCTTTTGATCAGTTCGTGATCCTGTGGAAGTGCACGGTCAGCCAGCTGGTTTACAGGCACGCGATTTCTACTGTTCTCCCGTCTCGCCCGGttt            <  1:201415/150‑1 (MQ=255)
                                                                                                                                     tCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCtca        >  2:341394/1‑150 (MQ=255)
                                                                                                                                           ttCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAg  <  1:339326/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
TGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAG  >  CP009273/4390031‑4390319

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: