Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,474,268 G→T 11.7% T423K (ACG→AAG) ‡ PP_3964 ← transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,474,2680GT11.7% 138.7 / 7.3 49T423K (ACG→AAG) ‡PP_3964transposase
Reads supporting (aligned to +/- strand):  ref base G (18/25);  new base T (3/3);  total (21/28)
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.24e-01

GATCAAACTCATGATTGCCGCAGCCCGCTTACCGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTC  >  NC_002947/4474123‑4474416
                                                                                                                                                 |                                                                                                                                                    
gATCAAACTCATGATTGCCGCAGCCCGCTTACCGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGt                                                                                                                                                     >  3:69263/1‑147 (MQ=255)
    aaaCTCATGATTGCCGCAGCCCGCTTACCGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAgcgc                                                                                                                                               <  4:250464/149‑1 (MQ=255)
        tCATGATTGCCGCAGCCCGCTTACCGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCTGACGCGCTAcc                                                                                                                                           >  3:71631/1‑149 (MQ=255)
        tCATGATTGCCGCAGCCCGCTTACCGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTAcc                                                                                                                                           >  8:34816/1‑149 (MQ=255)
             aTTGCCGCAGCCCGCTTACCGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCTGACGCGCTAcctgga                                                                                                                                       <  4:71631/148‑5 (MQ=18)
             aTTGCCGCAGCCCGCTTACCGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCg                                                                                                                                      >  2:163818/1‑149 (MQ=255)
             aTTGCCGCAGCCCGCTTACCGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCg                                                                                                                                      <  5:33590/149‑1 (MQ=255)
                  cGCAGCCCGCTTACCGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCa                                                                                                                                 <  5:212193/149‑1 (MQ=255)
                   gCAGCCCGCTTACCGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGTCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAg                                                                                                                                <  7:192412/149‑1 (MQ=255)
                        ccGCTTACCGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCGGGTAGCGCTGACGCGCTAcctgga                                                                                                                                       <  8:116761/137‑5 (MQ=18)
                        ccGCTTACCGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCGGGTAGCGCTGACGCGCTAcctgga                                                                                                                                       >  7:116761/1‑133 (MQ=18)
                          gCTTACCGCTGCGCGGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTaa                                                                                                                         <  3:146850/149‑1 (MQ=255)
                              aCCGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGAGCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATc                                                                                                                       <  2:171489/147‑1 (MQ=255)
                               ccGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGa                                                                                                                     <  1:163818/148‑1 (MQ=255)
                               ccGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAg                                                                                                                    <  6:165154/149‑1 (MQ=255)
                                   tGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTAGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCt                                                                                                                >  6:268906/1‑149 (MQ=255)
                                           gATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGt                                                                                                        <  1:34876/149‑1 (MQ=255)
                                               ccGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCGGGTAGCGCTGACGCGCTAcctgga                                                                                                                                       <  7:228526/114‑5 (MQ=14)
                                               ccGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCGGGTAGCGCTGACGCGCTAcctgga                                                                                                                                       >  8:228526/1‑110 (MQ=14)
                                                 gCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGa                                                                                                  >  5:196930/1‑149 (MQ=255)
                                                 gCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGa                                                                                                  <  8:158664/149‑1 (MQ=255)
                                                  acaacaaccaggtcgaaaaccagatacggccatgggcactcgggcgctcgaactggttGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTg                                                                                                       <  4:14783‑M2/85‑1 (MQ=255)
                                                  acaacaaccaggtcgaaaaccagatacggccatgggcactcgggcgctcgaactggttGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTg                                                                                                       >  3:14783‑M2/59‑143 (MQ=255)
                                                   aaaTAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTCGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATc                                                                                                <  5:268906/149‑1 (MQ=255)
                                                     aTAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCa                                                                                              >  2:100014/1‑149 (MQ=255)
                                                      caaccaggtcgaaaaccagatacggccatgggcactcgggcgctcgaactggttGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCAt                                                                                             >  4:17623‑M2/55‑149 (MQ=255)
                                                        ACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCAtt                                                                                            >  4:69263‑M1/1‑148 (MQ=255)
                                                         CCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGt                                                                                                                                                     >  4:236404‑M1/1‑90 (MQ=255)
                                                         CCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGt                                                                                                                                                     <  3:236404‑M1/90‑1 (MQ=255)
                                                                           tGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAg                                                                        <  7:34816/149‑1 (MQ=255)
                                                                              ccATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGcc                                                                      >  7:71768/1‑148 (MQ=255)
                                                                                aTGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAAc                                                                   <  1:100014/149‑1 (MQ=255)
                                                                                     cGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATc                                                                >  3:20760/1‑147 (MQ=255)
                                                                                             gTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATg                                                        <  4:20760/147‑1 (MQ=255)
                                                                                                    aCCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGtttttt                                               <  8:71768/149‑1 (MQ=255)
                                                                                                      ctggttGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGa                                             <  6:196930‑M2/143‑1 (MQ=255)
                                                                                                      ctggttGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGa                                             <  1:249590‑M2/143‑1 (MQ=255)
                                                                                                           tGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATc                                        >  5:6270‑M2/2‑149 (MQ=255)
                                                                                                             ttgtCGATGGGCACAGGCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCg                                       <  7:229544/148‑1 (MQ=255)
                                                                                                                tCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGtttttt                                               <  5:73290/137‑1 (MQ=255)
                                                                                                                tCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGtttttt                                               >  6:73290/1‑137 (MQ=255)
                                                                                                                  gATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCAcc                                  <  8:82684/148‑1 (MQ=255)
                                                                                                                  gATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCg                                 >  3:209745/1‑149 (MQ=255)
                                                                                                                    tGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGc                                <  4:209745/148‑1 (MQ=255)
                                                                                                                         acaGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTc                          >  7:161987/1‑149 (MQ=255)
                                                                                                                         acaGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTc                          >  3:49560/1‑149 (MQ=255)
                                                                                                                              cccATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTc                     <  4:49560/149‑1 (MQ=255)
                                                                                                                                       cAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTc                                                           >  7:225345/1‑102 (MQ=255)
                                                                                                                                       cAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTc                                                           <  8:225345/102‑1 (MQ=255)
                                                                                                                                              cgcgTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGa                                             <  4:63578/109‑1 (MQ=255)
                                                                                                                                              cgcgTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGa                                             >  3:63578/1‑109 (MQ=255)
                                                                                                                                                 gTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCtt   >  5:158004/1‑148 (MQ=255)
                                                                                                                                                 gTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTtc  <  3:17623/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GATCAAACTCATGATTGCCGCAGCCCGCTTACCGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTC  >  NC_002947/4474123‑4474416

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: