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

GGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACA‑GTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGTTTTT  >  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)
                        tgtGGTGAATGGTATTAAGGTGGCAGGGCAAATCGAGGCTTTGGATCAGTTCGTGATCCTGTTGAAAAACACGGTCATCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACA‑GTa                                                                                                               <  1:463064/141‑1 (MQ=255)
                               aagggTATTAGGCTGGGAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCCGGTGAGAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACA‑GTAACAACGc                                                                                                        <  1:202002/140‑1 (MQ=255)
                                          gCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACA‑GTAACAACGCCGGTGGCGGTACcagca                                                                                       <  2:269208/149‑1 (MQ=255)
                                            tGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACA‑GTAACAACGCCGGTGGCGGTAcc                                                                                           <  1:303371/143‑1 (MQ=255)
                                              cAAGGGCAAATCGAGTCTTTTGATTAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACA‑GTAACAACGCCGGTGGCGGTACcag                                                                                         <  1:211912/143‑1 (MQ=255)
                                              cAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACA‑GTAACAACGCCGGTGGCGGTACcac                                                                                         >  1:115201/1‑142 (MQ=255)
                                                  ggCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACA‑GTAACAACGCCGGTGGCGGTACCAGCAGt                                                                                     >  1:211262/1‑143 (MQ=255)
                                                          gAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACA‑GTAACAACGCCGGTGGCGGTACCAGCAGTAACTCCCATCATga                                                                       >  2:557453/1‑148 (MQ=255)
                                                            gTCTTTTGATCATTTTGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACA‑GTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTa                                                                     <  2:398295/149‑1 (MQ=255)
                                                                   gATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACA‑GTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTag                                                                    <  1:419355/143‑1 (MQ=255)
                                                                     tCAGTTCGTGATCCTGTTGAAAAACACGGTCAGGCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTTTCGCCCTGTTTCTCATCACA‑GTAACAACGCCGGTTGCGGTACCAGCAGTAACTACCATCATGGTGGCAgcgcgg                                                            >  2:74995/1‑148 (MQ=255)
                                                                           cGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCA‑CATAGGTTACATTGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATa                                                      >  2:464370/1‑149 (MQ=255)
                                                                           cGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACA‑GTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAgcgcgc                                                            >  1:406399/1‑143 (MQ=255)
                                                                             tGATCCTGTTGAAACACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCGCATCACA‑GTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACt                                                    <  2:100381/149‑1 (MQ=255)
                                                                               aTCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACA‑GTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGaa                                                        >  1:25873/1‑143 (MQ=255)
                                                                                  cTGTTGAAACATCCGGGTAGTCAGATGGTTTACACTCACGCGATTTCTAATGTTTTCCCGCCTCGCCCGGTTTCTCAGCACA‑GTAACAACGCCGGTTGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCcgc                                               <  2:294371/149‑1 (MQ=255)
                                                                                    ggtGAAAGACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACA‑GTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACtt                                                   <  1:461450/141‑1 (MQ=255)
                                                                                     ttGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACA‑GTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTc                                                  >  1:213950/1‑143 (MQ=255)
                                                                                        aaaaaCACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACA‑GTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACa                                         <  2:29541/149‑1 (MQ=255)
                                                                                        aaaaaCACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACA‑GTAACAACGCCGGTGGCGGGACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTCCCGCCCAACa                                         >  2:569782/1‑149 (MQ=255)
                                                                                          aaaCACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCAGGTTTCTCATCACA‑GTAACAACGCCGGTGGCGGTACCAGCAGTAAGTACCATCAGGGTAGCAGCGCGCAGAATACTTCCgcgc                                             >  1:120651/1‑143 (MQ=255)
                                                                                            acacGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACA‑GTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCaa                                           >  1:397210/1‑143 (MQ=255)
                                                                                                     gCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACA‑GTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACCGCgaagac                            >  2:296849/1‑148 (MQ=255)
                                                                                                      ccAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACA‑GTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGaaa                           >  2:111482/1‑149 (MQ=255)
                                                                                                      ccAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACA‑GTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCg                                 >  1:170297/1‑143 (MQ=255)
                                                                                                      ccAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACA‑GTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGaaa                           >  2:9290/1‑149 (MQ=255)
                                                                                                      ccAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACA‑GTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGaaa                           >  2:97043/1‑149 (MQ=255)
                                                                                                           tGGTTTACAAGCACGCGATTTCTAATTTTGGCACGTCTCGCCCGGTTTCTCATCACA‑TTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAATAGGACAGCGAAGAAACCGaa                      <  2:611480/149‑1 (MQ=255)
                                                                                                                    cgCACGCGATTTCTAATGTTGCCGCGGCTCGGCCGGTTTCTCATCACA‑GTAACAACGCCGTGAGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATaa                   <  1:188522/142‑1 (MQ=255)
                                                                                                                          cgATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACA‑GTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACCGGACCGCCAAGAAACCGAATAAGGCTTc             >  1:531391/1‑143 (MQ=255)
                                                                                                                             tttCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACA‑GTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGttt    <  2:168765/149‑1 (MQ=255)
                                                                                                                               tCTAATGTTTTCCCGTCTCGCCCGGTTTCTCATCACA‑GTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGttttt  <  2:238659/149‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                                 
GGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACA‑GTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGTTTTT  >  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: ‑