Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,324,596 C→T 7.1% A163V (GCG→GTG)  atoA → acetyl‑CoA:acetoacetyl‑CoA transferase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,324,5960CT 7.1% 177.2 / 3.5 56A163V (GCG→GTG) atoAacetyl‑CoA:acetoacetyl‑CoA transferase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (19/33);  new base T (2/2);  total (21/35)
Fisher's exact test for biased strand distribution p-value = 6.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

TGCCTGGGAAAATGGTGCCCGGTATGGGTGGCGCGATGGATCTGGTGACCGGGTCGCGCAAAGTGATCATCGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACG  >  NC_000913/2324449‑2324741
                                                                                                                                                   |                                                                                                                                                 
tGCCTGGGAAAATGGTGCCCGGTATGGGTGGCGCGATGGATCTGGTGACCGGGTCGCGCAAAGTGATCATCGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCg                                                                                                                                                  <  1:348017/149‑1 (MQ=255)
 gCCTGGGAAAATGGTGCCCGGTATGGGTGGCGCGATGGATCTGGTGACCGGGTCGCGCAAAGTGATCATCGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCgg                                                                                                                                                 <  3:318733/149‑1 (MQ=255)
  ccTGGGAAAATGGTGCCCGGTATGGGTGGCGCGATGGATCTGGTGACCGGGTCGCGCAAAGTGATCATCGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGt                                                                                                                                                <  4:177070/149‑1 (MQ=255)
    tGGGAAAATGGTGCCCGGTATGGGTGGCGCGATGGATCTGGTGACCGGGTCGCGCAAAGTGATCATCGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGc                                                                                                                                              <  1:189134/149‑1 (MQ=255)
    tGGGAAAATGGTGCCCGGTATGGGTGGCGCGATGGATCTGGTGACCGGGTCGCGCAAAGTGATCATCGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGc                                                                                                                                              >  1:318871/1‑149 (MQ=255)
          aaTGGTGCCCGGTATGGGTGGCGCGATGGATCTGGTGACCGGGTCGCGCAAAGTGATCATCGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGc                                                                                                                                        <  3:206730/149‑1 (MQ=255)
             ggTGCCCGGTATGGGTGGCGCGATGGATCTGGTGACCGGGTCGCGCAAAGTGATCATCGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTgg                                                                                                                                     <  1:492082/149‑1 (MQ=255)
             ggTGCCCGGTATGGGTGGCGCGATGGATCTGGTGACCGGGTCGCGCAAAGTGATCATCGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTgg                                                                                                                                     <  1:455355/149‑1 (MQ=255)
                gCCCGGTATGGGTGGCGCGATGGATCTGGTGACCGGGTCGCGCAAAGTGATCATCGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTa                                                                                                                                  <  4:522629/149‑1 (MQ=255)
                 cccGGTATGGGTGGCGCGATGGATCTGGCGACCGGGTCGCGCAAAGTGATCATCGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTAc                                                                                                                                 <  3:343944/149‑1 (MQ=255)
                      tATGGGTGGCGCGATGGATCTGGTGACCGGGTCGCGCAAAGTGATCATCGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAAc                                                                                                                            <  2:69053/149‑1 (MQ=255)
                           gTGGCGCGATGGATCTGGTGACCGGGTCGCGCAAAGTGATCATCGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCt                                                                                                                       <  3:420746/149‑1 (MQ=255)
                           gTGGCGCGATGGATCTGGTGACCGGGTCGCGCAAAGTGATCATCGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCt                                                                                                                       <  1:368180/149‑1 (MQ=255)
                                     ggATCTGGTGACCGGGTCGCGCAAAGTGATCATCGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGTGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGtt                                                                                                             >  2:402124/1‑149 (MQ=255)
                                     ggATCTGGTGACCGGGTCGCGCAAAGTGATCATCGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGTGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGtt                                                                                                             >  2:158408/1‑149 (MQ=255)
                                        tCTGGTGACCGGGTCGCGCAAAGTGATCATCGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTa                                                                                                          >  1:431654/1‑149 (MQ=255)
                                         cTGGTGACCGGGTCGCGCAAAGTGATCATCGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTAt                                                                                                         <  3:441021/149‑1 (MQ=255)
                                            gTGACCGGGTCGCGCAAAGTGATCATCGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCACAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGa                                                                                                      <  2:540652/149‑1 (MQ=255)
                                            gTGACCGGGGCGCGCTAAGTGATTAACGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGa                                                                                                      <  4:78041/149‑1 (MQ=255)
                                             tGACCGGGTCGCGCAAAGTGATCATCGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGAc                                                                                                     <  4:559265/149‑1 (MQ=255)
                                                         gcAAAGTGATCATCGCCATGGATCATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTGACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTgg                                                                                         <  3:347218/149‑1 (MQ=255)
                                                         gcAAAGTGATCATCGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTgg                                                                                         <  2:147234/149‑1 (MQ=255)
                                                         gcAAAGTGATCATCGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTgg                                                                                         >  4:22482/1‑149 (MQ=255)
                                                           aaaGTGATCATCGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCt                                                                                       >  1:544164/1‑149 (MQ=255)
                                                                    atcGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAAt                                                                              <  3:558148/149‑1 (MQ=255)
                                                                    atcGCCATCGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAAt                                                                              <  2:503742/149‑1 (MQ=255)
                                                                           tGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGAc                                                                       >  3:92186/1‑149 (MQ=255)
                                                                            ggAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACg                                                                      >  2:369931/1‑149 (MQ=255)
                                                                                   tGCGCCAAAGATGGTTCAACAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGa                                                               <  4:253213/149‑1 (MQ=255)
                                                                                       ccAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTa                                                           <  3:110870/149‑1 (MQ=255)
                                                                                         aaaGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGc                                                         >  2:409782/1‑149 (MQ=255)
                                                                                             aTGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCAcc                                                     >  2:475875/1‑149 (MQ=255)
                                                                                                 ttCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACcgtgc                                                 >  4:222247/1‑149 (MQ=255)
                                                                                                      cAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGcc                                            >  2:205608/1‑149 (MQ=255)
                                                                                                        aaaaTTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCaa                                          >  4:140363/1‑149 (MQ=255)
                                                                                                        aaaaTTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCaa                                          >  1:449325/1‑149 (MQ=255)
                                                                                                        aaaaTTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCaa                                          >  1:17416/1‑149 (MQ=255)
                                                                                                          acttttGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCaaaa                                        <  3:339031/147‑1 (MQ=255)
                                                                                                          aaTTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGTGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCaaaa                                        <  1:158408/149‑1 (MQ=255)
                                                                                                          aaTTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGTGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCaaaa                                        <  1:402124/149‑1 (MQ=255)
                                                                                                           aTTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCaaagc                                       <  3:309924/149‑3 (MQ=255)
                                                                                                                    cgcTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCg                              <  1:525675/149‑1 (MQ=255)
                                                                                                                      cTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGt                            >  3:152985/1‑149 (MQ=255)
                                                                                                                       tGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGtt                           <  3:190713/149‑1 (MQ=255)
                                                                                                                       tGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGtt                           <  4:529155/149‑1 (MQ=255)
                                                                                                                          tccATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTCTTGTCGACGGCATAATGTGGCACACCGAAATTGCCGACGGGTGTGATTCAGTCACCCTGCGTGCCAAAACAGAAGCTCGGTTTGa                        <  2:417562/148‑1 (MQ=255)
                                                                                                                           ccATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGaa                       <  3:439205/149‑1 (MQ=255)
                                                                                                                           ccATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGaa                       <  2:274015/149‑1 (MQ=255)
                                                                                                                               gCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTcg                   >  1:100066/1‑149 (MQ=255)
                                                                                                                                ccACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTcgc                  >  1:66442/1‑149 (MQ=255)
                                                                                                                                  actcactGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTAGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTcgcct                >  3:49951/1‑148 (MQ=255)
                                                                                                                                   ctcactGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTcgccgc               >  1:397869/1‑149 (MQ=255)
                                                                                                                                         gcgcAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCt         <  2:66442/149‑1 (MQ=255)
                                                                                                                                         gcgcAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCt         <  1:409782/149‑1 (MQ=255)
                                                                                                                                             aaCATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAAt     <  4:152985/149‑1 (MQ=255)
                                                                                                                                                aTGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACg  >  3:210912/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TGCCTGGGAAAATGGTGCCCGGTATGGGTGGCGCGATGGATCTGGTGACCGGGTCGCGCAAAGTGATCATCGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACG  >  NC_000913/2324449‑2324741

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑