Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,379,626 T→C 7.3% H189R (CAT→CGT)  tnpS ← cointegrate resolution protein S

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,379,6260TC 7.3% 157.5 / 2.4 55H189R (CAT→CGT) tnpScointegrate resolution protein S
Reads supporting (aligned to +/- strand):  ref base T (22/29);  new base C (2/2);  total (24/31)
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.52e-01

TGATCAGGCTGCCGGGATGCAAGGCCTCGTCGCTCAGGTGCCCCCAGCGATCCAGGCGACAGAACACTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAGGCTCATGCCCGCACCGGCTTCGGCCTGGATATGTTCGACCTGCAGCCGACAGAGTTCGTCACTACGAAAACCCCGCCAGAAGCCGATCAGCAGCAACG  >  NC_002947/3379477‑3379773
                                                                                                                                                      |                                                                                                                                                   
tGATCAGGCTGCCGGGATGCAAGGCCTCGTCGCTCAGGTGCCCCCAGCGATCCAGGCGACAGAACACTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAAt                                                                                                                                                     >  1:462579/1‑150 (MQ=255)
    cAGGCTGCCGTGATTCAAGTCCTCGTCGCTCAGGTGCCCCCAGCGATCCAGACGTGCCCAACAATTCGTACGGTTAATGCCGGCGGTACTGTGCCAGCCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGt                                                                                                                                                 <  1:252372/150‑1 (MQ=255)
       gCTGCCGGGATGCAAGGCCTCGTCGCTCAGGTGCCCCCAGCGATCCAGGCGACAGAACACTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCGCATCAAGGCTGGTACGTTATGAggcat                                                                                                                                              >  2:906584/1‑146 (MQ=255)
            cGGGATGCAAGGCCTCGTCGCTCAGGTGCCCCCAGCGATCCAGGCGACAGAACACTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGa                                                                                                                                         <  2:278734/150‑1 (MQ=255)
            cGGGATGCAAGGCCTCGTCGCTCAGGTGCCCCCAGCGATCCAGGCGACAGAACACTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAACGAGTCGTGCCGa                                                                                                                                         >  2:579829/1‑150 (MQ=255)
                            gtcgCTCAGGTGCCCCCAGCGATCCAGGCGACAGAACACTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCAcc                                                                                                                         <  1:906584/150‑1 (MQ=255)
                                 tCAGGTGCCCCCAGCGATCCAGGCGACAGAACACTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTAc                                                                                                                    >  1:133422/1‑150 (MQ=255)
                                 tCAGGTGCCCCCAGCGATCCAGGCGACAGAACACTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTAc                                                                                                                    >  1:133657/1‑150 (MQ=255)
                                      tGCCCCCAGCGATCCAGGCGACAGAACACTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTtgt                                                                                                               <  2:58022/150‑1 (MQ=255)
                                         ccccAGCGATCCAGGCGACAGAACACTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCt                                                                                                                               >  1:790257/1‑131 (MQ=255)
                                         ccccAGCGATCCAGGCGACAGAACACTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCt                                                                                                                               <  2:790257/131‑1 (MQ=255)
                                          cccAGCGATCCAGGCGACAGAACACTTCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCAccc                                                                                                                        <  1:804211/137‑1 (MQ=255)
                                          cccAGCGATCCAGGCGACAGAACACTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCAccc                                                                                                                        >  2:804211/1‑137 (MQ=255)
                                          cccAGCGATCCAGGCGACAGAACACTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCa                                                                                                           >  2:871883/1‑150 (MQ=255)
                                             aGCGATCCAGGCGACAGAACACTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGa                                                                                                        <  1:147763/150‑1 (MQ=255)
                                              gCGATCCAGGCGACAGAACACTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGaa                                                                                                       >  1:370679/1‑150 (MQ=255)
                                               cGATCCAGGCGACAGAACACTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAACGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAAc                                                                                                      <  1:579829/150‑1 (MQ=255)
                                                gATCCAGGCGACAGAACACTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTtgtg                                                                                                              >  2:752143/1‑141 (MQ=255)
                                                gATCCAGGCGACAGAACACTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTtgtg                                                                                                              <  1:752143/141‑1 (MQ=255)
                                                   ccAGGCGACAGAACACTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGTCCATAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAGGc                                                                                                  <  1:747580/150‑1 (MQ=255)
                                                         gACAGAACACTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAGGCTCATGc                                                                                            <  1:967607/150‑1 (MQ=255)
                                                         gACAGAACACTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAACGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAGGCTCATGc                                                                                            >  2:752463/1‑150 (MQ=255)
                                                             aGAACACTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTc                                                                                                                                                <  2:100896/95‑1 (MQ=255)
                                                             aGAACACTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTc                                                                                                                                                >  1:100896/1‑95 (MQ=255)
                                                             aGAACACTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAgg                                                                                                   >  2:1043255/1‑140 (MQ=255)
                                                             aGAACACTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAgg                                                                                                   <  1:1043255/140‑1 (MQ=255)
                                                                acacTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGa                                                                                                                                         >  2:713439/1‑99 (MQ=255)
                                                                acacTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGa                                                                                                                                         <  1:713439/99‑1 (MQ=255)
                                                                  acTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAg                                                                                                          <  2:72042/128‑1 (MQ=255)
                                                                  acTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAg                                                                                                          >  1:72042/1‑128 (MQ=255)
                                                                  acTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAGGCTCATGCCCGCACCg                                                                                    <  2:53392/150‑1 (MQ=255)
                                                                       accacGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAGGCTCATGCCCGCACCGGCTTc                                                                               <  2:772504/150‑1 (MQ=255)
                                                                              gCAATGCGGGCGGCACTGATCCAGTCGAGATAAGCGTGCACCGGGCAGAGGAACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAGGCTCATGCCCGCACCGGCTTCGGCCTgg                                                                        <  1:828346/150‑1 (MQ=255)
                                                                                   gCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAGGCTCATGGCGGCACCGTCTTCGGCCTGGATATg                                                                   >  2:9977/1‑150 (MQ=255)
                                                                                      ggcggcACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAGGCTCATGCCCGCACCTGTTTCGGCCTGGATATGTTc                                                                <  2:462579/150‑1 (MQ=255)
                                                                                              tGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACCGGCTCCTGACCGCACCGGCTTCGGCGTGGAAATGTTCGACCTGcc                                                        >  2:154957/1‑149 (MQ=255)
                                                                                              tGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAGGCTCATGCCCGCACCGGCTTCGGCCTGGATATGTTCGACCTGCa                                                        >  2:382010/1‑150 (MQ=255)
                                                                                                 tCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAGGCTCATGCCCGCACCGGCTTCGGCCTGGa                                                                       <  2:4303/132‑1 (MQ=255)
                                                                                                 tCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAGGCTCATGCCCGCACCGGCTTCGGCCTGGa                                                                       >  1:4303/1‑132 (MQ=255)
                                                                                                 tCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAGGCTCATGCCCGCACCGGCTTCGGCCTGGATATGTTCGACCTGCAGcc                                                     >  1:463675/1‑150 (MQ=255)
                                                                                                     gTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAGGCTCATGCCCGCACCGGCTTCGGCCTGGATATGTTCGACCTGCAGCCGACa                                                 <  2:370679/150‑1 (MQ=255)
                                                                                                           aTAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAGGCTCATGCCCGCACCGGCTTCGGCCTGGATATGTTCGACCTGCAGCCGACAGAGTTc                                           <  1:648917/150‑1 (MQ=255)
                                                                                                                ccTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAGGc                                                                                                  >  2:234213/1‑90 (MQ=255)
                                                                                                                ccTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAGGc                                                                                                  <  1:234213/90‑1 (MQ=255)
                                                                                                                         ggCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAGGCTCATGCCCGCACCGGCTTCGGCCTGGATATGTTCGACCTGCAGCCGACAGAGTTCGTCACTACGAAAAc                             <  1:9977/150‑1 (MQ=255)
                                                                                                                          gCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAGGCTCATGCCCGCa                                                                                       <  2:1086762/91‑1 (MQ=255)
                                                                                                                          gCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAGGCTCATGCCCGCa                                                                                       >  1:1086762/1‑91 (MQ=255)
                                                                                                                            agagGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAGGCTCATGCCCGCACCGGCTTCGGCCTGGATATGTTCGACCTGCAGCCGACAGAGTTCGTCACTACGAAAAcccc                          <  2:463675/150‑1 (MQ=255)
                                                                                                                                gCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAGGCTCATGCCCGCACCGGCTTCGGCCTGGATATGTTCGACCTGCAGCCGACAGAGTTCGTCACTACGAAAACCCCGCCa                      <  2:747897/150‑1 (MQ=255)
                                                                                                                                gCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAGGCTCATGCCCGCACCGGCTTCGGCCTGGATATGTTCGACCTGCAGCCGACAGAGTTCGTCACTACGAAAACCCCGCCa                      <  2:133422/150‑1 (MQ=255)
                                                                                                                                gCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAGGCTCATGCCCGCACCGGCTTCGGCCTGGATATGTTCGACCTGCAGCCGACAGAGTTCGTCACTACGAAAACCCCGCCa                      <  2:133657/150‑1 (MQ=255)
                                                                                                                                       aaGGCTGGTACGTAACGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAGGCTCATGCCCGCACCGGCTTCGGCCTGGATATGTTCGACCTGCAGCCGACAGAGTTCGTCACTACGAAAACCCCGCCAGAAGCCg               <  1:752463/150‑1 (MQ=255)
                                                                                                                                                  gTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAGGCTCATGCCCGCACCGGCTTCGGCCTGGATATGTTCGACCTGCAGCCGACAGAGTTCGTCACTACGAAAACCCCGCCAGAAGCCGATCAGCAGCaa    >  2:973784/1‑150 (MQ=255)
                                                                                                                                                  gTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAGGCTCATGCCCGCACCGGCTTCGGCCTGGATATGTTCGACCTGCAGCCGACAGAGTTCGTCACTACGAAAACCCCGCCAGAAGCCGATCAGCAGCaa    >  2:380165/1‑150 (MQ=255)
                                                                                                                                                   tAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAGGCTCATGCCCGCACCGGCTTCGGCCTGGATATGTTCGACCTGCAGCCGACAGAGTTCGTCACTACGAAAACCCCGCCAGAAGCCGATCAGCAGCAAc   <  1:319863/150‑1 (MQ=255)
                                                                                                                                                    aaTGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAGGCTCATGCCCGCACCGGCTTCGGCCTGGATATGTTCGACCTGCAGCCGACAGAGTTCGTCACTACGAAAACCCCGCCAGAAGCCGATCAGCAGCAAcg  >  1:183178/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
TGATCAGGCTGCCGGGATGCAAGGCCTCGTCGCTCAGGTGCCCCCAGCGATCCAGGCGACAGAACACTGCACCACGGGCAATGCCGGCGGCACTGATCCAGTCGAGATAAGCCTGCACCGGGCAGAGGCACTTCAAGGCTGGTACGTAATGAGTCGTGCCGAGGTTGTCCCTGTCACCCTTACTTTGTGGCAGGAACAGGCTCATGCCCGCACCGGCTTCGGCCTGGATATGTTCGACCTGCAGCCGACAGAGTTCGTCACTACGAAAACCCCGCCAGAAGCCGATCAGCAGCAACG  >  NC_002947/3379477‑3379773

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