Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 815,815 G→A A278V (GCC→GTC)  ybhK ← putative transferase YbhK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913815,8150GA80.8% 54.8 / ‑4.9 26A278V (GCC→GTC) ybhKputative transferase YbhK
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (13/8);  minor base C (0/5);  total (13/13)
Fisher's exact test for biased strand distribution p-value = 3.91e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCTGGGCGATTTCCTTCAGCAGCAGAATTGGCATCAGGC  >  NC_000913/815577‑816061
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      
ggAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGACCTCCAGTAcct                                                                                                                                                                                                                                            <  2:357396/251‑1 (MQ=255)
                    aCGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGACCTCCAGTACCTCCTGGATCACAATCCGctct                                                                                                                                                                                                                        >  2:80241/1‑251 (MQ=255)
                                           cgcggcaaaattttaattttgtgggTTTGGTAGGCCGGAAAAGGGGTTAACGCCGCATCCGGAATAAACAACGCGCTCTTTGTCTTAAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGAACCGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTAATGACGATACGGAATATCGTTGACCTCCAGTCCCTCCGGGAGCACAATCCGCTCTTTCACCGCCGAGACATCCACtt                                                                                                                                                                                                  <  2:253495/228‑1 (MQ=255)
                                                        gAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGACCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATg                                                                                                                                                                                    >  1:231827/1‑251 (MQ=255)
                                                                 cccattcttAGGCCGGATAGGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACGGGCGGTCATGACGATACGGAATATCGCTGACCTCCAGTCCCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATc                                                                                                                                                                            <  1:137707/243‑1 (MQ=255)
                                                                         tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGACCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATc                                                                                                                                                                            >  2:137707/1‑243 (MQ=255)
                                                                         gagccggaaaaagggtttCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCGTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCAAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACACCTGGCGGTCATGACGATACGGAATATCGCTGCCCTCCAGTCCCTCCTGGATCACAATCCCCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCt                                                                                                                                                                    <  2:123596/235‑1 (MQ=255)
                                                                                    aGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGACCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGCGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAAc                                                                                                                                                          >  2:171959/1‑250 (MQ=255)
                                                                                              gccgcATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGACCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTccc                                                                                                                                               >  1:291488/1‑250 (MQ=255)
                                                                                                                    cgcACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGACCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTcc                                                                                                                           >  1:205569/1‑249 (MQ=255)
                                                                                                                     gcACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGACCTCCAGTACCTCCTGGATCACAATcc                                                                                                                                                                                                                             <  2:275145/150‑1 (MQ=255)
                                                                                                                     gcACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGACCTCCAGTACCTCCTGGATCACAATcc                                                                                                                                                                                                                             >  1:275145/1‑150 (MQ=255)
                                                                                                                      cACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGGAACCCCTTTTCCAGCGCGTTATGTAACAACGGGCGGTCATGACGATACGGAATATCGCTGACCTCCAGTACCTCCGGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCt                                                                                                                       <  1:171959/251‑1 (MQ=255)
                                                                                                                        ttttGTCTTCAATTGAAACCCGGATTTTTCCTTCGGGTTTACGTTAACCTAAAGCCTTTACCGCCTTTTCCCGCCCTTTATGTAACAACTGGCGGTCATGCCGATACGAAATATCGCTGCCCTCAAGTCCCTCCTGGATCACATCCCCCTCTTTCACCGCCGATACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCt                                                                                                                       <  2:227919/248‑1 (MQ=255)
                                                                                                                            gTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTACCCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGACCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAAt                                                                                                                 <  1:21743/251‑1 (MQ=255)
                                                                                                                              cTTCAATAGGAAGACGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGACCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCCCCATAATTGCCAGCTTGCTTTCCAGCTTCAAAtt                                                                                                                >  2:253044/1‑250 (MQ=255)
                                                                                                                                                                                      gCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGACCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTgg                                                       >  1:10462/1‑251 (MQ=255)
                                                                                                                                                                                               aGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGACCTCCAGTAACTCATGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCACATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTgcggcg                                              >  1:170077/1‑251 (MQ=255)
                                                                                                                                                                                                       aTGTAACACCGGGCGGCCATGACGATCCGGAATATCGCTGCCCTCCAGTACCTCCTGAATCACAATCCCCTCTTTCCCCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCt                                      <  2:132395/251‑1 (MQ=255)
                                                                                                                                                                                                       aTGTAACAACTGGCGGTCATGACGATACGAAATATCGCTGCCCTCCAGTCCCTCCGGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCt                                      <  1:80741/251‑1 (MQ=255)
                                                                                                                                                                                                          aaacaacTGGGGGTCATGACGATACGGAATATCGCTGACCTCCAGTACCTCCTGGATCACAACCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCTggg                                   <  1:232971/250‑1 (MQ=255)
                                                                                                                                                                                                                 tGGCGGTCATGACGATACGGAATATCGCTGACCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCTGGGCGATTTc                            >  1:235816/1‑251 (MQ=255)
                                                                                                                                                                                                                         aTGACGATACGGAATATCGCTGACCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACAGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCTGGGCGATTTCCTTcagc                     >  1:180359/1‑250 (MQ=255)
                                                                                                                                                                                                                                     aaTATCGCTGACCTCCAGTCCCTCCTGGATCCCAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCTGGGCGATTTCCTTCAGCAGCAGAATTGGCa        <  1:91718/251‑1 (MQ=255)
                                                                                                                                                                                                                                          cGCTGACCTCAAGTACCTCCTGGATCACATTCTGCTCTTTCACCGCAGAGCCATCCACTTTTGGCCCGACGATGCCCGCATCAATGACATTTTTACCAAAATACTGATCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGACGCAGGTACACTCAACTCAAGCCCTAGATTGCCGATATAAACCATCGGTGCTGGCGTGCGGCGTAATGCCTGGGCGGTTACATTCAGCAACAGTATAGGCAtacgg   >  2:351067/1‑247 (MQ=255)
                                                                                                                                                                                                                                           gCTGCCCTCCAGTCCCTCCTGTATCACAATCCGCTCTTCCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCTGGGCTATTTCCTTCAGCAGCAGAATTTGCATCAGGc  <  1:102354/251‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      
GGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCTGGGCGATTTCCTTCAGCAGCAGAATTGGCATCAGGC  >  NC_000913/815577‑816061

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: