Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 4,390,205 G→T G34V (GGG→GTG)  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,2050GT97.9% 146.3 / NA 47G34V (GGG→GTG) hfqglobal sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (27/19);  total (27/20)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAG  >  CP009273/4390060‑4390346
                                                                                                                                                 |                                                                                                                                             
aTCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCa                                                                                                                                            >  1:360744/1‑149 (MQ=255)
aTCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTGCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCa                                                                                                                                            >  1:533728/1‑149 (MQ=255)
 tCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGGCTTCCAGTTTCTATGTATTTGGTGAATGGTATTAAGCTGCAAGTGCaa                                                                                                                                           >  1:567129/1‑149 (MQ=255)
      aGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGa                                                                                                                                      <  1:151593/149‑1 (MQ=255)
      aGGTTCAAAGAACAAATAAGCATATAAGGAAAAGATAGAATGGCCAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGa                                                                                                                                      <  1:387427/149‑1 (MQ=255)
          tCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCAAGTTTCTATTTATTTGGTGACTGGTATTAAGCTGCAAGTGCAAATc                                                                                                                                        >  2:353233/1‑143 (MQ=255)
                tACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCACGTGCAAATCGAGTCt                                                                                                                                  >  2:432096/1‑143 (MQ=255)
                tACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCt                                                                                                                                  >  2:623532/1‑143 (MQ=255)
                tACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCt                                                                                                                                  >  2:782411/1‑143 (MQ=255)
                     aTAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGtt                                                                                                                       >  1:699983/1‑149 (MQ=255)
                      tAAGCATATAAGGCAAAGAGAGAATGGCTAAGGGGCAACCTTTACAAGGTCCGTTCCTGAACGCACTGCGTCGGGAAGGTGTTCCAGGTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGAt                                                                                                                            <  2:627076/143‑1 (MQ=255)
                          cATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAACCGAGTCTTTTGATCAGTTCGTGa                                                                                                                  <  1:600979/149‑1 (MQ=255)
                           atatAAGGCAAAGAGAGCATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCGGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGtt                                                                                                                       <  2:687212/143‑1 (MQ=255)
                                aGGGAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGa                                                                                                                  <  2:408239/143‑1 (MQ=255)
                                aGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGt                                                                                                            <  1:691798/149‑1 (MQ=255)
                                      agagagAAGGGCTAAGGGGCAACCTTTACAAGATCCGTTCCTGAACGCACTGCGTGCGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGCGTCTTTTGATCAGTTCGTGATCCTGTTGaaaa                                                                                                      <  1:101851/149‑1 (MQ=255)
                                          agaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAacac                                                                                                  <  1:194520/149‑1 (MQ=255)
                                           gaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGTTCAGTTCGTGATCCTGTTGAAAAACACg                                                                                                 <  1:491850/149‑1 (MQ=255)
                                            aaTGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACgg                                                                                                >  1:194105/1‑149 (MQ=255)
                                                  tAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGcc                                                                                          <  1:522994/149‑1 (MQ=255)
                                                       ggCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAACCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGc                                                                                           >  2:750067/1‑143 (MQ=255)
                                                              tttACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACa                                                                              <  1:737476/149‑1 (MQ=255)
                                                                          gTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGCAAAACACGGTCAGCCCGATGGTTTCCAAGCCcg                                                                        >  2:816671/1‑143 (MQ=255)
                                                                               tGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACt                                                             <  1:853020/149‑1 (MQ=255)
                                                                                gAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTGTTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGACTGATGTTTTACAAGCACGCGAttt                                                                  >  2:584229/1‑143 (MQ=255)
                                                                                 aaggCACTGCGTCGGGACCGTGTTCCAGTTGCTATGTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTAAAACAAAAAAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTc                                                                 <  2:224057/140‑1 (MQ=255)
                                                                                      aCTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTcc                                                      <  1:94895/149‑1 (MQ=255)
                                                                                                    gtgtTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGtt                                        >  1:438731/1‑149 (MQ=255)
                                                                                                       ttCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGACAGATGGTTTACAAGCCCGCGATTTCTACTGTTGTCCCGTCTCGCCCg                                           >  2:153902/1‑143 (MQ=255)
                                                                                                          cAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCtcat                                  <  1:695992/149‑1 (MQ=255)
                                                                                                             tttCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGAGTTCTACTGTTGTCCCGTCTCGCCCGGTTtct                                     >  2:310074/1‑143 (MQ=255)
                                                                                                                  atttatttGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTaa                          >  1:835260/1‑149 (MQ=255)
                                                                                                                   tttatttGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATcac                               >  2:260640/1‑143 (MQ=255)
                                                                                                                   tttatttGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTaac                         <  1:306178/149‑1 (MQ=255)
                                                                                                                        ttGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGc                    >  1:160916/1‑149 (MQ=255)
                                                                                                                            tGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTaacaac                      >  2:683384/1‑143 (MQ=255)
                                                                                                                            tGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGt                >  1:620193/1‑149 (MQ=255)
                                                                                                                            tGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGAATTCTACTGTTGTCCCGTCTCACCCGGGATCTCATCACAGTAACAAGGCCGGt                >  1:891235/1‑149 (MQ=255)
                                                                                                                               aTGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGcc                   >  2:100432/1‑143 (MQ=255)
                                                                                                                                tGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCg                  >  2:93060/1‑143 (MQ=255)
                                                                                                                                   tATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTCACAACGCCGGTGGCGGTa         >  1:536823/1‑149 (MQ=255)
                                                                                                                                    aTTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTAc        <  1:118951/149‑1 (MQ=255)
                                                                                                                                    aTTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGACGGCTTCCAAGCACGCGACTCCTGCGATTGTCCCTACGCGCACGGTGGCTCCTCGCCGTAACCCCGCCggccagggaca        >  1:191860/1‑140 (MQ=255)
                                                                                                                                     ttAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGc             >  2:827044/1‑143 (MQ=255)
                                                                                                                                      tAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCg            <  2:637635/143‑1 (MQ=255)
                                                                                                                                       aaGCTGCAAGTGCAAATCGAGTCTTTTGATCCACTCGATGCCCTGGTGAAAAACACGGTCAGACAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACcag     <  1:622491/149‑1 (MQ=255)
                                                                                                                                          cTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACcagcag  >  1:173889/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
ATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAG  >  CP009273/4390060‑4390346

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: ‑