Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 551,137 G→T T233T (ACG→ACT ybcF → putative carbonate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913551,1370GT98.1% 163.4 / NA 53T233T (ACG→ACTybcFputative carbonate kinase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (26/26);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

ACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACGCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGAT  >  NC_000913/551006‑551268
                                                                                                                                   |                                                                                                                                   
acggatgacggaGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCgcagc                                                                                                                              >  1:1081995/1‑139 (MQ=255)
    atgacggaGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACg                                                                                                                          >  1:748678/1‑139 (MQ=255)
    atgacggaGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACg                                                                                                                          >  2:1090813/1‑139 (MQ=255)
     tgacggaGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGt                                                                                                                         >  1:1009515/1‑139 (MQ=255)
         ggaGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCa                                                                                                                     <  1:1090813/139‑1 (MQ=255)
             cAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCACCGTGccttacg                                                                                                                 <  1:708102/139‑6 (MQ=255)
             cAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCg                                                                                                                 <  1:845989/139‑1 (MQ=255)
                     gAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCca                                                                                                         <  1:1073359/139‑1 (MQ=255)
                           gTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCgg                                                                                                   <  1:976141/139‑1 (MQ=255)
                                 gatAAAGATGTCGCCGCTTCGGTGCTCCCCGAGTAGATGAATGTTGATCGACTGTGGAGCCTCACCGGTGCTGATGCGGTATATGCAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGt                                                                                             >  1:488042/1‑139 (MQ=255)
                                      aGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAgcg                                                                                        >  1:524173/1‑139 (MQ=255)
                                        aTCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGcc                                                                                      <  2:748678/139‑1 (MQ=255)
                                          ctcGCCGCTGTGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATg                                                                                                <  2:760460/127‑1 (MQ=255)
                                          ctcGCCGCTGTGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATg                                                                                                >  1:760460/1‑127 (MQ=255)
                                             gccgcTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTg                                                                                 >  1:983239/1‑139 (MQ=255)
                                             gccgcTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTg                                                                                 >  1:246651/1‑139 (MQ=255)
                                             gccgcTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTg                                                                                 >  1:618388/1‑139 (MQ=255)
                                                  tGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCaaa                                                                            >  2:898625/1‑139 (MQ=255)
                                                  tGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCaaa                                                                            >  1:19553/1‑139 (MQ=255)
                                                     gTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGcc                                                                         <  1:898625/139‑1 (MQ=255)
                                                            gCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGtt                                                                  <  2:19553/139‑1 (MQ=255)
                                                                aGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGAt                                                              <  2:1009515/139‑1 (MQ=255)
                                                                               aTGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCg                                                                                                                 >  2:560106/1‑73 (MQ=255)
                                                                               aTGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCg                                                                                                                 <  1:560106/73‑1 (MQ=255)
                                                                               aTGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCg                                                                                                                 <  2:453401/73‑1 (MQ=255)
                                                                               aTGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCg                                                                                                                 >  1:453401/1‑73 (MQ=255)
                                                                               aTGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAAc                                               >  2:335707/1‑139 (MQ=255)
                                                                                   aCTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAAcggcg                                           <  2:618388/139‑1 (MQ=255)
                                                                                    cTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCaa                                                                             >  1:592637/1‑104 (MQ=255)
                                                                                    cTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCaa                                                                             <  2:592637/104‑1 (MQ=255)
                                                                                       gTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGa                                       <  1:335707/139‑1 (MQ=255)
                                                                                           tCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTgg                                   >  1:374110/1‑139 (MQ=255)
                                                                                            ccTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAAcg                                              <  2:807162/127‑1 (MQ=255)
                                                                                            ccTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAAcg                                              >  1:807162/1‑127 (MQ=255)
                                                                                             cTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATg                                                                                                <  1:991756/76‑1 (MQ=255)
                                                                                             cTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATg                                                                                                >  2:991756/1‑76 (MQ=255)
                                                                                             cTCACCGATGCTGATGCGGTACATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCCGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGtt                                 <  1:710702/139‑1 (MQ=255)
                                                                                                  cGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGa                                       >  1:38606/1‑128 (MQ=255)
                                                                                                  cGATGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGa                                       <  2:38606/128‑1 (MQ=255)
                                                                                                    aTGCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAg                          >  1:24004/1‑139 (MQ=255)
                                                                                                      gCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGa                                                               >  1:780820/1‑100 (MQ=255)
                                                                                                      gCTGATGCGGTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGa                                                               <  2:780820/100‑1 (MQ=255)
                                                                                                        tGATGCGGTATATGAAAACTGGGGAACGCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTGATGTCagcagc                      >  2:788559/1‑139 (MQ=255)
                                                                                                               gTATATGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTa               >  1:999234/1‑139 (MQ=255)
                                                                                                                 atatGAAAACTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTaaa             <  2:24004/139‑1 (MQ=255)
                                                                                                                      aaaaCTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCcgc        <  2:999234/139‑1 (MQ=255)
                                                                                                                       aaaCTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCcgcg       <  2:246651/139‑1 (MQ=255)
                                                                                                                          cTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCAt                                                                                    >  1:342909/1‑59 (MQ=255)
                                                                                                                          cTGGGGAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCAt                                                                                    <  2:342909/59‑1 (MQ=255)
                                                                                                                            ggggAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGt      <  1:853752/135‑1 (MQ=255)
                                                                                                                            ggggAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGt      >  2:853752/1‑135 (MQ=255)
                                                                                                                            ggggAACTCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGAt  >  2:937782/1‑139 (MQ=255)
                                                                                                                                aaCTCCGCAGCAACGTGCCATTCGCCATGCCacac                                                                                                      <  1:434323/35‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                   
ACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACGCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGAT  >  NC_000913/551006‑551268

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