Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,171,207 T→A 7.8% intergenic (+128/+55) PP_4548 → / ← fadD‑I oxidoreductase/long‑chain‑fatty acid‑‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,171,2070TA 7.8% 141.6 / 2.3 51intergenic (+128/+55)PP_4548/fadD‑Ioxidoreductase/long‑chain‑fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (23/24);  new base A (2/2);  total (25/26)
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.48e-01

GAAGCTGCGCTGAGTTGCGCCTACAGGGACTCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGAGCTCTTCATCACGCAGTTCGCGGCGCAGGATCTTGCCCACGTTGGTGGTCGGCAGCGAAT  >  NC_002947/5171067‑5171353
                                                                                                                                            |                                                                                                                                                  
gaAGCTGCGCTGAGTTGCGCCTACAGGGACTCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGccc                                                                                                                                            >  3:204211/1‑149 (MQ=255)
  aGCTGCGCTGAGTTGCGCCTACAGGGACTCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCTGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGACGCGGCCCtg                                                                                                                                          <  3:235065/149‑2 (MQ=255)
      gcgcTGAGATGCGCCTACAGGGACTCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCgg                                                                                                                                               >  8:99637/1‑140 (MQ=255)
      gcgcTGAGATGCGCCTACAGGGACTCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCgg                                                                                                                                               <  7:99637/140‑1 (MQ=255)
       cgcTGAGTTGCGCCTACAGGGACTCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCGTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTcg                                                                                                                                       <  2:63416/147‑1 (MQ=255)
       cgcTGAGTTGCGCCTACAGGGACTCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTcgc                                                                                                                                      <  1:343018/148‑1 (MQ=255)
              ttGCGCCTACAGGGACTCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCgg                                                                                                                                               >  4:137835/1‑132 (MQ=255)
              ttGCGCCTACAGGGACTCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCgg                                                                                                                                               <  3:137835/132‑1 (MQ=255)
              ttGCGCCTACAGGGACTCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGacac                                                                                                                                 <  2:222915/146‑1 (MQ=255)
              ttGCGCCTACAGGGACTCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGacac                                                                                                                                 >  1:222915/1‑146 (MQ=255)
                gcgcCTACAGGGACTCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCAGACGAAATCACACGGCCGCGTCGCGGCCCATTCGCGACACGAgga                                                                                                                            >  2:151962/1‑148 (MQ=255)
                          ggACTCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGc                                                                                                                                              >  1:71553/1‑121 (MQ=255)
                          ggACTCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGc                                                                                                                                              <  2:71553/121‑1 (MQ=255)
                              tCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACt                                                                                                                >  7:146658/1‑147 (MQ=255)
                              tCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGa                                                                                                              >  4:185319/1‑149 (MQ=255)
                                acaATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGc                                                                                                            <  8:22173/149‑1 (MQ=255)
                                acaATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGc                                                                                                            >  7:22173/1‑149 (MQ=255)
                                  aaTGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGt                                                                                                          >  8:315928/1‑149 (MQ=255)
                                      tCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCa                                                                                                       <  3:322057/148‑1 (MQ=255)
                                          aaGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCtt                                                                                                  >  3:134544/1‑149 (MQ=255)
                                          aaGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCG‑CCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTcgc                                                                                                                                      >  2:185088/1‑112 (MQ=255)
                                          aaGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCG‑CCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTcgc                                                                                                                                      <  1:185088/112‑1 (MQ=255)
                                                      cATCTTTTTGTTGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCGAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGGCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGa                                                                                      <  4:134544/149‑1 (MQ=255)
                                                            aTTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCACGCGAtcttc                                                                                 <  1:151962/148‑1 (MQ=255)
                                                                    gCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTtgtg                                                                                               >  4:242031/1‑126 (MQ=255)
                                                                    gCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTtgtg                                                                                               <  3:242031/126‑1 (MQ=255)
                                                                         gtgtCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGc                                                                     <  3:185319/147‑1 (MQ=255)
                                                                            tCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGcttctt                                                                >  1:170291/1‑149 (MQ=255)
                                                                                 aaaGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGAGc                                                            >  5:238574/1‑148 (MQ=255)
                                                                                        cgcAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCGAGGGCCGCGTCGCGGCCCTCTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTTGTCTGGCGATCTTCTTCAAGCCTTTCTTCTTGAGCTCTtcatc                                                    <  7:189161/149‑1 (MQ=255)
                                                                                        cgcAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGAGCTCTtcatc                                                    <  8:215542/149‑1 (MQ=255)
                                                                                                   cccAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTCAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGAGCTCTTCATCACGCAGTTcgc                                         >  5:165363/1‑149 (MQ=255)
                                                                                                    ccAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGAGCTCTTCATCACGCAGTTcgcg                                        <  5:72815/149‑1 (MQ=255)
                                                                                                      aGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTg                                                               <  5:226066/124‑1 (MQ=255)
                                                                                                      aGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTg                                                               >  6:226066/1‑124 (MQ=255)
                                                                                                            tGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTTCTTCTTGAGCTCTTCATCACGCAGTTCGCGGCGCAgg                                 <  6:195257/148‑1 (MQ=255)
                                                                                                               tGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGAGCTCTTCATCACGCAGTTCGCGGCGCAGGATc                              <  4:204211/148‑1 (MQ=255)
                                                                                                               tGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGAGCTCTTCATCACGCAGTTCGCGGCGCAGGATCt                             >  8:86026/1‑149 (MQ=255)
                                                                                                                 cATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGAGCTCTTCATCACGCAGTTCGCGGCGCAGGATCTTg                           <  2:170291/149‑1 (MQ=255)
                                                                                                                      ccgccgAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACGCGAGGCCGCCCCCACACTGa                                                                                                              >  4:217925/1‑61 (MQ=255)
                                                                                                                      ccgccgAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACGCGAGGCCGCCCCCACACTGa                                                                                                              <  3:217925/61‑1 (MQ=255)
                                                                                                                      ccgccgAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGAGCTCTTCATCACGCAGTTCGCGGCGCAGGATCTTGCCCAc                      >  7:69527/1‑149 (MQ=255)
                                                                                                                             aaTCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGAGCTCTTCATCACGCAGTTCGCGGCGCAGGATCTTGCCCACGTtggtg               <  8:69527/149‑1 (MQ=255)
                                                                                                                              gcgaaaGGGCCGCGACGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGAGCTCTTCATCACGCAGTTCGCGGCGCAGGATCTTGCCCACGTtggtgg              >  5:105221/6‑149 (MQ=255)
                                                                                                                              gcgaaaGGGCCGCGACGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGAGCTCTTCATCACGCAGTTCGCGGCGCAGGATCTTGCCCACGTtggtg               >  6:352976/6‑148 (MQ=255)
                                                                                                                              gcgaaaGGGCCGCGACGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGAGCTCTTCATCACGCAGTTCGCGGCGCAGGATCTTGCCCACGTtggtg               <  5:352976/143‑1 (MQ=255)
                                                                                                                                cccAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGAGCTCTTCATCACGCAGTTCGCGGCGCAGGATCTTGCCCACGTTGGTGGTc            <  7:315928/149‑1 (MQ=255)
                                                                                                                                      gCCGCGTCGCGGCCCTTTCGCGACACCAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGAGCTCTTCATGACGCAGTTCGCGGCGCAGGATCTTGCCCGCGTTGGTGGTCGGCAGc      >  6:195166/1‑149 (MQ=255)
                                                                                                                                        cgcgTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGAGCTCTTCATCACGCAGTTCGCGGCGCAGGATCTTGCCCACGTTGGTGGTCGGCAGCg     >  6:25441/1‑148 (MQ=255)
                                                                                                                                         gcgTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGAGCTCTTCATCACGCAGTTCGCGGCGCAGGATCTTGCCCACGTTGGTGGTCGGCAGCGaa   <  5:25441/149‑1 (MQ=255)
                                                                                                                                          cgTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGAGCTCTTCATCACGCAGTTCGCGGCGCAGGATCTTGCCCACGTTGGTGGACGGCAGCGAAt  >  4:285424/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
GAAGCTGCGCTGAGTTGCGCCTACAGGGACTCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGAGCTCTTCATCACGCAGTTCGCGGCGCAGGATCTTGCCCACGTTGGTGGTCGGCAGCGAAT  >  NC_002947/5171067‑5171353

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