Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,172,360 G→A 11.3% G216G (GGC→GGT)  fadD‑I ← long‑chain‑fatty acid‑‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,172,3600GA11.3% 115.7 / 5.6 44G216G (GGC→GGT) fadD‑Ilong‑chain‑fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base G (21/18);  new base A (3/2);  total (24/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 = 9.94e-01

TGGTACAGCGGCAGCGGGGTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATCTTCTTCACGTACTTGATCACGC  >  NC_002947/5172223‑5172508
                                                                                                                                         |                                                                                                                                                    
tgGTACAGCGGCAGCGGGGTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTgc                                                                                                                                          <  1:977975/150‑1 (MQ=255)
                gggTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCAcgtcgt                                                                                                                          >  2:13712/1‑150 (MQ=255)
                          gATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCgg                                                                                                                <  1:1025301/150‑1 (MQ=255)
                              tCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTg                                                                                                            <  1:1143269/150‑1 (MQ=255)
                                gCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGc                                                                                                          <  1:1090557/150‑1 (MQ=255)
                                    ccTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGTAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCg                                                                                                      >  2:698110/1‑150 (MQ=255)
                                                   gAGCCCATCAGTGCACGGCACTACCGCACGTTCGCCAAAAGGTTGGAGAGCCTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTAGTTGGCTGGCGGGGTGGCTTCGGTCACGGGCTGGccc                                                                                       >  1:1114308/1‑150 (MQ=25)
                                                              tGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCtt                                                                                     >  1:1025505/1‑141 (MQ=255)
                                                              tGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCtt                                                                                     <  2:1025505/141‑1 (MQ=255)
                                                              tGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAgcgc                                                                            >  2:758323/1‑150 (MQ=255)
                                                                   ggCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGc                                                                                             >  2:900400/1‑128 (MQ=255)
                                                                   ggCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGc                                                                                             <  1:900400/128‑1 (MQ=255)
                                                                    gCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGc                                                                      <  1:862412/150‑1 (MQ=255)
                                                                      aCTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCgg                                                                                                     <  1:469300/117‑1 (MQ=255)
                                                                      aCTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCgg                                                                                                     >  2:469300/1‑117 (MQ=255)
                                                                        tGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCa                                                                  >  2:278226/1‑150 (MQ=255)
                                                                              aTGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGaaa                                                            <  1:758323/150‑1 (MQ=255)
                                                                                 ttGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGc                                                         >  2:272730/1‑150 (MQ=255)
                                                                                      caccaGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTg                                                                                                                                           <  2:95435/63‑1 (MQ=255)
                                                                                      caccaGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTg                                                                                                                                           >  1:95435/1‑63 (MQ=255)
                                                                                         caGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACg                                                      >  1:252001/1‑145 (MQ=255)
                                                                                         caGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACg                                                      <  2:252001/145‑1 (MQ=255)
                                                                                             ttGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCAAGTGGTTCGTTTGCTTTTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGc                                             <  2:947547/150‑1 (MQ=38)
                                                                                                       tCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGcc                                                                                  >  2:861018/1‑103 (MQ=255)
                                                                                                       tCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGcc                                                                                  <  1:861018/103‑1 (MQ=255)
                                                                                                       tCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAgg                                   >  1:456188/1‑150 (MQ=255)
                                                                                                           cATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCgg                               <  1:13712/150‑1 (MQ=255)
                                                                                                                   cTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCAtct                       >  1:431438/1‑150 (MQ=255)
                                                                                                                         gtacacCGGCGGTACCACCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGAACTTGCCCAGAGACAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATCTTCTTCa                 >  1:142665/3‑150 (MQ=255)
                                                                                                                         ccACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAg                                                                               <  1:772576/88‑1 (MQ=255)
                                                                                                                         ccACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAg                                                                               >  2:772576/1‑88 (MQ=255)
                                                                                                                          tacacCGGCGGTACCACCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATCTTCTTCAc                >  1:834635/2‑150 (MQ=255)
                                                                                                                          cacaCCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATCTTCTTCAc                >  2:641437/1‑150 (MQ=255)
                                                                                                                           acacCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATCTTCTTCACg               >  1:243707/1‑150 (MQ=255)
                                                                                                                           acacCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATCTTCTTCACg               >  1:809770/1‑150 (MQ=255)
                                                                                                                             acCGGCGGTACCACCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGc                                                         >  2:1159722/1‑106 (MQ=255)
                                                                                                                             acCGGCGGTACCACCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGc                                                         <  1:1159722/106‑1 (MQ=255)
                                                                                                                                 gtggtACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATCTTCTTCACGTACTTg         >  2:500478/1‑150 (MQ=255)
                                                                                                                                   ggtACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATCTTCTTCACGTACTTGAt       <  1:1080750/150‑1 (MQ=255)
                                                                                                                                    gtACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATcttctt                   <  1:1121441/137‑1 (MQ=255)
                                                                                                                                    gtACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATcttctt                   >  2:1121441/1‑137 (MQ=255)
                                                                                                                                    gtACCACCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATCTTCTTCACGTACTTGATc      <  2:834635/150‑1 (MQ=255)
                                                                                                                                     tACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTa                                     >  1:238597/1‑118 (MQ=255)
                                                                                                                                     tACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTa                                     <  2:238597/118‑1 (MQ=255)
                                                                                                                                        cgccgGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATCTTCTTCACGTACTTGATCACGc  <  1:698110/150‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                    
TGGTACAGCGGCAGCGGGGTGATGAGGATTTCGCAGCCTTCGTGCAGGTTGGAGCCCATCAGTGCACGGCACTGCAGCATGTTGGCCACCAGGTTGCGGTGGGTCAGCATGGCGCCTTTGGCCACACCGGTGGTACCGCCGGTGTACTGCAGCACCGCCACGTCGTTGGCCTGCGGGTTGGCTTCGGTCACGGGCTGGCCCTTGCCCAGCGCCAGGGCGTCATTGAAACGCACGGCCTGCGGCAGGTTGTAGGCCGGCACCATCTTCTTCACGTACTTGATCACGC  >  NC_002947/5172223‑5172508

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