Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 200,642 C→T 11.6% A54V (GCC→GTC)  skp → periplasmic chaperone Skp

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913200,6420CT11.6% 69.7 / 2.4 35A54V (GCC→GTC) skpperiplasmic chaperone Skp
Reads supporting (aligned to +/- strand):  ref base C (17/13);  new base T (2/2);  total (20/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

GCTGCAGGTCTCGGTTTAGCACTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCAGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTAAGCTGGAAAAAGACGTGATGGCTCAGCGCCAGACTTTTGCTCAGAAAGC  >  NC_000913/200500‑200774
                                                                                                                                               |                                                                                                                                    
gCTGCAGGTCTCGGTTTAGCTCTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCAGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGc                                                                                                                                      <  1:309283/143‑1 (MQ=255)
  ggCAGGTCTCGGTCTAGCACTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCAGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCa                                                                                                                                    <  2:350489/142‑1 (MQ=255)
        tctcGGTTTAGCACTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCAGGTAGAGCCGAACACCGGTGTCTCTAACACGCTGGAAACGGAGCTCAAGGGCCGTGCCGGCGacc                                                                                                                              >  2:508502/1‑141 (MQ=255)
          tcGGTGTAGCACTGGCGTCTTCTGGTCAGGCGGCTGACAAGATGGCGATCGGCGACATGTGCAGGCTGGTCCGGCAGGGAGCGCAGAAAACCGGTGTTGCTAACACGGTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTg                                                                                                                            <  2:407963/143‑1 (MQ=255)
              tttAGCGCTGGCAACCTCTGCGCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCAGGTAGCGCAGAAACCCGGTGTTTCTAACACGCTGGCAAATGATTTCAATGGCCGTGCCAGCGAACTGCAGc                                                                                                                        <  2:61543/143‑1 (MQ=255)
                   cACTGGCAACTTCTGCTCAGGCTGCTCCGCAGAAGTCAATCGTCAACATGGGCAGCCTGTTCCAGCAGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAAGGAGTTCAAAGGCCGTGCCAGCGAACCGCCGCGTATg                                                                                                                   >  2:604059/1‑143 (MQ=255)
                      tGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCAGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTAGGGaa                                                                                                                >  1:244172/1‑143 (MQ=255)
                            cTTGTGCTCAGGCGGCTGACAAAAGTGCAATCGTCAACATGGGCAGCCTGGTCCAGCAGGTAGCGCAGAAAACCGGTGTTGCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGAt                                                                                                          <  1:297654/143‑1 (MQ=255)
                                  cTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCAGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGCAg                                                                                                    >  1:490362/1‑143 (MQ=255)
                                     cggcggcGTACAAAAGTGCAATCGTCAGCATGGGCAGCCTGTTCCAGTAGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAGGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGCAGGCt                                                                                                 <  2:39814/142‑1 (MQ=255)
                                       gcggcTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCAGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGTGAGCGAACTGCAGCGTATGGAAACCGATCTGCAGGCTaa                                                                                               >  1:23009/1‑143 (MQ=255)
                                       gcggcTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCAGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGTCAGCGAACTGCAGCGTATGGAAACCGATCTGCAGGCTac                                                                                               >  1:23057/1‑142 (MQ=255)
                                                 aaaaTTGCAATCGTCAACATGGGCAGCCTGTTCCAGCAGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGCAGGCTAAAATGaaaaa                                                                                      >  1:322351/1‑143 (MQ=255)
                                                 aaaaTTGCAATCGTCAACATGGGCAGCCTGTTCCAGCAGGTAGCGCAGAAAACCGGTGTATCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACGAATCTGGAGGCTGAAAAGaaaaa                                                                                      >  1:777355/1‑143 (MQ=255)
                                                     ttGCAATCGTCAACATGGGCAGCCTGTTCCAGCAGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTg                                                                                  <  2:471052/143‑1 (MQ=255)
                                                     ttGCAATCGTCAACATGGGCAGCCTGGTCCAGCAGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTg                                                                                  <  1:278278/143‑1 (MQ=255)
                                                              tCAACATGGGCAGCCTGTTCCAGCAGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCCGTCCagg                                                                         >  2:283677/1‑141 (MQ=255)
                                                                        cAGCCTGTTCCAGCAGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGCAGGCTCAAATGAAAACGCTGCAGTCCATCACAGCGGGCa                                                               >  2:530014/1‑143 (MQ=255)
                                                                        cAGCCTGTACCAGCAGGGAGCGCAGAAAACCGGTGTTTCTAAAACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCa                                                               <  1:455514/143‑1 (MQ=255)
                                                                             tGTTCTAGCAGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGTCAGCGAt                                                          <  1:499431/143‑1 (MQ=255)
                                                                              gTTCCAGCAGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGCAGGCAAAAATGAAAAAGATGCAGGCCAAGGAGGGGGGGAGCGacc                                                         >  2:4643/1‑141 (MQ=255)
                                                                                 ccAGCAGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGAACTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCa                                                      >  2:353364/1‑143 (MQ=255)
                                                                                      agGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTAAg                                                 >  1:362238/1‑143 (MQ=255)
                                                                                        gTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTAAGCt                                               >  1:693637/1‑143 (MQ=255)
                                                                                          aGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTAAGCTgg                                             <  1:47973/143‑1 (MQ=255)
                                                                                            cgcAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCACGAAAGCGGGCAGCGATCGCACTACCCTGGaa                                           >  2:410541/1‑143 (MQ=255)
                                                                                                               tAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCTATCGCACTAAGCTGGAAAAAGACGTGATGGCTCAgc                        <  1:235978/143‑1 (MQ=255)
                                                                                                                   acGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTAAGCTGGAAAAAGACGTGATGGCTCAGCGCCa                    <  1:698496/143‑1 (MQ=255)
                                                                                                                     ggtGGAAAATGAGTTCAAAGGCCGTGTCAGCGAACTGCAGCGTATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTAAGCTGGAAAAAGACGTGATGGCTCAGCGCCAGa                  <  2:23057/141‑1 (MQ=255)
                                                                                                                     gCTGGAAAATGAGTTCAAAGGCCGTGTCAGCGAACTGCAGCGTATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTAAGCTGGAAAAAGACGTGATGGCTCAGCGCCAGa                  <  2:23009/143‑1 (MQ=255)
                                                                                                                        ggAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTAAGCTGGAAAAAGACGTGATGGCTCAGCGCCAGACtt               >  1:366084/1‑143 (MQ=255)
                                                                                                                         gAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCAATGAAAGCGGGCAGCGATCGCACTAAGCTGGAAACAGGCGGGATGACTCAGCGCCCGACttt              >  2:277265/1‑143 (MQ=255)
                                                                                                                         gAAAATGAGTTCAAAGGCCGTGACAGCGAACTGCAGCGTAAGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGAACTAAGCAGGAAGAAGACGTGAGGGCTCGACGCCAGagctt              >  2:656193/1‑139 (MQ=255)
                                                                                                                                 gTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGCAGGCTAAACTGAAAACGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTCAGCTGGAAAAAGACGGGATGGCTCAGCGCCAGGCGTTTGCTCAGa      >  1:562315/1‑143 (MQ=255)
                                                                                                                                   gcagcggCCGTGCCAGCGAACGGCAGCGTATGGAAGCCGATCTGCAGGCGAAAATGAAGAAGCTGCAGTCCATGAAAGCCGGCAGCGATCGCACTAAGCTGGAAAAAGACGTGATGGCTCAGCGCCAGACTTTTGCTCAGaaa    <  2:650130/138‑1 (MQ=255)
                                                                                                                                     aaaGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTGCATGAAAGCGGGCAGCGATCGCACTAAGCTGGAAAAAGACGTGATGGCTGAGCGCCAGCCTTTTGCTCAGcaacg  >  1:79790/1‑141 (MQ=255)
                                                                                                                                               |                                                                                                                                    
GCTGCAGGTCTCGGTTTAGCACTGGCAACTTCTGCTCAGGCGGCTGACAAAATTGCAATCGTCAACATGGGCAGCCTGTTCCAGCAGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTAAGCTGGAAAAAGACGTGATGGCTCAGCGCCAGACTTTTGCTCAGAAAGC  >  NC_000913/200500‑200774

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