Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,066,392 T→A 7.7% K310N (AAA→AAT) ‡ 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,3920TA 7.7% 134.4 / 4.0 52K310N (AAA→AAT) ‡insH1IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base T (25/23);  new base A (2/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 = 7.71e-01

TCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAGG  >  NC_000913/2066318‑2066537
                                                                          |                                                                                                                                                 
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)
                                                                          |                                                                                                                                                 
TCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAGG  >  NC_000913/2066318‑2066537

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: