Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,611,932 G→C 10.2% A63G (GCC→GGC)  PP_1412 ← 3‑hydroxyisobutyryl‑CoA hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,611,9320GC10.2% 108.9 / 3.2 39A63G (GCC→GGC) PP_14123‑hydroxyisobutyryl‑CoA hydrolase
Reads supporting (aligned to +/- strand):  ref base G (17/18);  new base C (2/2);  total (19/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.77e-01

GCACAACAACGGTTTGGGGTAGGTGTGCAGGGCAAAGTCCAGGCGGTATTCGGCGGCGAAGAAGGAGGCCGCCAGCGGTGGGACACTGCCGGGGTGGTCGC‑GGCAGGCTTGCGCCAGCGCACGGACATCGCCGCCTGCGCAAAAGGCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGACGATACCCGGGTCGCGGGCCCAGGCATGCAGTTGTTCGCCCAGCACCTCGATCATCGGCAGGTTCAGGGCATTGAGCGCCTTGGGCGCATCCAGGGTGGCGATGCCGATGCGG  >  NC_002947/1611787‑1612078
                                                                                                                                                  |                                                                                                                                                  
gCACAACAACGGTTTGGGGTAGGTGTGCAGGGCAAAGTCCAGGCGGTATTCGGCGGCGAAGAAGGAGGCCGCCAGCGGTGGGACACTGCCGGGGTGGTCGC‑GGCAGGCTTGCGCCAGCGCACGGACATCGCCGCCTGCGCAAAAGGCttt                                                                                                                                                >  2:27898/1‑150 (MQ=255)
gCACAACAACGGTTTGGGGTAGGTGTGCAGGGCAAAGTCCAGGCGGTATTCGGCGGCGAAGAAGGAGGCCGCCAGCGGTGGGACACTGCCGGGGTGGTCGC‑GGCAGGCTTGCGCCAGCGCACGGACATCGCCGCCTGCGCAAAAGGCttt                                                                                                                                                <  1:136698/150‑1 (MQ=255)
   caacaaCGGTTTGGGGTAGGTGTGCAGGGCAAAGTCCAGGCGGTATTCGGCGGCGAAGAAGGAGGCCGCCAGCGGTGGGACACTGCCGGGGTGGTCGC‑GGCAGGCTTGCGCCAGCGCACGGACATCGCCGCCTGCGCAAAAGGCTTTTGc                                                                                                                                             >  1:1037465/1‑150 (MQ=255)
      caacGGTTTGGGGTAGGTGTGCAGGGCAAAGTCCAGGCGGTATTCGGCGGCGAAGAAGGAGGCCGCCAGCGGTGGGACACTGCCGGGGTGGTCGC‑GGCAGGCTTGCGCCAGCCCACGGACATCGCCGCCTGCGCAAAAGGCTTTTGcccc                                                                                                                                          <  1:977017/150‑1 (MQ=255)
       aacGGTTTGGGGTAGGTGTGCAGGGCAAAGTCCAGGCGGTATTAGGCGGCGAAGAAGGAGGCCGCCAGCGGTGGGACACTGCCGGGGTGGTCGC‑GGCAGGCTTGCGCCAGCGCACGGACATCGCCGCCTGCGCAAAAGGCTTTTGCCCCg                                                                                                                                         >  2:876226/1‑150 (MQ=255)
            tttGGGGTAGGTGTGCAGGGCAAAGTCCAGGCGGTATTCGGCGGCGAAGAAGGAGGCCGCCAGCGGTGGGACACTGCCGGGGTGGTCGC‑GGCAGGCTTGCGCCAGCGCACAGACATCGCCGCCTGCGCAAAAGGCTTTTGCCCCGTTGcc                                                                                                                                    <  2:1087544/150‑1 (MQ=255)
                        gtgCAGGGCAAAGTCCAGGCAGTATTCGGCGGCGAAGAAGGAGGCCGCCAGCGGTGGGACACTGCCGGGGTGGTCGC‑GGCAGGCTTGCGCCAGCGCACGGACATCGCCGCCTGCGCAAAAGGCTTTTGCCCCGTTGCCGCGCAGCAGCAc                                                                                                                        <  1:7445/150‑1 (MQ=255)
                             gggCCAAGTCCAGGCTGTAGTCGGCGGCGAAGAAGGAGGCCGCCAGCGGTGGGACACTGCCGGGGTGGTCGC‑GGCAGGCTTGCGCCAGCGCACGGACATCGCCGCCTGCGCAAAAGGCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGa                                                                                                                   <  2:973208/150‑1 (MQ=255)
                             gggCAAAGTCCAGGCGGTATTCGGCGGCGAAGAAGGAGGCCGCCAGCGGTGGGACACTGCCGGGGTGGTCGC‑GGCAGGCTTGCGCCAGCGCACGGACATCGCCGCCTGCGCAAAAGCCTTTTGcgcag                                                                                                                                         <  2:80316/128‑5 (MQ=255)
                             gggCAAAGTCCAGGCGGTATTCGGCGGCGAAGAAGGAGGCCGCCAGCGGTGGGACACTGCCGGGGTGGTCGC‑GGCAGGCTTGCGCCAGCGCACGGACATCGCCGCCTGCGCAAAAGCCTTTTGcgcag                                                                                                                                         >  1:80316/1‑124 (MQ=255)
                              ggCAAAGTCCAGGCGGTATTCGGCGGCGAAGAAGGAGGCCGCCAGCGGTGGGACACTGCCGGGGTGGTCGC‑GGCAGGCTTGCGCCAGCGCACGGACATCGCCGCCTGCGCAAAAGGCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGAc                                                                                                                  <  2:1049720/150‑1 (MQ=255)
                              ggCAAAGTCCAGGCGGTATTCGGCGGCGAAGAAGGAGGCCGCCAGCGGTGGGACACTGCCGGGGTGGTCGC‑GGCAGGCTTGCGCCAGCGCACGGACATCGCCGCCTGCGCAAAAGGCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGAc                                                                                                                  <  2:747958/150‑1 (MQ=255)
                                       cAGGCGGTATTCGGCGGCGAAGAAGGAGGCCGCCAGCGGTGGGACACTGCCGAGGTGAACAAATGAAACCGCCTCCCCGCG‑ACGGACAGCGCCGCCTGCGCAAAAGGCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGACGATACCCgg                                                                                                         <  1:85533/150‑1 (MQ=255)
                                        aGGCGGTATTCGGCGGCGAAGAAGGAGGCCGCCAGCGGTGGGACACTGCCGGGGTGGTCGC‑GGCAGGCTTGCGCCAGCGCACGGACATCGCCGCCTGCGCAAAAGGCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGACGATAc                                                                                                             <  2:831543/145‑1 (MQ=255)
                                        aGGCGGTATTCGGCGGCGAAGAAGGAGGCCGCCAGCGGTGGGACACTGCCGGGGTGGTCGC‑GGCAGGCTTGCGCCAGCGCACGGACATCGCCGCCTGCGCAAAAGGCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGACGATAc                                                                                                             >  1:831543/1‑145 (MQ=255)
                                         ggCGGTATTCGGCGGCGAAGAAGGAGGCCGCCAGCGGTGGGACACTGCCGGGGTGGTCGC‑GGCAGGCTTGCGCCAGCGCACGGACATCGCCGCCTGCGCAAAAGGCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGACGATACCCGGGt                                                                                                       >  1:353700/1‑150 (MQ=255)
                                                       gcgAAGAAGGAGGCCGCCAGCGGTGGGACACTGCCGGGGTGGTCAT‑GGCAAGCTTGCGCCAGCGCACGGACATCGCCGCCTGCGCAAAAGGCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGACGATACCCGGGTCGCGGGCCCAGGCa                                                                                         >  1:929914/1‑150 (MQ=255)
                                                             aagGAGGCCGCCAGCGGTGGGACACTGCCGGGGTGGTCGC‑GGCAGGCTTGCGCCAGCGCACGGACATCGCCGCCTGCGCAAAAGGCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGACGATACCCGGGTCGCGGGCCCAGGCATGCAgt                                                                                   <  1:475035/150‑1 (MQ=255)
                                                             aagGAGGCCGCCAGCGGTGGGACACTGCCGGGGTGGTCGC‑GGCAGGCTTGCGCCAGCGCACGGACATCGCCGCCTGCGCAAAAGGCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGACGATACCCGGGTCGCGGGCCCAGGCATGCAgt                                                                                   >  1:554114/1‑150 (MQ=255)
                                                               ggaggCCGCCAGCGGTGGGACACTGCCGGGGTGGTCGC‑GGCAGGCTTGCGCCAGCGCACGGACATCGCCGCCTGCGCAAAAGGCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGACGATACCCGGGTCGCGGGCCCAGGCATGCAgttg                                                                                 >  2:1166180/1‑150 (MQ=255)
                                                                  ggCCGCCAGCGGTGGGACACTGCCGGGGTGGTCGC‑GGCAGGCTTGCGCCAGCGCACGGACATCGCCGCCTGCGCAAAAGGCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGACGATACCCGGGTCGCGGGCCCAGGCATGCAGTTGTTc                                                                              <  1:1166180/150‑1 (MQ=255)
                                                                  ggCCGCCAGCGGTGGGACACTGCCGGGGTGGTCGC‑GGCAGGCTTGCGCCAGCGCACGGACATCGCCGCCTGCGCAAAAGGCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGACGATACCCGGGTCGCGGGCCCAGGCATGCAGTTGTTc                                                                              <  2:1134732/150‑1 (MQ=255)
                                                                         aGCGGTGGGACACTGCCGGGGTGGTCGC‑GGCAGGCTTGCGCCAGCGCACGGACATCGCCGCCTGCGCAAAAGGCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGACGATACCCGGGTCGCGGGCCCAGGCATGCAGTTGTTCGCCCAGc                                                                       >  2:795385/1‑150 (MQ=255)
                                                                                ggACACTGCCGGGGTGGTCGC‑GGCAGGCTTGCGCCAGCGCACGGACATCGCCGCCTGCGCAAAAGGCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGACGATACCCGGGTCGCGGGCCCAGGCATGCAGTTGTTCGCCCAGCACCTCGa                                                                >  1:573216/1‑150 (MQ=255)
                                                                                  acacTGCAGGGGGGGTCGC‑GGCAGGCTTGCGCCAGAACACGGAGATCACCGGCTGCGCAACAGGCTTTTGACCCGTTGCAGCGCAGCAGCACTCAGACGATACCCGGGATGCGGGCCCAGGCATGCAGTTGTTCGCCCAGCACCTCGatc                                                              >  2:515604/1‑150 (MQ=255)
                                                                                        ccGGGGTGGTCGC‑GGCAGGCTTGCGCCAGCGCACGGACATCGCCGCCTGCGCAAAAGGCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGACGATACCCGGGTCGCGGGCCCAGGCATGCAGTTGTTCGCCCAGCACCTCGATCATCGGc                                                        >  2:1048500/1‑150 (MQ=255)
                                                                                                    c‑gGCAGGCTTGCGCCAGCGCACGGACATCGCCGCCTGCGCAAAAGGCTTTTGCCCCGTTGCCCCGCAGCAGCACGCAGACGATACCCGGGTCGCGGGCCCAGGCATGCAGTTGTTCGCCCAGCACCTCGATCATCGGCAGGTTCAGGGCa                                            <  1:44111/150‑1 (MQ=255)
                                                                                                                cgcCAGCGCACGGACATCGCCGCCTGCGCAAAAGGCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGACGATACCCGGGTCGCGGGCCCAGGCATGCAGTTGTTCGCCCAGCACCTCGATCATCGGCAGGTTCAGGGCATTGAGCGCCtt                                 >  1:788299/1‑150 (MQ=255)
                                                                                                                  ccAGCGCACGGACATCGCCGCCTGCGCAAAAGGCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGACGATACCCGGGTCGCGGGCCCAGGCATGCAGTTGTTCGCCCAGCACCTCGATCATCGGCAGGTTCAgg                                               <  2:1195946/134‑1 (MQ=255)
                                                                                                                  ccAGCGCACGGACATCGCCGCCTGCGCAAAAGGCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGACGATACCCGGGTCGCGGGCCCAGGCATGCAGTTGTTCGCCCAGCACCTCGATCATCGGCAGGTTCAgg                                               >  1:1195946/1‑134 (MQ=255)
                                                                                                                  ccAGCGCACGGACATCGCCGCCTGCGCAAAAGGCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGACGATACCCGGGTCGCGGGCCCAGGCATGCAGTTGTTCGCCCAGCACCTCGATCATCGGCAGGTTCAGGGCATTGAGCGCCTTgg                               <  2:613093/150‑1 (MQ=255)
                                                                                                                   cAGCGCACGGACATCGCCGCCTGCGCAAAAGGCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGACGATACCCGGGTCGCGGGCCCAGGCATGCAGTTGTTCGCCCAGCACCTCGATCATCGGCAGGTTCAGGGCATTGa                                        >  2:692203/1‑140 (MQ=255)
                                                                                                                   cAGCGCACGGACATCGCCGCCTGCGCAAAAGGCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGACGATACCCGGGTCGCGGGCCCAGGCATGCAGTTGTTCGCCCAGCACCTCGATCATCGGCAGGTTCAGGGCATTGa                                        <  1:692203/140‑1 (MQ=255)
                                                                                                                         aCGGACATCGCCGCCTGCGCAAAAGGCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGACGATACCCGGGTCGCGGGCCCAGGCATGCAGTTGTTCGCCCAGCACCTCGATCATCGGCAGGTTCAGGGCATTGAGCGCCTTGGGCGCATc                        >  2:356791/1‑150 (MQ=255)
                                                                                                                           ggACATCGCCGCCTGCGCAAAAGGCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGACGATACCCGGGTCGCGGGCCCAGGCATGCAGTTGTTCGCCCAGCACCTCGATCATCGGCAGGTTCAGGGCATTGAGCGCCTTGGGCGCATCCa                      >  2:184174/1‑150 (MQ=255)
                                                                                                                                 cgccgcCTGCGCAAAAGGCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGACGATACCCGGGTCGCGGGCCCAGGCATGCAGTTGTTCGCCCAGCACCTCGATCATCGGCAGGTTCAGGGCATTGAGCGCCTTGGGCGCATCCAGGGTgg                <  1:356791/150‑1 (MQ=255)
                                                                                                                                    cgccTGCGCAAAAGGCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGACGATACCCGGGTCGCGGGCCCAGGCATGCAGTTGTTCGCCCAGCACCTCGATCATCGGCAGGTTCAGGGCATTGAGCGCCTTGGGCGCATCCAGGGTGGcga             <  2:554114/150‑1 (MQ=255)
                                                                                                                                          ggcAAAAGCCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGACGATACCCGGGTCGCGGGCCCAGGCATGCAGTTGTTCGCCCAGCACCTCGATCATCGGCAGGTTCAGGGCATTGAGCGCCTTGGGCGCa                          <  2:159348/130‑1 (MQ=255)
                                                                                                                                          ggcAAAAGCCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGACGATACCCGGGTCGCGGGCCCAGGCATGCAGTTGTTCGCCCAGCACCTCGATCATCGGCAGGTTCAGGGCATTGAGCGCCTTGGGCGCa                          >  1:159348/2‑131 (MQ=255)
                                                                                                                                                  gCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGACGATACCCGGGTCGCGGGCCCAGGCATGCAGTTGTTCGCCCAGCACCTCGATCATCGGCAGGTTCAGGGCATTGAGCGCCTTGGGCGCATCCAGGGTGGCGATGCCGATGCgg  >  2:1081603/1‑147 (MQ=255)
                                                                                                                                                  gCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGACGATACCCGGGTCGCGGGCCCAGGCATGCAGTTGTTCGCCCAGCACCTCGATCATCGGCAGGTTCAGGGCATTGAGCGCCTTGGGCGCATCCAGGGTGGCGATGCCGATGCgg  <  1:1081603/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCACAACAACGGTTTGGGGTAGGTGTGCAGGGCAAAGTCCAGGCGGTATTCGGCGGCGAAGAAGGAGGCCGCCAGCGGTGGGACACTGCCGGGGTGGTCGC‑GGCAGGCTTGCGCCAGCGCACGGACATCGCCGCCTGCGCAAAAGGCTTTTGCCCCGTTGCCGCGCAGCAGCACGCAGACGATACCCGGGTCGCGGGCCCAGGCATGCAGTTGTTCGCCCAGCACCTCGATCATCGGCAGGTTCAGGGCATTGAGCGCCTTGGGCGCATCCAGGGTGGCGATGCCGATGCGG  >  NC_002947/1611787‑1612078

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