Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,376,949 A→T 7.5% noncoding (111/117 nt)
Q23H (CAA→CAT) 		 PP_mr38 →
tnpB →		 ALIL regulatory region
transposase subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,376,9490AT 7.5% 153.3 / 3.0 53noncoding (111/117 nt)
Q23H (CAA→CAT) 		PP_mr38
tnpB		ALIL regulatory region
transposase subunit B
Reads supporting (aligned to +/- strand):  ref base A (26/23);  new base T (2/2);  total (28/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 = 8.21e-01

AAGCGAGCCAACGAAATCCTCAAGGTCGCCAGCGCTTTTTCCGCCCAGGCGGAGCTCGACCGCCGTTTGAAGTGATCTGGCAGCTGATCGACAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGAGGAACTTGCTGTACAAATTCACCGTGTCTGGCAGGAAAACTACGAGGTCTACGGGGCTAGGAAAGTC  >  NC_002947/3376810‑3377096
                                                                                                                                           |                                                                                                                                                   
aaGCGAGCCAACGAAATCCTCAAGGTCGCCAGCGCTTTTTCCGCCCAGGCGGAGCTCGACCGCCGTTTGAAGTGATCTGGCAGCTGATCGACAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGtctc                                                                                                                                              <  2:213463/147‑1 (MQ=255)
     aGCCAACGAAATCCTCAAGGTCGCCAGCGCTTTTTCCGCCCAGGCGGAGCTCGACCGCCGTTTGAAGTGATCTGGCAGCTGATCGACAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGc                                                                                                                                       <  1:170815/149‑1 (MQ=255)
         aaCGAAATCCTCAAGGTCGCCAGCGCTTTTTCCGCCCAGGCGGAGCTCGACCGCCGTTTGAAGTGATCTGGCAGCTGATCGACAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCa                                                                                                                                            <  4:18641/140‑1 (MQ=255)
         aaCGAAATCCTCAAGGTCGCCAGCGCTTTTTCCGCCCAGGCGGAGCTCGACCGCCGTTTGAAGTGATCTGGCAGCTGATCGACAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCa                                                                                                                                            >  3:18641/1‑140 (MQ=255)
         aaCGAAATCCTCAAGGTCGCCAGCGCTTTTTCCGCCCAGGCGGAGCTCGACCGCCGTTTGAAGTGATCTGGCAGCTGATCGACAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTAc                                                                                                                                   <  2:50634/149‑1 (MQ=255)
              aaTCCTCAAGGTCTCCAGCGCTTTTTCCGCCCAGGCGGAGCTCGACCGCCGTTTGAAGTGATCTGGCAGCTGATCGACAAGCACCGGCATAGCTTTGGAGTCGAGCCACTGTTCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGgcg                                                                                                                              <  1:221893/149‑1 (MQ=255)
                   tCAAGGTCGCCAGCGCTTTTTCCGCCCAGGCGGAGCTCGACCGCCGTTTGAAGTGATCTGGCAGCTGATCGACAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCATGTGTCTCCATCGGCTTACCGGCGCCAtg                                                                                                                         >  8:342872/1‑149 (MQ=255)
                             cAGCGCTTTTTCCGCCCAGGCGGAGCTCGACCGCCGTTTGAAGTGATCTGGCAGCTGATCGACAAGCACTGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCt                                                                                                               <  5:321385/149‑1 (MQ=255)
                               gcgcTTTTTCCGCCCAGGCGGAGCTCGACCGCCGTTTGAAGTGATCTGGCAGCTGATCGACAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCaa                                                                                                                                                     <  1:134106/109‑1 (MQ=255)
                               gcgcTTTTTCCGCCCAGGCGGAGCTCGACCGCCGTTTGAAGTGATCTGGCAGCTGATCGACAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCaa                                                                                                                                                     >  2:134106/1‑109 (MQ=255)
                                  cTTTTTCCGCCCATGCGGAGCTCGACCGCCGTTTGAAGTGATCTGGCAGCTGATCGACAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGt                                                                                                           >  5:81178/1‑148 (MQ=255)
                                  cTTTTTCCGCCCAGGCGGAGCTCGACCGCCGTTTGAAGTGATCTGGCAGCTGATCGACAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGgcg                                                                                                                              >  8:217599/1‑129 (MQ=255)
                                  cTTTTTCCGCCCAGGCGGAGCTCGACCGCCGTTTGAAGTGATCTGGCAGCTGATCGACAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGgcg                                                                                                                              <  7:217599/129‑1 (MQ=255)
                                   tttttCCGCCCAGGCGGAGCTCGACCGCCGTTTGAAGTGATCTGGCAGCTGATCGACAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGa                                                                                                         >  3:15494/1‑149 (MQ=255)
                                    ttttCCGCCCAGGCGGAGCTCGACCGCCGTTTGAAGTGATCTGGCAGCTGATCGACAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGAc                                                                                                        >  5:154527/1‑149 (MQ=255)
                                           cccAGGCGGAGCTCGACCGCCGTTTGAAGTGATCTGGCAGCTGATCGACAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCac                                                                                                 >  5:365971/1‑149 (MQ=255)
                                                gCGGAGCTCGACCGCCGTTTGAAGTGATCTGGCAGCTGATCGACAAGCACCGGCAAAGCTTTGGAGTCGACCCACTGTGCCGCGTGTAGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGAACGGCTACGTGACCGCTCaagactc                                                                                            >  5:304918/1‑147 (MQ=255)
                                                         gACCGCCGTTTGAAGTGATCTGGCAGCTGATCGACAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGa                                                                                   <  7:276298/149‑1 (MQ=255)
                                                            cgccgTTTGAAGTGATCTGGCAGCTGATCGACAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAg                                                                                  <  7:144716/147‑1 (MQ=255)
                                                            cgccgTTTGAAGTGATCTGGCAGCTGATCGACAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAg                                                                                  >  8:144716/1‑147 (MQ=255)
                                                              ccgTTTGAAGTGATCTGGCAGCTGATCGACAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGa                                                                                                         <  3:403993/122‑1 (MQ=255)
                                                              ccgTTTGAAGTGATCTGGCAGCTGATCGACAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGa                                                                                                         >  4:403993/1‑122 (MQ=255)
                                                                    gAAGTGATCTGGCAGCTGATCGACAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTg                                                                         <  6:154527/148‑1 (MQ=255)
                                                                    gAAGTGATCTGGCAGCTGATCGACAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGa                                                                        >  2:76117/1‑149 (MQ=255)
                                                                       gTGATCTGGCAGCTGATCGACAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGa                                                                     <  3:230455/149‑1 (MQ=255)
                                                                        tGATCTGGCAGCTGATCGACAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCATGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGAg                                                                    >  4:63768/1‑149 (MQ=255)
                                                                           tCTGGCAGCTGATCGACAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGAGGaa                                                                 >  4:362158/1‑149 (MQ=255)
                                                                               gCAGCTGATCGACAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGAGGAACTTg                                                             >  6:161839/1‑149 (MQ=255)
                                                                                 aGCTGATCGACAAGCACCGGCATAGCTTTGGAGTCGATCCACTTTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGTCGCCATGCTGCACGGCTACGTGAGCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGAGGAACTTGCt                                                           <  6:81178/149‑1 (MQ=255)
                                                                                      aTCGACAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGAGGAACTTGCTGTa                                                        <  6:365971/147‑1 (MQ=255)
                                                                                       tCGACAAGCACCGGCAAAGCTTTTGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGa                                                                        >  7:68680/1‑130 (MQ=255)
                                                                                       tCGACAAGCACCGGCAAAGCTTTTGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGa                                                                        <  8:68680/130‑1 (MQ=255)
                                                                                       tCGACAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGAGGAACTTGCTGTACaa                                                     <  1:180124/149‑1 (MQ=255)
                                                                                         gACAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGAGGAACTTGCTGTACAAAt                                                   >  6:416460/1‑149 (MQ=255)
                                                                                          aCAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGAGGAACTTGCTGTACAAAtt                                                  >  8:363325/1‑149 (MQ=255)
                                                                                              gCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGAGGAACTTGCTGTACAAATTCAcc                                              <  5:161839/149‑1 (MQ=255)
                                                                                              gCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGAGGAACTTGCTGTACAAATTCAcc                                              >  8:16920/1‑149 (MQ=255)
                                                                                                aCCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGAGGAACTTGCTGTACAAATTCACCgt                                            <  7:363325/149‑1 (MQ=255)
                                                                                                     cAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGAGGAACTTGCTGTACAAATTCACCGTGTCTg                                       >  6:180403/1‑149 (MQ=255)
                                                                                                      aaaGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCATGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGAGGAACTTGCTGTACAAATTCACCGTGTCTgg                                      <  3:63768/149‑1 (MQ=255)
                                                                                                         gCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGAGGAACTTGCTGTACAAATTCACCGTGTCTGGCAg                                   <  1:76117/149‑1 (MQ=255)
                                                                                                            ttGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGAGGAACTTGCTGTACAAATTCACCGTGTCTGGCAGGa                                 >  2:145347/1‑148 (MQ=255)
                                                                                                                 gTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGAGGAACTTGCTGTACAAATTCACCGTGTCTGGCAGGAAAACTa                           >  6:105431/1‑149 (MQ=255)
                                                                                                                    gAGCCACTGTGCCGCGTGTTGCATGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGAGGAACTTGCTGTACAAATTCACCGTGTCTGGCAGGAAAACTACGa                        <  7:342872/149‑1 (MQ=255)
                                                                                                                    gAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGAGGAACTTGCTGTACAAATTCACCGTGTCTGGCAGGAAAACTAc                          >  3:120883/1‑147 (MQ=255)
                                                                                                                    gAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGAGGAACTTGCTGTACAAATTCACCGTGTCTGGCAGGAAAACTACGa                        <  5:105431/149‑1 (MQ=255)
                                                                                                                       ccACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGAGGAACTTGCTGTACAAATTCACCGTGTCTGGCAGGAAAACTACGAGGt                     <  6:346930/149‑1 (MQ=255)
                                                                                                                         aCTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGa                                                                     >  5:105027/1‑99 (MQ=255)
                                                                                                                         aCTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGa                                                                     <  6:105027/99‑1 (MQ=255)
                                                                                                                         aCTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGAGGAACTTGCTGTACAAATTCa                                                <  2:230634/120‑1 (MQ=255)
                                                                                                                         aCTGTGCCGCGTGTTGCAAGTGTCTCCAACGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGAGGAACTTGCTGTACAAATTCa                                                >  1:230634/1‑120 (MQ=255)
                                                                                                                          cTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGAGGAACTTGCTGTACAAATTCACCGTGTCTGGCAGGAAAACTAC‑AGGTTTAc                 >  5:128864/1‑149 (MQ=255)
                                                                                                                              gCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGAGGAACTTGCTGTACAAATTCACCGTGTCTGGCAGGAAAACTACGAGGTCTACggg              >  7:333417/1‑149 (MQ=255)
                                                                                                                                     gtTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGAGGAACTTGCTGTACAAATTCACCGTGTCTGGCAGGAAAACTACGAGGTCTACGGGGCTAgg        >  4:67117/1‑148 (MQ=255)
                                                                                                                                         cAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGAGGAACTTGCTGTACAAATTCACCGTGTCTGGCAGGAAAACTACGAGGTCTACGGGGCTAGGAAAGt   >  5:222704/1‑149 (MQ=255)
                                                                                                                                          aaGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGAGGAACTTGCTGTACAAATTCACCGTGTCTGGCAGGAAAACTACGAGGTCTACGGGGCTAGGAAAGTc  >  7:351747/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
AAGCGAGCCAACGAAATCCTCAAGGTCGCCAGCGCTTTTTCCGCCCAGGCGGAGCTCGACCGCCGTTTGAAGTGATCTGGCAGCTGATCGACAAGCACCGGCAAAGCTTTGGAGTCGAGCCACTGTGCCGCGTGTTGCAAGTGTCTCCATCGGCTTACCGGCGCCATGCTGCACGGCTACGTGACCGCTCACGACGCTCCGCCCGAGCGATCCGTGACGAGGAACTTGCTGTACAAATTCACCGTGTCTGGCAGGAAAACTACGAGGTCTACGGGGCTAGGAAAGTC  >  NC_002947/3376810‑3377096

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: ‑