Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,390,286 C→A 81.4% T61N (ACT→AAT)  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,2860CA81.4% 62.8 / 14.3 38T61N (ACT→AAT) hfqglobal sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication
Reads supporting (aligned to +/- strand):  ref base C (5/2);  new base A (15/16);  total (20/18)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.94e-01

GGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGTTTTT  >  CP009273/4390155‑4390430
                                                                                                                                   |                                                                                                                                                 
ggAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCGAATGTTGTCCCGt                                                                                                                                        >  1:364557/1‑143 (MQ=255)
        ttCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTtct                                                                                                                          <  2:248483/149‑1 (MQ=255)
        ttCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTtct                                                                                                                          <  2:353327/149‑1 (MQ=255)
         tCCATTTTCTATGTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTGACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCgg                                                                                                                               <  1:458582/143‑1 (MQ=255)
          ccAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCtca                                                                                                                        >  2:513339/1‑149 (MQ=255)
                   atttatttGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTCCAAGCAAGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCtcatca                                                                                                                     >  1:428392/1‑143 (MQ=255)
                        tgtGGTGAATGGTATTAAGGTGGCAGGGCAAATCGAGGCTTTGGATCAGTTCGTGATCCTGTTGAAAAACACGGTCATCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTa                                                                                                               <  1:463064/141‑1 (MQ=255)
                               aagggTATTAGGCTGGGAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCCGGTGAGAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGc                                                                                                        <  1:202002/140‑1 (MQ=255)
                                          gCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACcagca                                                                                       <  2:269208/149‑1 (MQ=255)
                                            tGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTAcc                                                                                           <  1:303371/143‑1 (MQ=255)
                                              cAAGGGCAAATCGAGTCTTTTGATTAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACcag                                                                                         <  1:211912/143‑1 (MQ=255)
                                              cAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACcac                                                                                         >  1:115201/1‑142 (MQ=255)
                                                  ggCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGt                                                                                     >  1:211262/1‑143 (MQ=255)
                                                          gAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTCCCATCATga                                                                       >  2:557453/1‑148 (MQ=255)
                                                            gTCTTTTGATCATTTTGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTa                                                                     <  2:398295/149‑1 (MQ=255)
                                                                   gATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTag                                                                    <  1:419355/143‑1 (MQ=255)
                                                                     tCAGTTCGTGATCCTGTTGAAAAACACGGTCAGGCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTTTCGCCCTGTTTCTCATCACAGTAACAACGCCGGTTGCGGTACCAGCAGTAACTACCATCATGGTGGCAgcgcgg                                                            >  2:74995/1‑148 (MQ=255)
                                                                           cGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCACATAGGTTACATTGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATa                                                      >  2:464370/1‑149 (MQ=255)
                                                                           cGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAgcgcgc                                                            >  1:406399/1‑143 (MQ=255)
                                                                             tGATCCTGTTGAAACACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCGCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACt                                                    <  2:100381/149‑1 (MQ=255)
                                                                               aTCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGaa                                                        >  1:25873/1‑143 (MQ=255)
                                                                                  cTGTTGAAACATCCGGGTAGTCAGATGGTTTACACTCACGCGATTTCTAATGTTTTCCCGCCTCGCCCGGTTTCTCAGCACAGTAACAACGCCGGTTGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCcgc                                               <  2:294371/149‑1 (MQ=255)
                                                                                    ggtGAAAGACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACtt                                                   <  1:461450/141‑1 (MQ=255)
                                                                                     ttGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTc                                                  >  1:213950/1‑143 (MQ=255)
                                                                                        aaaaaCACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACa                                         <  2:29541/149‑1 (MQ=255)
                                                                                        aaaaaCACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGGACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTCCCGCCCAACa                                         >  2:569782/1‑149 (MQ=255)
                                                                                          aaaCACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCAGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAAGTACCATCAGGGTAGCAGCGCGCAGAATACTTCCgcgc                                             >  1:120651/1‑143 (MQ=255)
                                                                                            acacGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCaa                                           >  1:397210/1‑143 (MQ=255)
                                                                                                     gCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACCGCgaagac                            >  2:296849/1‑148 (MQ=255)
                                                                                                      ccAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGaaa                           >  2:111482/1‑149 (MQ=255)
                                                                                                      ccAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCg                                 >  1:170297/1‑143 (MQ=255)
                                                                                                      ccAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGaaa                           >  2:9290/1‑149 (MQ=255)
                                                                                                      ccAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGaaa                           >  2:97043/1‑149 (MQ=255)
                                                                                                           tGGTTTACAAGCACGCGATTTCTAATTTTGGCACGTCTCGCCCGGTTTCTCATCACATTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAATAGGACAGCGAAGAAACCGaa                      <  2:611480/149‑1 (MQ=255)
                                                                                                                    cgCACGCGATTTCTAATGTTGCCGCGGCTCGGCCGGTTTCTCATCACAGTAACAACGCCGTGAGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATaa                   <  1:188522/142‑1 (MQ=255)
                                                                                                                          cgATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACCGGACCGCCAAGAAACCGAATAAGGCTTc             >  1:531391/1‑143 (MQ=255)
                                                                                                                             tttCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGttt    <  2:168765/149‑1 (MQ=255)
                                                                                                                               tCTAATGTTTTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGttttt  <  2:238659/149‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                                 
GGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGTTTTT  >  CP009273/4390155‑4390430

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: