Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,716,692 G→A 8.3% K469K (AAG→AAA)  cti → esterified fatty acid cis/trans isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,716,6920GA 8.3% 133.2 / 2.3 48K469K (AAG→AAA) ctiesterified fatty acid cis/trans isomerase
Reads supporting (aligned to +/- strand):  ref base G (22/22);  new base A (2/2);  total (24/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.88e-01

AAGCGCAACGACTACGAGAAACTGCGCCGTGAAGCCTATGCCGAAATGCCGGCACCGAGATGGTCGACGCTGTGGGCCGGTAACGACAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTACCAGCTGGCGGTCAACTTCGATGTGTATGGCAACGTTTCGCACCAGTTGCAGACGCGCCTGTATTTCGACCTGATCC  >  NC_002947/2716555‑2716837
                                                                                                                                         |                                                                                                                                                 
aaGCGCAACGACTACGAGAAACTGCGCCGTGAAGCCTATGCCGAAATGCCGGCACCGAGATGGTCGACGCTGTGGGCCGGTAACGACAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTggg                                                                                                                                        <  8:115395/149‑1 (MQ=255)
    gcAACGACTACGAGAAACTGCGCCGTGAAGGCTATCCCGAAATGCCGGCACCGAGATGGGCGACGCTGTGGGCCGGTAACGACAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGAt                                                                                                                                    <  6:172025/149‑1 (MQ=255)
            tACGAGAAACTGCGCCGTGAAGCCTATGCCGAAATGCCGGCACCGAGATGGTCGACGCTGTGGGCCGGTAACGACAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCt                                                                                                                            >  4:36618/1‑149 (MQ=255)
            tACGAGAAACTGCGCCGTGAAGCCTATGCCGAAATGCCGGCACCGAGATGGTCGACGCTGTGGGCCGGTAACGACAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCt                                                                                                                            >  3:322388/1‑149 (MQ=255)
               gagaAACTGCGCCGTGAAGCCTATGCCGAAATGCCGGCACCGAGATGGTCGACGCTGTGGGCCGGTAACGACAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGGGACCAAGGGCCTGGTGGGGGATGTGCCGCTGAc                                                                                                                         >  7:228386/1‑149 (MQ=255)
                  aaaCTGCGCCGTGAAGCCTATGCCGAAATGCCGGCACCGAGATGGTCGACGCTGTGGGCCGGTAACGACAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCgt                                                                                                                      >  3:281090/1‑149 (MQ=255)
                            gTGAAGCCTATGCCGAAATGCCGTCACCGAGATGGTCGACGCTGTGGGCCGGTAACGACAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTc                                                                                                            <  6:256651/149‑1 (MQ=255)
                               aaGCCTATGCCGAAATGCCGGCACCGAGATGGTCGACGCTGTGGGCCGGTAACGACAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGa                                                                                                          <  2:129929/148‑1 (MQ=255)
                               aaGCCTATGCCGAAATGCCGGCACCGAGATGGTCGACGCTGTGGGCCGGTAACGACAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAAGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGAc                                                                                                         >  8:198148/1‑149 (MQ=255)
                               aaGCCTATGCCGAAATGCCGGCACCGAGATGGTCGACGCTGTGGGCCGGTAACGACAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAAGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGAc                                                                                                         >  4:312290/1‑149 (MQ=255)
                                aGCCTATGCCGAAATGCCGGCACCGAGATGGTCGACGCTGTGGGCCGGTAACGACAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGAc                                                                                                         >  8:39271/1‑148 (MQ=255)
                                 gCCTATGCCGAAATGCCGGCACCGAGATGGTCGACGCTGTGGGCCGGTAACGACAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCATGGGCTTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTa                                                                                                       >  7:220560/1‑149 (MQ=255)
                                   cTATGCCGAAATGCCGGCACCGAGATGGTCGACGCTGTGGGCCGGTAACGACAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTAc                                                                                                      >  5:298277/1‑148 (MQ=255)
                                         cGAAATGCCGGCACCGAGATGGTCGACGCTGTGGGCCGGTAACGACAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCgttgt                                                                                               >  6:70583/1‑149 (MQ=255)
                                           aaaTGCCGGCACCGAGATGGTCGACGCTGTGGGCCGGTAACGACAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTc                                                                                             >  1:331394/1‑149 (MQ=255)
                                             aTGCCGGCACCGAGATGGTCGACGGTGTGGGCCGGTAACGACAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGa                                                                                           <  7:305496/149‑1 (MQ=255)
                                                   gCACCGAGATGGTCGACGCTGTGGGCCGGTAACGACAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCAc                                                                                     <  3:94417/149‑1 (MQ=255)
                                                     aCCGAGATGGTCGACGCTGTGGGCCGGTAACGACAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGt                                                                                   <  1:173956/149‑1 (MQ=255)
                                                         agaTGGTCGACGCTGTGGGCCGGTAACGACAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCgtg                                                                                                                     >  1:95871/1‑111 (MQ=255)
                                                         agaTGGTCGACGCTGTGGGCCGGTAACGACAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCgtg                                                                                                                     <  2:95871/111‑1 (MQ=255)
                                                         agaTGGTCGACGCTGTGGGCCGGTAACGACAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTa                                                                               <  3:36618/149‑1 (MQ=255)
                                                         agaTGGTCGACGCTGTGGGCCGGTAACGACAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGATCGAGCGCACGTATTa                                                                               <  5:70583/149‑1 (MQ=255)
                                                              gTCGACGCTGTGGGCCGGTAACGACAATGCGCTGCTTAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTAcc                                                                             <  8:261758/146‑1 (MQ=255)
                                                              gTCGACGCTGTGGGCCGGTAACGACAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTAcc                                                                             >  7:261758/1‑146 (MQ=255)
                                                                  aCGCTGTGGGCCGGTAACGACAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTACCAGCTGGc                                                                      <  4:281090/149‑1 (MQ=255)
                                                                   cGCTGTGGGCCGGTAACGACAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTACCAGCTGGCg                                                                     <  6:3005/149‑1 (MQ=255)
                                                                        tgGGCCGGTAACGACAATGCGCTGCTGAGCATCTTCCGTGACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTACCAGCTGGCGGTCaa                                                                <  5:323965/149‑1 (MQ=255)
                                                                          ggCCGGTAACGGCAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTACCAGCTGGCGGTCAACt                                                              >  8:272686/1‑149 (MQ=255)
                                                                          ggCCGGTAACGACAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTACCAGCTGGCGGTCAACt                                                              >  4:279040/1‑149 (MQ=255)
                                                                              ggTAACGACAATGCGCTGCTGAGCATCTTGCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTACCAGCTGGCGGTCAACTTCGa                                                          <  4:322388/149‑1 (MQ=255)
                                                                                    gACAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTAc                                                                              >  8:267659/1‑123 (MQ=255)
                                                                                    gACAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTAc                                                                              <  7:267659/123‑1 (MQ=255)
                                                                                     aCAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAAGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTACCAGCTGGCGGTCAACTTCGATGTGTAt                                                   <  7:198148/149‑1 (MQ=255)
                                                                                     aCAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAAGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTACCAGCTGGCGGTCAACTTCGATGTGTAt                                                   <  3:312290/149‑1 (MQ=255)
                                                                                          gcgcTGCTGAGCATCTTCCGTCACTTCGACAGCGAATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTACCAGCTGGCGGTCAACTTCGATGTGTATGGCaa                                              >  1:223020/1‑149 (MQ=255)
                                                                                             ctgctgAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTACCAGCTGGCGGTCAACTTCGATGTGTATGGCAACGt                                           <  7:272686/149‑1 (MQ=255)
                                                                                                 tgAGGATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTACCAGCTGGCGGTCAACTTCGa                                                          <  6:44908/130‑1 (MQ=255)
                                                                                                 tgAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTACCAGCTGGCGGTCAACTTCGa                                                          >  5:44908/1‑130 (MQ=255)
                                                                                                      aTCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTACCAGCTGGCGGTCAACTTCGATGTGTATGGCAACGTTTCGCACCa                                  <  2:331394/149‑1 (MQ=255)
                                                                                                        cTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTACCAGCTGGCGGTCAACTTCGATGTGTATGGCAACGTTTCGCACCAGt                                >  7:270976/1‑149 (MQ=255)
                                                                                                              tCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTACCAGCTGGCGGTCAACTTCGATGTGTATGGCAACGTTTCGCACCAGTTGCAGa                          >  7:305297/1‑149 (MQ=255)
                                                                                                                 cTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTACCAGCTGGCGGTCAACTTCGATGTGTATGGCAACGTTTCGCACCAGTTGCAGAcgc                       <  7:39271/149‑1 (MQ=255)
                                                                                                                 cTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTACCAGCTGGCGGTCAACTTCGATGTGTATGGCAACGTTTCGCACCAGTTGCAGAcgc                       <  8:270976/149‑1 (MQ=255)
                                                                                                                                ggTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTACCAGCTGGCGGTCAACTTCGATGTGTATGGCAACGTTTCGCACCAGTTGCAGACGCGCCTGTATTTCGAc         <  8:305297/148‑1 (MQ=255)
                                                                                                                                     ccAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTACCAGCTGGCGGTCAACTTCGATGTGTATGGCAACGTTTCGCACCAGTTGCAGACGCGCCTGTATTTCGACCTGATc   >  6:264285/1‑149 (MQ=255)
                                                                                                                                      cAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTACCAGCTGGCGGTCAACTTCGATGTGTATGGCAACGTTTCGCACCAGTTGCAGACGCGCCTGTATTTCGACCTGATcc  >  6:71138/1‑149 (MQ=255)
                                                                                                                                      cAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTACCAGCTGGCGGTCAACTTCGATGTGTATGGCAACGTTTCGCACCAGTTGCAGACGCGCCTGTATTTCGACCTGATcc  >  4:41666/1‑149 (MQ=255)
                                                                                                                                       aaGGGCCTGGTGGGGGATGTGCCGCTGCCCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTACCAGCTGGAGGTCAACTTCGATGTGTATGGCAACGTTTCGCACCAGTTGCAGACGCGCCTGTATTTCTACCTGATcc  >  3:206970/1‑148 (MQ=255)
                                                                                                                                         gggCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTACCAGCTGGCGGTCAACTTCGATGTGTATGGCAAc                                             <  6:87009/103‑1 (MQ=255)
                                                                                                                                         gggCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTACCAGCTGGCGGTCAACTTCGATGTGTATGGCAAc                                             >  5:87009/1‑103 (MQ=255)
                                                                                                                                         |                                                                                                                                                 
AAGCGCAACGACTACGAGAAACTGCGCCGTGAAGCCTATGCCGAAATGCCGGCACCGAGATGGTCGACGCTGTGGGCCGGTAACGACAATGCGCTGCTGAGCATCTTCCGTCACTTCGACAGCGCATCGGTGACCAAGGGCCTGGTGGGGGATGTGCCGCTGACCGTGTGGCTGTTCGACTACCCGTTGTTCGAGCGCACGTATTACCAGCTGGCGGTCAACTTCGATGTGTATGGCAACGTTTCGCACCAGTTGCAGACGCGCCTGTATTTCGACCTGATCC  >  NC_002947/2716555‑2716837

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