Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 659,205 A→T 7.3% intergenic (+55/+46) tatE → / ← lipA TatABCE protein translocation system subunit/lipoate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913659,2050AT 7.3% 186.2 / 6.4 55intergenic (+55/+46)tatE/lipATatABCE protein translocation system subunit/lipoate synthase
Reads supporting (aligned to +/- strand):  ref base A (25/26);  new base T (2/2);  total (27/28)
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.97e-01

AGGGTTCAAGAAGGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCAGCGCT  >  NC_000913/659057‑659347
                                                                                                                                                    |                                                                                                                                              
aGGGTTCAAGAAGGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCa                                                                                                                                                >  3:415044/1‑149 (MQ=255)
   gTTCAAGAAGGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCaaaa                                                                                                                                             >  3:243607/1‑149 (MQ=255)
      cAAGAAGGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGc                                                                                                                                          <  2:288056/149‑1 (MQ=255)
       aagaagGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCaaaaa                                                                                                                                            >  3:296578/1‑146 (MQ=255)
       aagaagGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCaaaaa                                                                                                                                            <  4:296578/146‑1 (MQ=255)
        agaagGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCgg                                                                                                                                        >  3:504923/1‑149 (MQ=255)
            ggCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGttt                                                                                                                                    <  1:231951/149‑1 (MQ=255)
                 tgaatgaTGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTAt                                                                                                                               >  2:632608/1‑149 (MQ=255)
                       atgaCGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATcagaca                                                                                                                         >  1:210596/1‑149 (MQ=255)
                         gaCGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATcagacaga                                                                                                                       >  3:569057/1‑149 (MQ=255)
                               gctgcGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGCTTTTTATCAGACAGATgtaag                                                                                                                 <  4:243607/149‑1 (MQ=255)
                                 tgcGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATgtaagta                                                                                                               >  4:769940/1‑149 (MQ=255)
                                  gcggcgAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATgtaagtaa                                                                                                              >  3:349128/1‑149 (MQ=255)
                                        aaaaaaGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAAttatta                                                                                                        >  3:2809/1‑149 (MQ=255)
                                          aaaaGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACa                                                                                                      <  1:632608/149‑1 (MQ=255)
                                           aaaGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAg                                                                                                     <  4:415044/149‑1 (MQ=255)
                                                  cAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACt                                                                                              <  3:79390/149‑1 (MQ=255)
                                                  cAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACt                                                                                              <  1:254014/149‑1 (MQ=255)
                                                        ttGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACtt                                                                                        <  4:504923/149‑1 (MQ=255)
                                                         tGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTc                                                                                       <  3:392289/149‑1 (MQ=255)
                                                                  ggCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCttt                                                                              >  2:320709/1‑149 (MQ=255)
                                                                    cTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCg                                                                            >  4:765481/1‑149 (MQ=255)
                                                                     taaaaaGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATgtaagtaa                                                                                                              >  1:644042/3‑114 (MQ=255)
                                                                     tGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATgtaagtaa                                                                                                              <  2:644042/114‑1 (MQ=255)
                                                                        aaaGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAGAAGCGGGtt                                                                                                                                     >  1:13131/1‑88 (MQ=255)
                                                                        aaaGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGtt                                                                                                                                     <  2:13131/88‑1 (MQ=255)
                                                                         aaGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGc                                                                                                                                          <  4:338853/82‑1 (MQ=255)
                                                                         aaGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGc                                                                                                                                          >  3:338853/1‑82 (MQ=255)
                                                                           gCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCaa                                                                     >  4:510943/1‑149 (MQ=255)
                                                                             tctctcATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAAt                                                                   <  3:769940/149‑1 (MQ=255)
                                                                             tctctcATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAAt                                                                   <  1:320709/149‑1 (MQ=255)
                                                                                   aTAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGt                                                             >  4:558127/1‑149 (MQ=255)
                                                                                   aTAAAGAGTGACGTGGCGAGCAGGACACTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTTTTGAAGCGGGtttttt                                                                                                                                 <  2:627805/81‑1 (MQ=25)
                                                                                   aTAAAGAGTGACGTGGCGAGCAGGACACTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTTTTGAAGCGGGtttttt                                                                                                                                 >  1:627805/1‑81 (MQ=25)
                                                                                       agagTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTa                                                         >  4:323122/1‑149 (MQ=255)
                                                                                        gagTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTaa                                                        <  1:36933/149‑1 (MQ=255)
                                                                                           tGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTaagaa                                                     <  1:634950/149‑1 (MQ=255)
                                                                                           tGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTaagaa                                                     <  3:510943/149‑1 (MQ=255)
                                                                                                tGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCg                                                <  2:210596/149‑1 (MQ=255)
                                                                                                  gCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGa                                              <  4:349128/149‑1 (MQ=255)
                                                                                                      gCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACaaa                                          >  1:477632/1‑149 (MQ=255)
                                                                                                       cAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAAc                                         <  4:2809/149‑1 (MQ=255)
                                                                                                          gACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGa                                      <  4:209993/149‑1 (MQ=255)
                                                                                                                     aaTATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATgcagc                           >  1:345107/1‑149 (MQ=255)
                                                                                                                           ttGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGt                     <  3:323122/149‑1 (MQ=255)
                                                                                                                                   tACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCa             <  3:558127/149‑1 (MQ=255)
                                                                                                                                      aaaaaaCCCGCTTTTTGAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATgg                       >  1:548215/1‑136 (MQ=255)
                                                                                                                                      aaaaaaCCCGCTTTTTGAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATgg                       <  2:548215/136‑1 (MQ=255)
                                                                                                                                      aaaaaaCCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGa                                                                                                   <  4:797025/60‑1 (MQ=255)
                                                                                                                                      aaaaaaCCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGa                                                                                                   >  3:797025/1‑60 (MQ=255)
                                                                                                                                      aaaaaaCCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTGGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCg          >  1:799976/1‑149 (MQ=255)
                                                                                                                                      aaaaaaCCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACgg                                       <  1:778680/120‑1 (MQ=255)
                                                                                                                                      aaaaaaCCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACgg                                       >  2:778680/1‑120 (MQ=255)
                                                                                                                                      aaaaaaCCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCg          >  3:317110/1‑149 (MQ=255)
                                                                                                                                      aaaaaaCCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCg          >  2:785172/1‑149 (MQ=255)
                                                                                                                                             ccGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCAgcgc   >  1:101494/1‑149 (MQ=255)
                                                                                                                                              cGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCAGCGCt  >  4:605897/1‑149 (MQ=255)
                                                                                                                                               gCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTaa                                                                                           >  2:184432/1‑59 (MQ=255)
                                                                                                                                               gCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTaa                                                                                           <  1:184432/59‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
AGGGTTCAAGAAGGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCAGCGCT  >  NC_000913/659057‑659347

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