Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,099,496:1 +A 11.4% coding (232/1470 nt) PP_2709 → long‑chain fatty acid‑‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,099,4961.A11.4% 114.7 / 11.8 35coding (232/1470 nt)PP_2709long‑chain fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base . (13/18);  new base A (2/2);  total (15/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.24e-01

GCGCAGCGTTACACCTACCGGCAACTGCTCAACGAAGTCGAATCCCGTGCGAGCTGCTTGCGCAGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCGTGTGTGAGGCGCAGTGGGCCGAACAGCTCATTGCGCTTGGCTTTGCCCGGCAGGCGCCAGCTGGC  >  NC_002947/3099349‑3099642
                                                                                                                                                    |                                                                                                                                                  
gcgcAGCGTTACACCTACCGGCAACTGCTCAACGAAGTCGAATCCCGTGCGAGCTGCTTGCGCAGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAAcgc                                                                                                                                                  <  2:68686/150‑1 (MQ=255)
  gcAGCGTTACACCTACCGGCAACTGCTCAACGAAGTCGAATCCCGTGCGAGCTGCTTGCGCAGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACgcgc                                                                                                                                                >  1:1131217/1‑150 (MQ=255)
            acCTACCGGCAACTGCTCAACGAAGTCGAATCCCGTGCGAGCTGCTTGCGCAGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACa                                                                                                                                      <  1:536821/150‑1 (MQ=255)
              cTACCGGCAACTGCTCAACGAAGTCGAATCCCGTGCGAGCTGCTTGCGCAGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGt                                                                                                                                    >  2:240340/1‑150 (MQ=255)
                  cGGCAACTGCTCAACGAAGTCGAATCCCGTGCGAGCTGCTTGCGCAGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTgccg                                                                                                                                >  2:65610/1‑150 (MQ=255)
                              aaCGAAGTCGAATCCCGTGCGAGCTGCTTGCGCAGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAgc                                                                                                                    >  2:45799/1‑150 (MQ=255)
                              aaCGAAGTCGAATCCCGTGCGAGCTGCTTGCGCAGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCGTTCGGCGAACAGCGccg                                                                                                                                              <  1:1065350/124‑2 (MQ=255)
                              aaCGAAGTCGAATCCCGTGCGAGCTGCTTGCGCAGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCGTTCGGCGAACAGCGccg                                                                                                                                              >  2:1065350/1‑123 (MQ=255)
                                  aaGTCGAATCCCGTGCGAGCTGCTTGCGCAGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGa                                                                                                                >  1:1068090/1‑150 (MQ=255)
                                     tCGAATCCCGTGCGAGCTGCTTGCGCAGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCGTTCGGCGAACAGCGccg                                                                                                                                              >  1:588198/1‑116 (MQ=255)
                                     tCGAATCCCGTGCGAGCTGCTTGCGCAGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCGTTCGGCGAACAGCGccg                                                                                                                                              <  2:588198/117‑2 (MQ=255)
                                                        cgcaagcAGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTgca                                                                                          <  2:1200984/145‑1 (MQ=255)
                                                             gcAGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCaaa                                                                                     <  1:240340/150‑1 (MQ=255)
                                                               aGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGc                                                                                   <  2:862566/150‑1 (MQ=255)
                                                                gCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACTCAGTTCCAGCACTGCTTGCAGCAAAGCg                                                                                  <  1:65610/150‑1 (MQ=255)
                                                                    accaccGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTCTTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAgg                                                                              >  1:125146/1‑149 (MQ=255)
                                                                      caccGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAg                                                                               <  1:471292/147‑1 (MQ=255)
                                                                      caccGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAg                                                                               >  2:471292/1‑147 (MQ=255)
                                                                        ccGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGAcg                                                                          <  2:934743/150‑1 (MQ=255)
                                                                         cGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGAcgg                                                                         >  1:127276/1‑150 (MQ=255)
                                                                           ggAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGAcggcg                                                                       <  2:269865/150‑1 (MQ=255)
                                                                               ccTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACTCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCg                                                                   <  2:1068090/150‑1 (MQ=255)
                                                                                cTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCgt                                                                  >  2:544565/1‑150 (MQ=255)
                                                                                cTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCACCACTGCTTGCAGCAAATCGGAGGGACGGCGTTCgt                                                                  >  2:315957/1‑150 (MQ=255)
                                                                                      cGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCGTGTGTGa                                                            <  2:452995/150‑1 (MQ=255)
                                                                                        tggtgCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCGTGTGTGAgg                                                          <  2:1131217/150‑1 (MQ=255)
                                                                                               ggATAACGGGCCCGAAGCGCTATTCCGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCGTGTGTGAGGCGCAGTg                                                   <  1:444553/150‑1 (MQ=255)
                                                                                                   aaCGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCGTGTGTGAGGCGCAGTGGGcc                                               <  1:222664/150‑1 (MQ=255)
                                                                                                          ccGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCGTGTGTGAGGCGCAGTggg                                                 <  1:690368/141‑1 (MQ=255)
                                                                                                          ccGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCGTGTGTGAGGCGCAGTggg                                                 >  2:690368/1‑141 (MQ=255)
                                                                                                          ccGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCGTGTGTGAGGCGCAGTGGGCCGAACAGc                                        <  2:975743/150‑1 (MQ=255)
                                                                                                                       cggggATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCTGTCGTGTGTGAGGCGCAGTGGGCCGAACAGCTCATTGCGCTTgg                           <  2:525688/148‑1 (MQ=255)
                                                                                                                           gATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCGTGTGTGAGGCGCAGTGGGCCGAACAGCTCATTGCGCTTGGCttt                       <  1:315957/150‑1 (MQ=255)
                                                                                                                                        tGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCGTGTGTGAGGCGCAGTGGGcc                                               <  2:291531/113‑1 (MQ=255)
                                                                                                                                        tGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCGTGTGTGAGGCGCAGTGGGcc                                               >  1:291531/1‑113 (MQ=255)
                                                                                                                                              ccGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCGTGTGTGAGGCGCAGTGGGCCGAACAGCTCATTGCGCTTGGCTTTGCCCGGCAGGCGCCAGCTg    >  1:1023137/1‑150 (MQ=255)
                                                                                                                                                gAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCGTGTGTGAGGCGCAGTGGGCCGAACAGCTCATTGCGCTTGGCTTTGCCCGGCAGGCGCCAGCTGGc  >  2:1182462/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCGCAGCGTTACACCTACCGGCAACTGCTCAACGAAGTCGAATCCCGTGCGAGCTGCTTGCGCAGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCGTGTGTGAGGCGCAGTGGGCCGAACAGCTCATTGCGCTTGGCTTTGCCCGGCAGGCGCCAGCTGGC  >  NC_002947/3099349‑3099642

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: