Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,631,141 C→G 9.5% L202L (CTC→CTG)  lon‑II → DNA‑binding ATP‑dependent protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,631,1410CG 9.5% 142.3 / 6.9 53L202L (CTC→CTG) lon‑IIDNA‑binding ATP‑dependent protease
Reads supporting (aligned to +/- strand):  ref base C (20/27);  new base G (2/3);  total (23/30)
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.70e-01

CAGCATCGACGAGCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGCTGAA  >  NC_002947/2630994‑2631287
                                                                                                                                                   |                                                                                                                                                  
cAGCATCGACGAGCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGAcctcc                                                                                                                                                   >  7:218/1‑149 (MQ=255)
    aTCGACGAGCCAGGGCGCCTGGTCGACACCATGGACGCGCACATGGCCCAGAAGATCGAGCAGAAGCAGGAAATTCTAGAAATCGTCGACTTGAATACCCGTGTCCAGCACGTCCTGGCGCTGCTGGAAGCAGAAATCCACCTCCTg                                                                                                                                                 >  8:196856/1‑147 (MQ=255)
       gacgaGCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGt                                                                                                                                            <  1:136086/149‑1 (MQ=255)
        acgaGCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTc                                                                                                                                           <  3:162564/149‑1 (MQ=255)
         cgaGCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCg                                                                                                                                          <  8:218/149‑1 (MQ=255)
            gCCAGGGCGCCTGGTCGACACCACGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTGCAGGAGGTCGatt                                                                                                                                       <  7:109575/149‑3 (MQ=255)
               aGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAgc                                                                                                                                    >  2:123831/1‑149 (MQ=255)
               aGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAAACGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTGCAGGAGGTCGatttcg                                                                                                                                    <  1:102798/149‑6 (MQ=255)
                  gcgcCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCa                                                                                                                                 >  3:254678/1‑149 (MQ=255)
                             acacCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCgg                                                                                                                      <  2:109973/149‑1 (MQ=255)
                             acacCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCgg                                                                                                                      <  5:61351/149‑1 (MQ=255)
                                 cATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCTCATTCTCGGCCGGGTc                                                                                                                   >  6:240139/1‑148 (MQ=255)
                                 cATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTc                                                                                                                   >  1:236432/1‑148 (MQ=255)
                                 cATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTc                                                                                                                   >  5:55564/1‑148 (MQ=255)
                                       cgcgCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGt                                                                                                                    <  1:60225/141‑1 (MQ=255)
                                       cgcgCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGt                                                                                                                    >  2:60225/1‑141 (MQ=255)
                                        gcgcACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCGGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCa                                                                                                           <  6:213548/149‑1 (MQ=255)
                                           cacaTGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGAt                                                                                                        <  4:218435/149‑1 (MQ=255)
                                             caTGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATgg                                                                                                      <  2:236432/149‑1 (MQ=255)
                                              aTGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATg                                                                                                       >  8:144596/1‑147 (MQ=255)
                                                  cccTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAgcgc                                                                                                 <  7:195220/149‑1 (MQ=255)
                                                   ccTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCCGCCgg                                                                                                                      <  6:322238/127‑1 (MQ=255)
                                                   ccTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCCCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCCGCCgg                                                                                                                      >  5:322238/1‑127 (MQ=255)
                                                         gATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGcagccagc                                                                                          <  8:256/149‑1 (MQ=255)
                                                            cGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAgcg                                                                                         >  6:274831/1‑147 (MQ=255)
                                                            cGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAgcg                                                                                         <  5:274831/147‑1 (MQ=255)
                                                                cAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCgg                                                                                                                      <  6:69264/114‑1 (MQ=255)
                                                                cAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCgg                                                                                                                      >  5:69264/1‑114 (MQ=255)
                                                                cAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGta                                                                                   >  6:263375/1‑149 (MQ=255)
                                                                cAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACCTGCTCGCGCTGCAGCATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATCGAGCGCAGCCAGAGCGAGt                                                                                    >  7:140916/1‑148 (MQ=255)
                                                                       aGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCt                                                                             >  4:22122/1‑148 (MQ=255)
                                                                         gAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAgcgc                                                                                        <  6:34142/135‑1 (MQ=255)
                                                                         gAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAgcgc                                                                                        >  5:34142/1‑135 (MQ=255)
                                                                          aaaTTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAAt                                                                         <  1:123831/149‑1 (MQ=255)
                                                                             ttCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAg                                                                      <  6:26756/149‑1 (MQ=255)
                                                                              tCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAg                                                                      >  1:207383/1‑148 (MQ=255)
                                                                                 tGAAATCGTCGACTTGACTACCCTTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCTAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAg                                                                   <  7:196856/148‑1 (MQ=255)
                                                                                        gtcgACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAgg                                                            <  7:144596/148‑1 (MQ=255)
                                                                                                     cccGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCagaaagaa                                              >  1:227128/1‑149 (MQ=255)
                                                                                                      ccGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAAt                                                                         >  8:241442/1‑121 (MQ=255)
                                                                                                      ccGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAAt                                                                         <  7:241442/121‑1 (MQ=255)
                                                                                                      ccGTGTCGAGCACGTGCTGGCGATGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCTGGTCAAGAAGCAGATGGAGCGCAGCCAGCGAGAGTACTACCTGAATGAGCAGATGAAGGCCAtt                                                       <  8:207940/139‑1 (MQ=255)
                                                                                                      cagtgtCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTACTGCAGGTCGAAAAGCGCATTCGCGGCCTGGGAAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCAtt                                                       >  7:207940/3‑139 (MQ=255)
                                                                                                         tgtCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTg                                          <  3:22122/149‑1 (MQ=255)
                                                                                                               gCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATg                                    <  4:254678/149‑1 (MQ=255)
                                                                                                                cACGTGCTGGCGCTTCTGGATGCCTTAATCGACCTCCTGCAGGTCGAAAAGCGCATTCTCGGCCTGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATgg                                   <  8:140916/149‑1 (MQ=255)
                                                                                                                cACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATgg                                   <  5:263375/149‑1 (MQ=255)
                                                                                                                   gTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGa                                <  4:249746/149‑1 (MQ=255)
                                                                                                                                   aTGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGaa                <  2:227128/149‑1 (MQ=255)
                                                                                                                                       gcttttCGACCTGCAGGAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTgg                                         <  3:84941/115‑1 (MQ=255)
                                                                                                                                       gcttttCGACCTGCAGGAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTgg                                         >  4:84941/6‑120 (MQ=255)
                                                                                                                                       cGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTc             >  8:269544/1‑148 (MQ=255)
                                                                                                                                         ttttCGACCTGCAGGAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCg            >  2:207863/4‑147 (MQ=255)
                                                                                                                                             cGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTAACTGAATGAGCAGATCAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAAAGAAGTCGAagag       >  2:291859/1‑148 (MQ=255)
                                                                                                                                                 ctcctGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGCTGa   >  3:162143/1‑148 (MQ=255)
                                                                                                                                                  tcctGCAGGTCGAAAAGCGCATTCGCCGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGCTGaa  >  6:158746/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CAGCATCGACGAGCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGCTGAA  >  NC_002947/2630994‑2631287

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑