Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,631,145 C→G 7.7% Q204E (CAG→GAG)  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,1450CG 7.7% 158.0 / 3.5 52Q204E (CAG→GAG) lon‑IIDNA‑binding ATP‑dependent protease
Reads supporting (aligned to +/- strand):  ref base C (24/24);  new base G (2/2);  total (26/26)
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.50e-01

CATCGACGAGCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGCTGAAAA  >  NC_002947/2630997‑2631289
                                                                                                                                                    |                                                                                                                                                 
cATCGACGAGCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGc                                                                                                                                                   <  4:201750/149‑1 (MQ=255)
             gggCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAgcg                                                                                                                                      >  5:108281/1‑149 (MQ=255)
             gggCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAgcg                                                                                                                                      >  6:170083/1‑149 (MQ=255)
                      gTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCgg                                                                                                                             <  5:170083/149‑1 (MQ=255)
                         gACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCg                                                                                                                          <  7:169215/149‑1 (MQ=255)
                           cacCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCgg‑g                                                                                                                       >  5:332416/1‑149 (MQ=255)
                              cATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTc                                                                                                                     >  1:175234/1‑148 (MQ=255)
                                                    aaGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGcagc                                                                                                >  1:230096/1‑147 (MQ=255)
                                                      gATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGa                                                                                                           >  7:21249/1‑134 (MQ=255)
                                                      gATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGa                                                                                                           <  8:21249/134‑1 (MQ=255)
                                                               gAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGtac                                                                                    >  4:250117/1‑148 (MQ=255)
                                                                 aGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGtactac                                                                                 <  2:230096/149‑1 (MQ=255)
                                                                  gCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTAcc                                                                                <  6:108281/149‑1 (MQ=255)
                                                                    aGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCt                                                                               >  5:280389/1‑148 (MQ=255)
                                                                         aTTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTc                                                                                                                     >  4:275925/1‑105 (MQ=255)
                                                                         aTTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTc                                                                                                                     <  3:275925/105‑1 (MQ=255)
                                                                          ttCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAg                                                                        >  5:179909/1‑149 (MQ=255)
                                                                               gAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCa                                                                      <  6:280389/146‑1 (MQ=255)
                                                                                 aaTCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTCCTGGATGCCGAAAACGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTCAATGAGCAGATg                                                                  >  3:70044/1‑148 (MQ=255)
                                                                                   tcgtcgACTTGACTACCCGTGTCGAGCCCGTGCGGGCGCGGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAg                                                               <  6:179909/149‑1 (MQ=255)
                                                                                      tcgACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGcagcc                                                                                               >  3:76701/1‑114 (MQ=255)
                                                                                      tcgACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGcagcc                                                                                               <  4:76701/114‑1 (MQ=255)
                                                                                          cTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTc                                                        <  6:332416/149‑1 (MQ=255)
                                                                                           ttGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAgc                                                                                                     <  8:154108/103‑1 (MQ=255)
                                                                                                   ccGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAAc                                               >  8:187066/1‑149 (MQ=255)
                                                                                                   ccGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAAc                                               <  5:132122/149‑1 (MQ=255)
                                                                                                        tCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGc                                                             >  7:39921/1‑130 (MQ=255)
                                                                                                        tCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGc                                                             <  8:39921/130‑1 (MQ=255)
                                                                                                        tCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCAtt                                                         >  4:207222/1‑134 (MQ=255)
                                                                                                        tCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCAtt                                                         <  3:207222/134‑1 (MQ=255)
                                                                                                          gAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATg                                                                                                         >  7:310074/1‑84 (MQ=255)
                                                                                                          gAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATg                                                                                                         <  8:310074/84‑1 (MQ=255)
                                                                                                          gAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATg                                                                                                         >  1:261083/1‑84 (MQ=255)
                                                                                                          gAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATg                                                                                                         <  2:261083/84‑1 (MQ=255)
                                                                                                          gAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCAtt                                                         <  5:236991/132‑1 (MQ=255)
                                                                                                          gAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCAtt                                                         >  6:236991/1‑132 (MQ=255)
                                                                                                          gAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGa                                        >  2:274954/1‑149 (MQ=255)
                                                                                                            gCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTACTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGAt                                       <  2:175234/148‑1 (MQ=255)
                                                                                                                gTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGa                                  >  4:210307/1‑149 (MQ=255)
                                                                                                                    tGGCGCTGGTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCTGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGATAGAACTTGGCGATGGTGACGaa                              <  4:70044/149‑1 (MQ=255)
                                                                                                                      gcgcTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAg                             <  3:210307/148‑1 (MQ=255)
                                                                                                                              ggATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCAtt                                                         >  6:334556/1‑112 (MQ=255)
                                                                                                                              ggATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCAtt                                                         <  5:334556/112‑1 (MQ=255)
                                                                                                                                  gCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGAGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCagaaagaa                                                >  1:73767/1‑118 (MQ=255)
                                                                                                                                  gCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCagaaagaa                                                <  2:73767/117‑1 (MQ=255)
                                                                                                                                    gcttttCGACCTGCAGGAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATg                                      <  1:52766/120‑1 (MQ=255)
                                                                                                                                    gcttttCGACCTGCAGGAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATg                                      >  2:52766/6‑125 (MQ=255)
                                                                                                                                    gcttttCGACCTGCAGGAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTc               >  4:83729/6‑148 (MQ=255)
                                                                                                                                    cGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATg                                      >  8:21392/1‑125 (MQ=255)
                                                                                                                                           gACCTGCAGGAGGTCTAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGCt       <  3:83729/149‑1 (MQ=255)
                                                                                                                                           gACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGc        <  1:274954/148‑1 (MQ=255)
                                                                                                                                           gACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGCt       >  1:142440/1‑149 (MQ=255)
                                                                                                                                                 ctGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GGCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGCTGaaaa  <  5:193136/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CATCGACGAGCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGCTGAAAA  >  NC_002947/2630997‑2631289

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