Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,474,270 C→A 11.4% L422L (CTG→CTT) ‡ PP_3964 ← transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,474,2700CA11.4% 153.3 / 5.5 50L422L (CTG→CTT) ‡PP_3964transposase
Reads supporting (aligned to +/- strand):  ref base C (18/26);  new base A (3/3);  total (21/29)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACTCATGATTGCCGCAGCCCGCTTACCGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTC  >  NC_002947/4474127‑4474416
                                                                                                                                               |                                                                                                                                                  
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                                                                                                  <  8:158664/149‑1 (MQ=255)
                                             gCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGa                                                                                                  >  5:196930/1‑149 (MQ=255)
                                              acaacaaccaggtcgaaaaccagatacggccatgggcactcgggcgctcgaactggttGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTg                                                                                                       >  3:14783‑M2/59‑143 (MQ=255)
                                              acaacaaccaggtcgaaaaccagatacggccatgggcactcgggcgctcgaactggttGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTg                                                                                                       <  4:14783‑M2/85‑1 (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)
                                                                       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                          >  3:49560/1‑149 (MQ=255)
                                                                                                                     acaGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTc                          >  7:161987/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                                             >  3:63578/1‑109 (MQ=255)
                                                                                                                                          cgcgTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGa                                             <  4:63578/109‑1 (MQ=255)
                                                                                                                                             gTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCtt   >  5:158004/1‑148 (MQ=255)
                                                                                                                                             gTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTtc  <  3:17623/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
AAACTCATGATTGCCGCAGCCCGCTTACCGCTGCGCAGGGATCCCGCAAATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTC  >  NC_002947/4474127‑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: ‑