Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 551,066 G→T A210S (GCC→TCC)  ybcF → putative carbamate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913551,0660GT97.1% 118.2 / NA 34A210S (GCC→TCC) ybcFputative carbamate kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (16/17);  minor base A (1/0);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGGCGGTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACGCCG  >  NC_000913/550988‑551140
                                                                              |                                                                          
ggcgggggTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCAcc                                                                          >  2:116114/1‑81 (MQ=255)
ggcggcggTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTcc                                                                          >  1:116107/1‑81 (MQ=255)
 gcggcggTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCg                                                                         <  2:116067/81‑1 (MQ=255)
 gcggcggTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCg                                                                         <  1:116066/81‑1 (MQ=255)
      ggTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAg                                                                       >  1:116112/1‑78 (MQ=255)
      ggTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAg                                                                       >  1:116108/1‑78 (MQ=255)
        tgtgCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGAt                                                                  >  2:116113/1‑81 (MQ=255)
                  acggatgacggaGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGAt                                                        <  1:116068/81‑1 (MQ=255)
                  acggatgacggaGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGAt                                                        <  1:116069/81‑1 (MQ=255)
                    ggatgacggaGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATgg                                                      <  2:116070/81‑1 (MQ=255)
                     gatgacggaGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGa                                                     >  2:116100/1‑81 (MQ=255)
                     gatgacggaGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGa                                                     >  1:116122/1‑81 (MQ=255)
                      atgacggaGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGAc                                                    >  1:116116/1‑81 (MQ=255)
                        gacggaGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTg                                                  >  1:116109/1‑81 (MQ=255)
                          cggaGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGt                                                <  2:116071/81‑1 (MQ=255)
                           ggaGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGt                                                >  1:116118/1‑80 (MQ=255)
                            gaGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGa                                              >  2:116105/1‑81 (MQ=255)
                                 gggAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTc                                         <  1:116072/81‑1 (MQ=255)
                                        aaGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGAt                                   <  1:116073/80‑1 (MQ=255)
                                        aaGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGAt                                   <  1:116074/80‑1 (MQ=255)
                                          gCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCt                                >  1:116111/1‑81 (MQ=255)
                                          gCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCt                                <  1:116075/81‑1 (MQ=255)
                                           cAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCt                                >  1:116130/1‑80 (MQ=255)
                                           cAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTg                               <  1:116076/81‑1 (MQ=255)
                                           cAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTg                               <  1:116077/81‑1 (MQ=255)
                                            aGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGa                              <  1:116078/81‑1 (MQ=255)
                                                attgatAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCg                          >  2:116119/1‑81 (MQ=255)
                                                     tAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGtata                     <  2:116079/81‑1 (MQ=255)
                                                               gccgcTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTgg           >  1:116131/1‑81 (MQ=255)
                                                               gccgcTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTg            <  1:116080/80‑1 (MQ=255)
                                                                 cgcTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTgggg         <  2:116081/81‑1 (MQ=255)
                                                                      cGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAAcgc    <  2:116082/81‑1 (MQ=255)
                                                                          gCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAAcgccg  >  1:116115/1‑79 (MQ=255)
                                                                          gCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAAcgccg  >  2:116120/1‑79 (MQ=255)
                                                                              |                                                                          
GGCGGCGGTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACGCCG  >  NC_000913/550988‑551140

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: