Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,440,489 T→C 9.6% intergenic (+54/‑22) pcaF‑II → / → PP_2138 beta‑ketoadipyl‑CoA thiolase subunit beta/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,440,4890TC 9.6% 106.2 / 2.6 42intergenic (+54/‑22)pcaF‑II/PP_2138beta‑ketoadipyl‑CoA thiolase subunit beta/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (18/19);  new base C (2/2);  total (21/21)
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.00e+00

AGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAGTGTTGCGCGGCACTCCGAAACCGAAGAGTGGATGGTGTTCTACCAATGCCTGTATGGTGATTA  >  NC_002947/2440353‑2440632
                                                                                                                                        |                                                                                                                                               
aGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCGTAGTGCCGCGGTTTTGTTTTTTCAGGGTTTTTTCagtag                                                                                                                                     >  8:224056/1‑149 (MQ=255)
aGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGttt                                                                                                                                               <  6:305347/139‑1 (MQ=255)
aGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGttt                                                                                                                                               >  5:305347/1‑139 (MQ=255)
aGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCaggag                                                                                                                                     <  8:229689/149‑1 (MQ=255)
 gCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGAGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGCTTTTTTCaggagg                                                                                                                                    >  2:164445/1‑149 (MQ=255)
 gCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGAGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGCTTTTTTCaggagg                                                                                                                                    >  2:169609/1‑149 (MQ=255)
    aaaTGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGgtgt                                                                                                                                 <  1:277471/149‑1 (MQ=255)
            gTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATg                                                                                                                         >  8:103047/1‑149 (MQ=255)
                gCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATg                                                                                                                         >  6:281402/1‑145 (MQ=255)
                  tGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGACCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATa                                                                                                                   >  1:8123/1‑149 (MQ=255)
                           cGACCATGGGCGTCGGCCCGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGc                                                                                                          <  8:91740/149‑1 (MQ=255)
                            gACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACACGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCg                                                                                                         >  7:193061/1‑149 (MQ=255)
                               cATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGta                                                                                                       >  5:199670/1‑148 (MQ=255)
                                   tgCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCg                                                                                                  >  6:28252/1‑149 (MQ=255)
                                     cGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCTAAGGCTTCTGATAGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCgg                                                                                                 <  5:336404/148‑1 (MQ=255)
                                     agTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGCTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCgg                                                                                                 <  1:164445/147‑1 (MQ=255)
                                     agTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGCTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCgg                                                                                                 <  1:169609/147‑1 (MQ=255)
                                      gTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCa                                                                                               <  1:253652/149‑1 (MQ=255)
                                         ggCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTa                                                                                            <  5:281402/149‑1 (MQ=255)
                                               ggCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAAc                                                                                                                            <  5:143196/111‑1 (MQ=255)
                                               ggCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAAc                                                                                                                            >  6:143196/1‑111 (MQ=255)
                                               ggCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGACATGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACaaa                                                                                        <  1:240610/147‑1 (MQ=255)
                                                gCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAggg                                                                                     >  6:136000/1‑149 (MQ=255)
                                                  cAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGcc                                                                                   >  5:204363/1‑149 (MQ=255)
                                                        aTCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCACCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATAACGGCATTACAAATGGACTGAGTa                                                                             >  5:164109/1‑149 (MQ=255)
                                                                tGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACaaa                                                                                        <  4:304970/130‑1 (MQ=255)
                                                                tGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGAAGATACAACCAGGCGTATACCGGCATTACaaa                                                                                        >  3:304970/1‑130 (MQ=255)
                                                                      cGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAgt                                                                <  7:103047/148‑1 (MQ=255)
                                                                      cGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAgt                                                                <  8:193607/148‑1 (MQ=255)
                                                                      cGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAg                                                                 <  4:225476/147‑1 (MQ=255)
                                                                       gAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAgtgt                                                              >  8:192928/1‑149 (MQ=255)
                                                                        aaCGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAGTGtt                                                             >  4:244954/1‑149 (MQ=255)
                                                                        aaCACGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAGTGtt                                                             <  8:193061/149‑1 (MQ=255)
                                                                              tCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAGTGTTGCGCgg                                                       <  6:204363/149‑1 (MQ=255)
                                                                                  aTCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAGTGTTGCGCGGCACt                                                   >  8:224544/1‑149 (MQ=255)
                                                                                               aCAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAGTGTTGCGCGGCACTCCGAAACCGaa                                        >  2:174940/1‑147 (MQ=255)
                                                                                                     gcCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGGAACATGCAGATACAACCAGGCGTATACCGGCGTTACAAAGGGCCTGAGTACCGTGTCTTCAGTGTTGCGCGGCACTCCGAAACCGAAGAGTGGAt                                >  1:259364/1‑149 (MQ=255)
                                                                                                         gggCCTTGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAGTGTTGCGCGGCACTCCGAAACCGAAGAGTGGATGgtg                            <  7:192928/149‑1 (MQ=255)
                                                                                                                  gCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAGTGtt                                                             >  5:52649/1‑107 (MQ=255)
                                                                                                                  gCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAGTGtt                                                             <  6:52649/107‑1 (MQ=255)
                                                                                                                             gttttTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAGTGTTGCGCGGCACTCCGAAACCGAAGAGTGGATGGTGTTCTACCAATGCCTTTATg         <  2:8123/148‑1 (MQ=255)
                                                                                                                              ttttttCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAGTGTTGCGCGGCACTCCGAAACCGAAGAGTGGATGGTGTTCTACCAATGCCTGTATGGt       <  5:136000/149‑1 (MQ=255)
                                                                                                                                 tttCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAGTGTTGCGCGGCACTCCGAAACCGAAGAGTGGATGGTGTTCTACCAATGCCTGTATGGTGa     <  6:199670/148‑1 (MQ=255)
                                                                                                                                   tCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAGTGTTGCGCGGCACTCCGAAACCGAAGAGTGGATGGTGTTCTACCAATGCCTGTATGGTGATTa  >  3:284080/1‑149 (MQ=255)
                                                                                                                                        |                                                                                                                                               
AGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAGTGTTGCGCGGCACTCCGAAACCGAAGAGTGGATGGTGTTCTACCAATGCCTGTATGGTGATTA  >  NC_002947/2440353‑2440632

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: ‑