Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,380,709 G→A 9.1% Q127Q (CAG→CAA) ‡ tnpT‑I → cointegrate resolution protein T

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,380,7090GA 9.1% 152.7 / 4.6 55Q127Q (CAG→CAA) ‡tnpT‑Icointegrate resolution protein T
Reads supporting (aligned to +/- strand):  ref base G (24/26);  new base A (3/2);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01

GCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCAGGAGGTGCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGAAAAGCATCAA  >  NC_002947/3380562‑3380856
                                                                                                                                                   |                                                                                                                                                   
gCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAgg                                                                                                                                                    >  4:134027/1‑149 (MQ=255)
           aTGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCCTGAAGCTGTTccgtga                                                                                                                                         <  4:187826/149‑6 (MQ=255)
              cgcAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAgcg                                                                                                                                      >  7:818/1‑149 (MQ=255)
               gcAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAgcgc                                                                                                                                     >  4:43979/1‑149 (MQ=255)
                 aGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTa                                                                                                                                   <  4:196341/149‑1 (MQ=255)
                      aCCGTGTCCCAGGAGCGCGAACAACTACAACGCGAACGGCTGTATTACCAGAAACAAGCTCGTCTTGCCGAAAGCCGAATCCAGCAGTTGTAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTAcagca                                                                                                                              <  3:42218/149‑1 (MQ=255)
                      aCCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTAcagca                                                                                                                              <  3:134027/149‑1 (MQ=255)
                      aCCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTAcagca                                                                                                                              <  6:259861/149‑1 (MQ=255)
                       ccGGGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTAcagcag                                                                                                                             <  8:279370/149‑1 (MQ=255)
                        cGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAgg                                                                                                                            <  5:251020/149‑1 (MQ=255)
                                    gcgcGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAAc                                                                                                                >  4:143229/1‑149 (MQ=255)
                                     cgcgAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACg                                                                                                               >  8:25139/1‑149 (MQ=255)
                                        gAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGcc                                                                                                             >  1:133946/1‑148 (MQ=255)
                                          acaacTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAAcagcag                                                                                                                    >  7:210544/1‑139 (MQ=255)
                                          acaacTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAAcagcag                                                                                                                    <  8:210544/139‑1 (MQ=255)
                                            aacTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCaa                                                                                                           >  2:243292/1‑146 (MQ=255)
                                                  aaCGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAg                                                                                                  >  8:278948/1‑149 (MQ=255)
                                                       gAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGtt                                                                                               >  5:300809/1‑147 (MQ=255)
                                                       gAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTga                                                                                             <  7:81740/149‑1 (MQ=255)
                                                                      cAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTAcag                                                                                                                                <  3:38496/99‑1 (MQ=255)
                                                                      cAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTAcag                                                                                                                                >  4:38496/1‑99 (MQ=255)
                                                                      cAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGgc                                                                              >  2:100660/1‑149 (MQ=255)
                                                                      cAGAACCAAGCTCGACTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTAcagca                                                                                                                              >  5:61329/1‑101 (MQ=255)
                                                                      cAGAACCAAGCTCGACTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTAcagca                                                                                                                              <  6:61329/101‑1 (MQ=255)
                                                                       aGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGTTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGgcg                                                                             >  7:136247/1‑149 (MQ=255)
                                                                         aaCCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCa                                                                           <  7:25139/149‑1 (MQ=255)
                                                                         aaCCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCa                                                                           <  8:818/149‑1 (MQ=255)
                                                                            cAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAGGc                                                                        <  3:43979/149‑1 (MQ=255)
                                                                               gCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGcc                                                                                      <  1:111685/132‑1 (MQ=255)
                                                                               gCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGcc                                                                                      >  2:111685/1‑132 (MQ=255)
                                                                                       ggCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGAt                                                              >  6:297/1‑148 (MQ=255)
                                                                                                               ggAAAGTCAGAGCAGCGGGCTCACTGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCaggag                                                        <  4:126472/130‑1 (MQ=255)
                                                                                                               ggAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCaggag                                                        >  3:126472/1‑130 (MQ=255)
                                                                                                               ggAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCAGGAGGTGCGCCTGCAGGATcgc                                      <  5:297/148‑1 (MQ=255)
                                                                                                                gAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCAGGAGGTGCGCCTGCAGGATCGCGa                                    <  3:143229/149‑1 (MQ=255)
                                                                                                                 aaaGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCAGGAGGTGCGCCTGCAGGATCGCGAc                                   >  7:273884/1‑149 (MQ=255)
                                                                                                                  aaGTCAGAGCAGCGGTCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCTCTTGGCGCAGGCTAACGGCGATCAGGAGGTGCACCTGCAGGATCGCGACa                                  <  1:243292/149‑1 (MQ=255)
                                                                                                                      cAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCAGGAGGTGCGCCTGCAGGATCGCGACAGcc                               >  2:51405/1‑148 (MQ=255)
                                                                                                                         agcagcGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCAGGAGGTgcgc                                                  <  7:36442/126‑1 (MQ=255)
                                                                                                                         agcagcGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCAGGAGGTgcgc                                                  >  8:36442/1‑126 (MQ=255)
                                                                                                                          gcagcGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCAGGAGGTGCGTCTGCAGGATCGCGACAGCCAAAtt                          <  6:300809/149‑1 (MQ=255)
                                                                                                                          gcagcGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCAGGAGGTGCGCCTGCAGGATCGCGACAGCCAAAtt                          <  1:51405/149‑1 (MQ=255)
                                                                                                                             gcGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAgg                                                                         >  4:42686/1‑99 (MQ=255)
                                                                                                                             gcGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAgg                                                                         <  3:42686/99‑1 (MQ=255)
                                                                                                                                     aCGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCAGGAGGTGCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTgga               <  8:136247/149‑1 (MQ=255)
                                                                                                                                     aCGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCAGGAGGTGCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTgga               <  2:133946/149‑1 (MQ=255)
                                                                                                                                     aCGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCAGGAGGTGCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTgg                <  1:100660/148‑1 (MQ=255)
                                                                                                                                     aCGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCAGGAGGTGCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTg                 >  8:133856/1‑147 (MQ=255)
                                                                                                                                      ctgAACAGCCTGAAGCTGTTCAGCAAGCGCTACAGCAGGTACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCAGGAGGTGCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTggag              >  5:124668/3‑149 (MQ=255)
                                                                                                                                      ctgAACAGCCTGAAGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCAGGAGGTGCGCCTGCAGGATCGCGACa                                  >  2:255873/3‑129 (MQ=255)
                                                                                                                                      ctgAACAGCCTGAAGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCAGGAGGTGCGCCTGCAGGATCGCGACa                                  <  1:255873/127‑1 (MQ=255)
                                                                                                                                       tgAACAGCCTGAAGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGGTGCAGGCTAACGGCGATCAGGAGGTGCGCCTGCAGGATCGCGCCAGCCAAATTCGCTCGCTggag              >  4:258423/2‑148 (MQ=255)
                                                                                                                                          aCAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCAGGAGGTGCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGaa           >  1:206934/1‑148 (MQ=255)
                                                                                                                                               ttCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCAGGAGGTGCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGAAAAGCAt     >  4:219314/1‑149 (MQ=255)
                                                                                                                                                tCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCAGGAGGTGCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGAAAAGCAt     >  5:183331/1‑148 (MQ=255)
                                                                                                                                                  aGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCAGGAGGTGCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGAAAAGCATCaa  <  7:133856/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCTAGAGGAGGATGCGCAGGCCACCGTGGCCCAGGAGCGCGAACAACTACAACGCGAACGGCTCGATTACCAGAACCAAGCTCGGCTGGCCGAAAGCCGAATCCAGCAGTTGGAAAGTCAGAGCAGCGGGCTCACGGAACAGCTTCAGGCTGTTCAGCAAGCGCTACAGCAGGAACAGCAGCAACGCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCAGGAGGTGCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGAAAAGCATCAA  >  NC_002947/3380562‑3380856

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: