Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 659,206 A→T 7.3% intergenic (+56/+45) 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,2060AT 7.3% 187.0 / 6.4 55intergenic (+56/+45)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 = 5.72e-01

GTTCAAGAAGGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCAGCGCTTCGGCTT  >  NC_000913/659060‑659354
                                                                                                                                                  |                                                                                                                                                    
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                              <  2:548215/136‑1 (MQ=255)
                                                                                                                                   aaaaaaCCCGCTTTTTGAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATgg                              >  1:548215/1‑136 (MQ=255)
                                                                                                                                   aaaaaaCCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGa                                                                                                          >  3:797025/1‑60 (MQ=255)
                                                                                                                                   aaaaaaCCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGa                                                                                                          <  4:797025/60‑1 (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)
                                                                                                                                                  aaaaGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCAGCGCTTCGGCtt  >  1:766383/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GTTCAAGAAGGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCAGCGCTTCGGCTT  >  NC_000913/659060‑659354

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: