Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,390,235 T→C 43.0% 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,2350TC43.0% 2.9 / 26.2 28I44T (ATC→ACC) hfqglobal sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication
Reads supporting (aligned to +/- strand):  ref base T (8/8);  new base C (6/6);  total (14/14)
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.83e-01

ATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGG  >  CP009273/4390089‑4390361
                                                                                                                                                  |                                                                                                                               
ataAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATcc                                                                                                                               >  2:194560/1‑149 (MQ=255)
            gagaATGGCTTAGGGGCAATCTTTTCAAGTTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAaca                                                                                                                   <  2:116800/149‑1 (MQ=255)
              gaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACg                                                                                                                 <  1:331158/149‑1 (MQ=255)
                       aGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAg                                                                                                        >  2:284026/1‑149 (MQ=255)
                           gCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATgg                                                                                                    >  1:208134/1‑149 (MQ=255)
                            cAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGt                                                                                                   <  1:198340/149‑1 (MQ=255)
                                      aaGAGCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCAc                                                                                         <  1:189559/149‑1 (MQ=255)
                                                 cTGACCGCACTGCGTCGGGAACGTTTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCGAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAc                                                                              <  1:126682/149‑1 (MQ=255)
                                                    aaCGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACtgt                                                                           >  1:118285/1‑149 (MQ=255)
                                                             cGTCGGGGACGTGGTCCAGTTTCTATTTCTTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCCGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTATACAAGCACGCGAGTTCTACAGTTGTCCCGtc                                                                  >  1:201696/1‑149 (MQ=255)
                                                                  ggAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATAGCATAAGTGTATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGcc                                                             <  2:171938/149‑1 (MQ=255)
                                                                       gtggTCGAGTGTGTATCTGTTTGGTGAAGGGCAGTGAGCTGCCGGGGCAAATCGAGTCGTTTGATCAGTGCGTGCCCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGCTTTCTACTGTTGTCCCGTCTCGCCCGGtt                                                        <  1:96471/149‑1 (MQ=255)
                                                                           tCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTtctc                                                    <  1:194560/149‑1 (MQ=255)
                                                                                     atttatttGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTaa                                         >  1:190900/1‑149 (MQ=255)
                                                                                     atttatttGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACCGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTaa                                         >  1:58699/1‑149 (MQ=255)
                                                                                     atgtatttGGTGAATGGTATTAAGCTTCAAGGGCAAATCTAGTCTTTTGATCAGTTCGTGATCCTGTTAAAAAACACGGTCAGCCAGATGGTTTACACGCAAGCGTTATCTACTTTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTaa                                         >  1:342171/4‑149 (MQ=255)
                                                                                        tatttGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTaacaa                                      <  2:50408/149‑1 (MQ=255)
                                                                                                 aaTGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGCTTTATTATGACAGTAACAACGCCGGTgg                             >  2:331789/1‑149 (MQ=255)
                                                                                                      tATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTa                        <  2:233270/149‑1 (MQ=255)
                                                                                                          aaGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACcag                    >  2:284294/1‑149 (MQ=255)
                                                                                                          aaGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTCCAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTCCcaa                    >  1:6382/1‑148 (MQ=255)
                                                                                                          aaGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCCTGCACGCTACAACGCCGGTGGCGGTACcag                    >  1:162678/1‑149 (MQ=255)
                                                                                                            gCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACcagca                  <  1:284026/149‑1 (MQ=255)
                                                                                                                    ggCAAACCGAGTCTTTTGATCAGTTCGTGACCGTGTTGAAAAGCACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTAc          <  1:201646/149‑1 (MQ=255)
                                                                                                                     gCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTAcc         >  2:92301/1‑149 (MQ=255)
                                                                                                                       aaaTCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACcat       <  2:74878/149‑1 (MQ=255)
                                                                                                                           cGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCCGTACCAGCAGTAACTACCATCATg   <  2:181092/149‑1 (MQ=255)
                                                                                                                            gAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATgg  >  2:308360/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                               
ATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGG  >  CP009273/4390089‑4390361

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: