Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,066,394 T→C 9.6% K310E (AAA→GAA) ‡ insH1 ← IS5 transposase and trans‑activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,066,3940TC 9.6% 129.3 / 6.3 52K310E (AAA→GAA) ‡insH1IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base T (24/23);  new base C (3/2);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.48e-01

TCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAGGTTC  >  NC_000913/2066318‑2066540
                                                                            |                                                                                                                                                  
tcCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCtt                                                                            >  2:103359/1‑149 (MQ=31)
  ccACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTgg                                                                          >  2:450654/1‑149 (MQ=32)
    aCTGACTTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGGCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGcc                                                                        <  3:72877/149‑1 (MQ=17)
         cGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGAt                                                                   <  1:103359/149‑1 (MQ=32)
            aTCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCt                                                                >  7:504002/1‑149 (MQ=32)
             tCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTg                                                               <  2:44145/149‑1 (MQ=32)
              cATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCGTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTGGGCCCGGATGCTg                                                               <  4:220660/148‑1 (MQ=17)
                tttGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCt                                                            >  8:195076/1‑149 (MQ=32)
                    gTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTc                                                         <  2:92745/148‑1 (MQ=32)
                     tCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCAt                                                       >  7:336515/1‑149 (MQ=32)
                     tCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCAt                                                       <  7:14194/149‑1 (MQ=32)
                        gCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCACCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGCTGCGAAATGGGTGCTCCACCTTGGCGCGGATGCTGGCTTTCATgga                                                    >  6:454402/1‑147 (MQ=12)
                        gCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTa                                                    >  3:240124/1‑149 (MQ=32)
                           cGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTc                                                         <  4:408096/141‑1 (MQ=32)
                           cGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTc                                                         >  3:408096/1‑141 (MQ=32)
                                   gTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGa                                         >  1:79014/1‑149 (MQ=32)
                                     tGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCATTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTTGCCAGGATGCTGGCTTTCATGTATTCGATGTTGAt                                        <  5:137060/148‑1 (MQ=12)
                                         cAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCg                                   >  2:65662/1‑149 (MQ=32)
                                               gAATAACATCGCCAGTTGGTTATCGTTATCCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCAt                                                       <  2:147896/123‑1 (MQ=17)
                                               gAATAACATCGCCAGTTGGTTATCGTTATCCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCAt                                                       <  2:154587/123‑1 (MQ=17)
                                               gAATAACATCGCCAGTTGGTTATCGTTATCCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCAt                                                       >  1:147896/1‑123 (MQ=17)
                                               gAATAACATCGCCAGTTGGTTATCGTTATCCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCAt                                                       >  1:154587/1‑123 (MQ=17)
                                               gAATAACATCGCCACTTGGGTATCGTTTTTCACCAACCCCTTGTATCTGGCTTTCACGAGGCCGAACTGTCGCTTGATGAAGCGCAATGCGTGCTCCACCTTGGCCCGGATCCTCGCTTTCATGTATTCGATGGTGATGGCCGGTTTGt                             <  2:101228/149‑1 (MQ=11)
                                                aaTAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGt                             <  8:37725/148‑1 (MQ=32)
                                                 aTAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTc                           >  1:220642/1‑149 (MQ=32)
                                                  tAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTg                                                               >  2:370819/1‑112 (MQ=32)
                                                  tAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTg                                                               <  1:370819/112‑1 (MQ=32)
                                                  tAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCt                          <  8:504002/149‑1 (MQ=32)
                                                    aCATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCtt                         <  3:427747/148‑1 (MQ=32)
                                                          ccAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTTGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTTCTCCTCCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGa                  <  2:439958/149‑1 (MQ=12)
                                                          ccAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGa                  >  6:157760/1‑149 (MQ=32)
                                                           cAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGAt                 <  1:450654/149‑1 (MQ=32)
                                                           cAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGAt                 >  7:142201/1‑149 (MQ=32)
                                                                 gTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGt            <  8:142201/148‑1 (MQ=32)
                                                                  ttATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGttt          <  4:240124/149‑1 (MQ=32)
                                                                  ttATCGTTTTCCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGttt          >  8:445506/1‑149 (MQ=17)
                                                                   tATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATg                                             >  2:259639/1‑113 (MQ=32)
                                                                   tATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATg                                             <  1:259639/113‑1 (MQ=32)
                                                                   tATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGttt          >  3:175469/1‑148 (MQ=32)
                                                                   tATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGttt          >  1:369373/1‑148 (MQ=32)
                                                                    aTCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTc         >  5:284386/1‑148 (MQ=32)
                                                                    aTCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTc         >  2:370921/1‑148 (MQ=32)
                                                                    aTCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTc         >  5:149893/1‑148 (MQ=32)
                                                                    aTCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTAGGCCCGGATGCTGGCTTTCATGAATTCGATGTTGATGGCCGTTTTGTTCTTGAGTGGATGCTGTTTc         >  2:4592/1‑148 (MQ=12)
                                                                    aTCGTTTTTCAGCAACCCCATGTATCAGGCTTTCACGAAGCCGAAATCTCCCATGAAGATCCGAAATGGGTGCTCCACCTTGCCCCTGAAGCTGGCTTACATGTATTCGAAGTTGATTCCCGTTTTGTTCTTGCGTTCATGCTGTTTc         >  4:25167/1‑148 (MQ=11)
                                                                     tCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTc         >  8:456891/1‑147 (MQ=32)
                                                                     tCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCaa       <  2:220642/149‑1 (MQ=32)
                                                                     tCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCaa       <  8:393821/149‑1 (MQ=32)
                                                                     tCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCaa       <  7:456891/149‑1 (MQ=32)
                                                                       gTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAgg     >  7:7202/1‑149 (MQ=32)
                                                                       gTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAgg     <  3:25167/149‑1 (MQ=32)
                                                                       gTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAg      <  1:361669/148‑1 (MQ=32)
                                                                        tttttCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTg                                                               >  2:100165/1‑90 (MQ=33)
                                                                        tttttCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTg                                                               <  1:100165/90‑1 (MQ=33)
                                                                           ttCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAGGTTc  <  8:336515/148‑1 (MQ=32)
                                                                            |                                                                                                                                                  
TCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAGGTTC  >  NC_000913/2066318‑2066540

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