Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 400,404 A→G 9.7% S13P (TCA→CCA)  tatC ← Tat family twin arginine targeting transporter TatC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730400,4040AG 9.7% 104.6 / 3.0 41S13P (TCA→CCA) tatCTat family twin arginine targeting transporter TatC
Reads supporting (aligned to +/- strand):  ref base A (18/19);  new base G (2/2);  total (20/21)
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.79e-01

ATGAAACAAAACGCGATGAAAAATATAATCATCACGTATATTTGTATCATTTCTGTGAATGAAAATGCATGTAAGGACACATGTGCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGAAACATATACGACGATGACCGTAATGACGAACGACAATAAAATTTTTACTAACCGATGGCGTAGTTCGCTAAAGTGAACCAGTAAGGTGTGGTCAGTGCTATTGCTCTCGCTGTTGTTTCGATTCCTTACTGGGTGTA  >  CP000730/400274‑400541
                                                                                                                                  |                                                                                                                                         
atgaaaCAAAACGCGATGAAAAATATAATCATCACGTATATTTGTATCATTTCTGTGAATGAAAATGCATGTAAGGACACATGTGCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGAAACATATacg                                                                                                                                 <  2:155436/141‑1 (MQ=255)
atgaaaCAAAACGCGATGAAAAATATAATCATCACGTATATTTGTATCATTTCTGTGAATGAAAATGCATGTAAGGACACATGTGCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGAAACATATacg                                                                                                                                 <  1:16844/141‑1 (MQ=255)
    aaCAAAACGCGATGAAAAATATAATCATCCCGTATATTTGTATCATTTCTGTGAATGAAAATGCATGTAAGGACACATGTGCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGAAACATATAcgacga                                                                                                                             <  1:384396/141‑1 (MQ=255)
    aaCAAAACGCGATGAAAAAGATAATCATCACGTATATTTGTATCATTTCTGTGAATGAAAATGCATGTAAGGACACATGTGCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGAAACATATAcgacga                                                                                                                             <  2:456035/141‑1 (MQ=255)
       aaaaCGCGATGAAAAATATAATCATCACGTATATTTGTATCATTTCTGTGAATGAAAATGCATGTAAGGACACATGTGCAAGGATAATATACATTATGAATGGTGTCATCAACCAAAATGATGAAACATATACGACGATGa                                                                                                                          <  2:304923/141‑1 (MQ=255)
         aaCGCGATGAAAAATATAATCATCACGTATATTTGTATCATTTCTGTGAATGAAAATGCATGTAAGGACACATGTGCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGAAACATATACGACGATGAcc                                                                                                                        >  2:178823/1‑141 (MQ=255)
         aaCGCGATGAAAAATATAATCATCACGTATATTTGTATCATTTCTGTGAATGAAAATGCATGTAAGGACACATGTGCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGAGGAAACATATACGACGATTAcc                                                                                                                        >  1:315504/1‑141 (MQ=255)
              gATGAAAAATATAATCATCACGTATATTTGTATCATTTCTGTGAATGAAAATGCATGTAAGGACACATGTGCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGGAACATATACGACGATGACCGTAAt                                                                                                                   >  1:336494/1‑141 (MQ=255)
              gATGAAAAATATAATCATCACGTATATTTGTATCATTTCTGTGAATGAAAATGCATGTAAGGACACATGTGCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGGAACATATACGACGATGACCGTAAt                                                                                                                   >  1:100166/1‑141 (MQ=255)
               aTGAAAAATATAATCATCACGTATATTTGTATCATTTCTGTGAATGAAAATGCATGTAAGGACACATGTGCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGAAACATATACGACGATGACCGTAATg                                                                                                                  >  1:424068/1‑141 (MQ=255)
                          aaTCATCACGTATATTTGTATCATTTCTGTGAATGAAAATGCATGTAAGGACACATGTGCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGAAACATATACGACGATGACCGTAATGACGAACGACaa                                                                                                       <  2:387138/141‑1 (MQ=255)
                          aaTCATCACGTATATTTGTATCATTTCTGTGAATGAAAATGCATGTAAGGACACATGTGCACGGGTAATATACGTTATGAACGGTGTCATCCACCAAAATGATGAAACATATACGACGATGACCGTAATGACGAACGACaa                                                                                                       >  1:182092/1‑141 (MQ=255)
                            tcatcaCGTATATTTGTATCATTTCTGTGAATGAAAATGCATGTAAGGACACATGTGCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGAAACATATACGACGATGACCGTAATGACGAACGACAATa                                                                                                     <  1:352952/141‑1 (MQ=255)
                                  cGTATATTTGTATCATTTCTGTGAATGAAAATGCATGTAAGGACACATGTGCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGAAACATATACGACGATGACCGTAATGACGAACGACAATAAAAttt                                                                                               >  1:338607/1‑141 (MQ=255)
                                        tttGTATCATTTCTGTGAATGAAAATGCATGTAAGGACACATGTGCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGAAACATATACGACGATGACCGTAATGACGAACGACAATAAAATTTTTACTa                                                                                         >  2:135667/1‑141 (MQ=255)
                                        tttGTATCATTTCTGTGAATGAAAATGCATGTAAGGACACATGTGCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGAAACATATACGACGATGACCGTAATGACGAACGACAATAAAATTTTTACTa                                                                                         >  1:311554/1‑141 (MQ=255)
                                           gTATCATTTCTGTGAATGAAAATGCATGTAAGGACACATGTGCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGAAACATATACGACGATGACCGTAATGACGAACGACAATAAAATTTTTACTAAcc                                                                                      <  2:316089/141‑1 (MQ=255)
                                                aTTTCTGTGAATGAAAATGCATGTAAGGACACATGTGCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGAAACATATACGACGATGACCGTAATGACGAACGACAATAAAATTTTTACTAACCGATgg                                                                                 <  1:14589/141‑1 (MQ=255)
                                                     tgtgAATGAAAATGCATGTAAGGACACATGTGCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGAAACATATACGACGATGACCGTAATGACGAACGACAATAAAATTTTTACTAACCGATGGCGTAg                                                                            >  1:145508/1‑141 (MQ=255)
                                                     tgtgAATGAAAATGCATGTAAGGACACATGTGCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGAAACATATACGACGATGACCGTAATGACGAACGACAATAAAATTTTTACTAACCGATGGCGTAg                                                                            <  1:252364/141‑1 (MQ=255)
                                                       tgaatgaaAATGCATGTAAGGACACATGTGCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATTAAACATGTACGACGATGACCGTAATGACGAACGACAATAAGATTTCTACTAACCGATGGCGTAGtt                                                                          <  2:318914/141‑1 (MQ=255)
                                                          atgaaAATGCATGTAAGGACACATGTGCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGAAACATATACGACGATGACCGTAATGACGAACGACAATAAAATTTTTACTAACCGATGGCGTAGTTCGc                                                                       >  1:61561/1‑141 (MQ=255)
                                                          atgaaAATGCATGTAAGGACACATGTGCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGAAACATATACGACGATGACCGTAATGACGAAAGACAATAAAAATTTTACTAACCGATAGCGTAGTTCGc                                                                       >  2:87739/1‑141 (MQ=255)
                                                                   cATGTAAGGACACATGTGCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGGAACATATACGACGATGACCGTAATGACGAACGACAATAAAATTTTTACTAACCGATGGCGTAGTTCGCTAAAGTGaa                                                              <  2:336494/141‑1 (MQ=255)
                                                                   cATGTAAGGACACATGTGCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGGAACATATACGACGATGACCGTAATGACGAACGACAATAAAATTTTTACTAACCGATGGCGTAGTTCGCTAAAGTGaa                                                              <  2:100166/141‑1 (MQ=255)
                                                                    aTGTAAGGACACATGTGCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGAAACATATACGACGATGACCGTAATGACGAACGACAATAAAATTTTTACTAACCGATGGCGTAGTTCGCTAAAGTGAAc                                                             <  2:403154/141‑1 (MQ=255)
                                                                    aTGTAAGGACACATGTGCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGAAACATATACGACGATGACCGTAATGACGAACGACAATAAAATTTTTACTAACCGATGGCGTAGTTCGCTAAAGTGAAc                                                             <  2:171384/141‑1 (MQ=255)
                                                                       tAAGGACACATGTGCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGAAACATATACGACGATGACTGTGACGACGAACGACAATAGTATTTTAACTAAACGATGACGTAGCTCGCTAAAATGAActag                                                          <  1:12508/141‑4 (MQ=255)
                                                                             cacaTGTGCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGAAACATATACGACGATGACCGTAATGACGAACGACAATAAAATTTTTACTAACCGATGGCGTAGTTCGCTAAAGTGAACCAGTAAGgt                                                    >  1:37280/1‑141 (MQ=255)
                                                                               caTGTGCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGAAACATATACGACGATGACCGTAATGACGAACGACAATAAAATTTTTACTAACCGATGGCGTAGTTCGCTAAAGTGAACCAGTAAGgtgt                                                  >  2:398434/1‑141 (MQ=255)
                                                                                 tgtgCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGAAACATATACGACGATGACCGTAATGACGAACGACAATAAAATTTTTACTAACCGATGGCGTAGTTCGCTAAAGTGAACCAGTAAGGTGTgg                                                >  2:230851/1‑141 (MQ=255)
                                                                                  gtgCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGAAACATATACGACGATGACTGTGACGACGAACGACAATAGTATTTTAACTAAACGATGACGTAGCTCGCTAAAATGAACTAGTAAGGTGTGGt                                               >  1:2890/1‑141 (MQ=255)
                                                                                  gtgCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGAAACATATACGACGATGACCGTAATGACGAACGACAATAAAATTTTTACTAACCGATGGCGTAGTTCGCTAAAGTGAACCAGTAAGGTGTGGt                                               >  2:322484/1‑141 (MQ=255)
                                                                                    gCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGAAACATATACGACGATGACCGTAATGACGAACGACAATAAAATTTTTACTAACCGATGGCGTAGTTCGCTAAAGTGAACCAGTAAGGTGTGGTCa                                             <  1:108898/141‑1 (MQ=255)
                                                                                      aCGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGAAACATATACGACGATGACCGTAATGACGAACGACAATAAAATTTTTACTAACCGATGGCGTAGTTCGCTAAAGTGAACCAGTAAGGTGTGGTCAGt                                           <  1:169103/141‑1 (MQ=255)
                                                                                           tAATATACGTTATGAATGGTGTCATCCACCAAAATGATGAAACATATACGACGATGACCGTAATGACGAACGACAATAAAATTTTTACTAACCGATGGCGTAGTTCGCTAAAGTGAACCAGTAAGGTGTGGTCAGTGCTAt                                      >  1:6331/1‑141 (MQ=255)
                                                                                                         aaTGGTGTCATCCACCAAAATGATGAAACATATACGACGATGACCGTAATGACGAACGACAATAAAATTTTTACTAACCGATGGCGTAGTTCGCTAAAGTGAACCAGTAAGGTGTGGTCAGTGCTATTGCTCTCGCtgttg                        >  1:225987/1‑141 (MQ=255)
                                                                                                         aaTGGTGTCATCCACCAAAATGATGAAACATATACGACGATGACCGTAATGACGAACGACAATAAAATTTTTACTAACCGATGGCGTAGTTCGCTAAAGTGAACCAGTAAGGTGTGGTCAGTGCTATTGCTCTCGCtgttg                        <  2:389298/141‑1 (MQ=255)
                                                                                                                  aTCCACCAAAATGATGAAACATATACGACGATGACCGTAATGACGAACGACAATAAAATTTTTACTAACCGATGGCGTAGTTCGCTAAAGTGAACCAGTAAGGTGTGGTCAGTGCTATTGCTCTCGCTGTTGTTTCGATTc               <  2:315413/141‑1 (MQ=255)
                                                                                                                  aTCCACCAAAATGATGAAACATATACGACGATGACCGTAATGACGAACGACAATAAAATTTTTACTAACCGATGGCGTAGTTCGCTAAAGTGAACCAGTAAGGTGTGGTCAGTGCTATTGCTCTCGCTGTTGTTTCGATTc               >  1:120801/1‑141 (MQ=255)
                                                                                                                   tCCACCAAAATGATGAAACATATACGACGATGACCGTAATGACGAACGACAATAAAATTTTTACTAACCGATGGCGTAGTTCGCTAAAGTGAACCAGTAAGGTGTGGTCAGTGCTATTGCTCTCGCTGTTGTTTCGATTcc              <  1:357668/141‑1 (MQ=255)
                                                                                                                               atgaAACATATACGACGATGACCGTAATGACGAACGACAATAAAATTTTTACTAACCGATGGCGTAGTTCGCTAAAGTGAACCAGTAAGGTGTGGTCAGTGCTATTGCTCTCGCTGTTGTTTCGATTCCTTACTGGGTGTa  >  1:160138/1‑141 (MQ=255)
                                                                                                                                  |                                                                                                                                         
ATGAAACAAAACGCGATGAAAAATATAATCATCACGTATATTTGTATCATTTCTGTGAATGAAAATGCATGTAAGGACACATGTGCACGGGTAATATACGTTATGAATGGTGTCATCCACCAAAATGATGAAACATATACGACGATGACCGTAATGACGAACGACAATAAAATTTTTACTAACCGATGGCGTAGTTCGCTAAAGTGAACCAGTAAGGTGTGGTCAGTGCTATTGCTCTCGCTGTTGTTTCGATTCCTTACTGGGTGTA  >  CP000730/400274‑400541

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